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ABSTRACT 


This  thesis  is  concerned  with  the  application  of  cross-correlation 
techniques  to  linear  antenna  arrays.  The  basic  cross-correlation  system, 
which  is  considered,  consists  of  two  linear  receiving  arrays  excited  by  a 
distribution  of  remote  radio  sources.  The  terminal  voltage  of  each  array  is 
passed  through  a  narrow-band  RF  filter  and  the  two  resulting  signals  are  cross- 
correlated.  It  is  demonstrated  that  this  system  can  measure  the  mutual 
coherence  function  of  the  source  distribution.  To  do  this  the  patterns  of 
the  two  antennas  must  be  scanned  independently.  A  Fourier  analysis  shows 
that  the  cross-correlation  system’s  output  is  a  filtered  version  of  the  mutual 
coherence  function.  From  this  output  a  three-dimensional  principal  solution 
can  be  deduced;  it  is  a  generalization  of  the  one-dimensional  principal 
solution  given  by  Bracewell  and  Roberts  in  connection  with  radio  astronomy. 

In  addition  to  cross-correlating  the  voltages  obtained  from  two  distinct 
arrays  one  can  perform  a  cross-correlation  of  signals  obtained  from  one  of  the 
arrays.  By  dividing  the  signal  from  each  element  into  two  parts  (with 
predetermined  weighting) ,  and  by  combining  additively  each  set  of  signal s,  one 
obtains  two  output  voltages  from  what  are  effectively  two  coincident  arrays. 

These  signals  are  then  cross-correlated  in  the  usual  fashion.  A  similar  cross- 
correlation  can  be  performed  on  the  voltages  obtained  from  the  other  array  of  the 
system.  Finally,  a  fourth  cross-correlation  function  results  when  the  signals  of 
the  two  arrays  are  combined  in  reverse  order  and  it  is  the  complex  conjugate  of 
the  original  cross-correlation  function.  These  four  distinct  outputs  are  the 
elements  of  the  2x2  correlation  matrix  of  the  system.  \  Fourier  analysis  of 
the  correlation  matrix  leads  to  a  more  general  principal  solution  for  the  system 
as  a  whole  which  yields  considerably  more  information  about  the  source 


ACKNOWLEDGMENT 


It  has  been  my  privilege  to  prepare  this  thesis  and  technical  report 
under  the  direction  of  Professor  George  A.  Deschamps.  It  is  with  a  deep 
sense  of  gratitude  that  I  take  this  opportunity  to  thank  him  for  his  sound 
advice  and  helpful  criticism  during  the  course  of  this  work,  I  also  want 
to  thank  Professor  Y.  T.  Lo  for  his  continued  interest  and  many  suggestions. 
During  my  stay  at  the  Antenna  Laboratory  I  was  able  to  discuss  my  work  with 
a  great  many  individuals  and  I  appreciate  very  much  their  comments  and 
criticism, 

I  also  would  like  to  express  my  appreciation  to  Miss  Nanita  Patterson, 
who  typed  my  numerous  rough  drafts,  and  to  Mrs.  Arlene  Albert  and  her  staff 
in  the  publications  office.  Finally  I  want  to  thank  the  sponsor  of  my 
work.  Aeronautical  Systems  Division,  Wright-Patterson  Air  Force  Base,  Ohio, 
for  their  continued  support. 


iv 

CONTENTS 

Page 

1.  Introduction  1 

2.  Linear  Array  Theory  7 

3.  Description  of  the  Source  Distribution  13 

3.1  Assumptions  13 

3.2  Coherence  Theory  15 

3.3  Fourier  Analysis  17 

3.3.1  Time  Domain  17 

3.3.2  Space  Domain  18 

3.3.3  Combined  Time  and  Space  Domain  20 

3.4  Special  Limiting  Cases  21 

3.4.1  The  Coherent  Limit  21 

3.4.2  The  Incoherent  Limit  22 

4.  The  General  Theory  of  Cross-Correlation  25 

4.1  The  Complex  Cross-Correlation  Function  26 

4.2  The  Cross-Correlator  28 

4.3  The  Autocorrelation  Functions  34 

4.4  The  Correlation  Matrix  35 

4.5  The  Effect  of  Finite  Averaging  Time  37 

5.  The  Two-Antenna  Cross-Correlation  System  44 

5.1  Description  of  the  Antenna  Voltages  44 

5.2  Cross-Correlation  of  the  Antenna  Voltages  46 

5.3  Fourier  Analysis  of  the  System  Output  48 

5.4  The  Principal  Solution  50 

5.5  Direct  Measurement  of  the  Principal  Solution  54 

5.6  Cross-Correlation  System  Outputs  for  Coherent  and 

Incoherent  Source  Distributions  58 

5.7  Simplification  of  the  Cross-Correlation  System  for  the 

Case  of  Incoherent  Sources  61 

6.  The  Correlation  Matrix  of  the  Antenna  System  68 

6.1  Cross-Correlation  of  Signals  from  Two  Coincident  Antennas  68 

6.2  The  Correlation  Matrix  of  the  Two-Antenna  System  69 

6.3  Principal  Solution  70 

7.  Comparison  of  the  Cross-Correlation  and  the  Conventional  Antenna 

Systems  74 

7.1  A  Linear  Antenna  with  Square-Law  Detection  74 

7.2  Disadvantages  of  the  Conventional  System  76 

7.3  Comparison  of  the  Two  Systems  for  the  Case  of  Dolph-Chebyshev  77 

Array  Synthesis  * 


V 


CON  TENTS  (continued) 


Page 

8.  Multiple  Cross-Correlation  86 

8.1  A  Multiple  Cross-Correlation  Antenna  System  87 

8.2  Generalized  Correlation  Matrix  91 

8.3  Sources  with  Gaussian  Statistics  93 

9.  The  Application  of  Cross-Correlation  Antenna  Systems  to  Radar  96 

9.1  The  Cross-Correlation  Radar  System  96 

9.2  Optimum  Dolph-Chebyshev  Design  100 

9a3  Noise  Suppression  105 

9.4  Suppression  of  Jamming  Signals  109 

10.  The  Mapping  of  a  Target  Distribution  with  a  Cross-Correlation 

Radar  System  H0 

10.1  Description  of  the  Target  Distribution  HO 

10.2  System  Output  iLifi  the  Presence  of  the  Target  Distribution  HI 

10.3  Fourier  Analysis  of  the  System  Output  H4 

11.  Conclusions  H9 

References  120 

Appendix  A  The  Case  of  Two-Dimensional  Source  Distributions  and 

Planar  Antennas  124 

Appendix  B  Output  of  the  Multiple  Cross-Correlation  System  137 

Appendix  C  Correlation  Radar  System  142 

Appendix  D  Radar  Mapping  of  a  Target  Distribution  156 


vi 


ILLUSTRATIONS 

Figure 

Number  Page 

1.  The  Two-Antenna  Cross-Correlation  System*  ® 

2.  The  Cross-Correlator.  29 

3*  Spatial  Frequency  Plane  of  the  Cross-Correlation  Antenna 

System  with  the  "Aperture'1  Shown  Hatched.  52 

4.  Combined  Spatio-Temporal  Frequency  Space  for  the  Cross- 

Correlation  System  in  which  the  Region  is  Shown  as  a 

Rectangular  Parallelpiped ,  53 

5.  Impulse  Responses  of  Realizable  and  Unrealizable  Time-Domain 

Filters.  57 

6.  Aperture  Distributions  and  Spatial  Frequency  Spectrum  (for 

Incoherent  Sources)  of  Two  Uniformly  Weighted  Apertures*  63 

7.  Aperture  Dis tributions  and  Spatial  Frequency  Spectrum  of  the 

Compound  Interferometer  in  the  Presence  of  Incoherent  Sources.  65 

8.  Spatio-Temporal  Frequency  Plane  for  t}ie  Case  of  an  Incoherent 

Source  Distribution  with  the  Region  Shown  Hatched.  66 

9.  Spatial  Frequency_Plane  on  which  the  Cross-Correlation  System’s 

"Apertures"  QAB,  QBA,  and  QgB  are  Shown  Hatched.  71 

10.  A  Conventional  Linear  Antenna  which  Uses  Square-Law  Detection.  75 

11.  System  Used  to  Obtain  the  Optimum  Factor  Patterns  for  the  Dolph- 

Chebyshev  Product  Array.  79 

12.  Sidelobe  Level  Improvement  of  the  Dolph-Chebyshev  Cross-Correlation 

Array  Over  the  Conventional  Array  when  Both  Have  the  Same  Beam- 
widths  Measured  to  the  First  Null.  82 

13.  System  for  Cross-Correlating  the  Voltages  of  N  Antennas,  88 

14.  Cross-Correlation  Radar  System.  97 

15.  The  Aperture  Weighting  Network  Associated  with  the  Four  Coincident  102 

Radar  Arrays . 

16.  Sidelobe  Level  Improvement  of  the  Mattingly  and  Cross-Correlation 

Radar  Patterns  as  a  Fun^cion  of  the  Conventional  Dolph-Chebyshev  ^ 

Sidelobe  Level.  ^ 


1.  INTRODUCTION 


The  concept  of  correlation  has  been  used  to  advantage  for  many  years  in 

various  branches  of  mathematical  science.  It  is  basically  a  statistical  concept 

1  2 

and  such  well-known  statisticians  as  Karl  Pearson  and  A.  A.  Chuprov  have 

employed  it  extensively  in  their  work.  For  example,  Pearson  defined  the 

correlation  coefficient  of  two  sets  of  numbers  x  ,x  . .  x  .  x  . 

-nJ  -n+1**  9  n-1 J  n* 


which  is  the  cosine  of  the  angle  in  2n+l-space  between  the  two  vectors  X  and  Y. 

If  the  elements  of  X  and  Y  are  samples  of  a  time  sequence,  it  is  of 
interest  to  obtain  a  correlation  coefficient  between  the  two  sequences 
when  one  of  them  is  shifted  in  time,  time  being  represented  by  the 
index  of  the  correlated  data.  Hence  we  can  write 

,  n 

T*  lim  1  y' 

k  ~  n-^oo  2n+l  ^ • 

If  the  values  x  and  y  are  complex,  we  modify  the  above  expression  as 
k  k 

follows 

-  1  y  * 

X'k  ~  n-^co  2n+l  a  “l4 

*  =~n 
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4t 

where  the  superscript  indicates  the  complex  conjugate  of  the  quantity. 

The  above  expression^  as  a  function  of  the  index  k,  is  the  cross-correlation 
between  the  two  time  series  X  and  Y  It  is  a  simple  step  to  the  generalization 
of  these  discrete  time  series  to  continuous  ones  and  the  resulting  cross¬ 
correlation  function  for  the  continuous  case  is 


R(T> 


lim 
T— *  oo 


x(t)  y*(t-T)  dt 


Where  T  is  the  time  delay  of  one  function  relative  to  the  other.  Norbert 
3  4 

Wiener  *  has  inade  a  thorough  study  of  stationary  time  series  and  has  shown 
that  the  methods  of  Fourier  analysis  can  be  carried  over  into  the  domain 
of  such  functions  in  spite  of  their  statistical  nature. 

In  the  case  of  antennas  the  time  series  represent  signals  which 
originate  from  remote  radio  sources  distributed  in  various  directions. 
Consequently  the  cross-correlation  function  of  a  pair  of  signals  should 
now  be  written  as 


T(T,u,v) 


lim 

oo 


1_ 

2T 


V(t,u) 


V*(t-T,v) 


dt 


There  are  now  two  additional  variables,.  u  and  v,  which  represent  the 
directions  from  which  the  signals  originate. 

There  is  a  considerable  literature  on  the  subject  of  such  spatio- 
temporal  cross-correlation  functions  or  mutual  coherence  functions  as 
they  are  usually  called,  most  of  it  having  originated  in  optical  studies. 
The  mutual  coherence  function  of  two  light  disturbances  is  analogous  to 
the  cross-correlation  function  of  two  voltages..  Complete  correlation  and 
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zero  correlation  correspond  to  coherence  and  incoherence  respectively. 

More  generally  one  speaks  of  partial  coherence,  historically,  the  theory 

5  6 

of  partial  coherence,  according  to  Wolf  ,  dates  hack  to  Verdet  .  More 

7  89  10  11  12 

recent  contributors  have  been  Von  Laue  ,  Van  Cittert  ,  Zernike  ,  Wolf  3  3  , 

13  14  15  16 

Blanc-Lapierre  and  Dumontet  ,  Beran  ,  and  Parrent  ’  .  Although  developed 

primarily  to  describe  optical  phenomena,  the  theory  of  partial  coherence  can 

also  be  used  to  describe  the  correlation  properties  of  radio  fields. 

It  is  the  role  of  the  receiving  antenna  to  act  as  a  device  which  responds 

to  the  fields  from  the  remote  radio  sources.  In  many  practical  cases  the 

antenna  consists  of  a  number  of  smaller  elementary  antennas,  which  are 

located  on  a  straight  line  in  space.  Only  such  composite  antennas,  called 

linear  arrays,  will  be  considered  in  this  thesis.  They  were  chosen  not  only 

because  they  are  widely  used  but  because  they  lend  themselves  readily  to 

mathematical  analysis.  In  this  connection  the  concept  of  the  linear 

17 

antenna  array  as  a  filter  of  spatial  frequencies  has  become  quite 
widespread  in  recent  years.  The  definitive  paper  on  the  subject,  by 
Bracewell  and  Roberts^,  appeared  in  1954  in  connection  with  radio 
astronomy.  A  remote  source  distribution  (e.g.,  the  extraterrestrial  radio 
sources),  can  be  analysed  as  a  function  of  direction  by  Fourier  techniques 
to  yield  a  spectrum  of  spatial  frequency  components.  One  can  show  that  a 
finite  antenna  acts  as  a  low-pass  spatial  frequency  filter.  In  the  time 
domain  there  is  usually  a  narrow-band  RF  filter  associated  with  the 
antenna.  Thus  the  antenna  can  be  thought  of  as  a  spatio-temporal  filter 
of  the  incoming  radio  signals. 

The  basic  cross-correlation  antenna  system  consists  of  two  antennas 
whose  terminal  voltages  are  cross-correlated.  There  are  several  techniques 
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19  20 

to  do  this  electrically  *  .  The  first  practical  example  of  such  a 

21 

system  appears  to  have  been  the  Millfs  Cross  which  consists  of  two  long 

arrays  intersecting  at  right  angles.  When  the  array  voltages  are  cross- 

correlated  with  T  s=  0,  the  power  pattern- of  the  system  is  the  product  of  the 

two  field  strength  patterns  of  the  arrays.  This  product  pattern  is 

practically  zero  except  in  the  region  where  the  fani-shaped  beams  of  the  two 

arrays  intersect  at  right  angles.  This  produces  a  pencil  beam  which  is  the 

same  as  that  of  the  field  strength  pattern  of  a  rectangular  planar  antenna 

whose  dimensions  are  given  by  the  lengths  of  the  two  arrays*  Another  well 

22 

known  example  of  this  type  of  system  is  the  Compound  Interferometer  in 
which  a  uniformly  weighted  linear  aperture  of  length  L  is  placed  adjacent 
to  a  simple  two-element  interferometer  of  the  same  length.  The  cross¬ 
correlation  output  (again  with  T  -  0)  yields  a  product  pattern  which  is 
the  same  as  the  field  strength  pattern  of  a  uniformly  weighted  linear  aperture 
of  length  4L.  Both  of  these  systems  are  used  in  radio  astronomy. 

Taking  inspiration  from  the  success  in  this  field,  many  workers  have 
extended  the  use  of  cross-correlation  techniques  to  radio  direction  finding, 

radar,  and  general  communication  antenna  systems.  Some  of  these 

23  24  25  26 

contributors  are:  Arsac  ,  Barber  9  Berman  and  Clay  ,  Drane  and 

27,28  29,30  31  32,33,34 

Parrent  9  9  Hanbury  Brown  and  Twiss  ;  ,  Linder  ,  MacPhie  9 

ije  qC  n^7 

Pedinoff  and  Ksienski  ,  Price  ,  Weisby  and  Tucker  j  White,  Ball  and 
38  39 

Deckett  ,  and  Young  .  Since  some  of  the  above  investigators  deal  with 
single  and  multiple  products  of  antenna  voltages  the  resulting  outputs  no 
longer  are  related  to  inputs  by  the  law  of  superposition.  Consequently,  the 
generic  terms  "nonlinear  antenna  systems"  or  "data  processing  antenna  systems" 
have  become  associated  with  these  various  schemes. 
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In  this  thesis  a  general  theory  of  mapping  the  mutual  coherence  function, 

T(T,u, v),  of  a  distribution  of  remote  radio  sources  is  presented  for  the 

first  time.  Previous  workers  considered  only  the  case  of  incoherent 
18 

sources  and  zero  time  shift  (T=0)  .  It  will  be  shown  that  T(t^u,v)  can 

be  measured  with  a  correlation  antenna  system  whose  two  linear  arrays  are 

scanned  independently.  Four  distinct  outputs  from  the  system  are  obtained. 

These  four  correlation  functions  are  the  elements  of  the  correlation  matrix 

of  the  system.  A  three-dimensional  Fourier  analysis  shows  that  each  of  the 

system  outputs  is  a  filtered  version  of  the  mutual  coherence  function  of 

the  source  distribution.  The  analysis  is  for  a  partially  coherent 

distribution,  but  the  limiting  cases  of  complete  coherence  and  incoherence 

are  considered  and  are  shown  to  give  rise  to  special  types  of  system  outputs. 

Consequently,  If  one  or  the  other  of  these  types  of  outputs  is  observed  in 

practice,  one  can  infer  that  the  source  distribution  is  completely  coherent 

or  Incoherent  aa  the  case  may  be.  The  problem  of  taking  a  finite  time- 

average  in  measuring  the  correlation  of  signals  is  also  considered. 

It  is  shown  that  when  one  compares  cross-correlation  with  conventional 

antenna  systems  the  comparison  should  be  made  between  the  cross-correlated 

output  of  the  former  and  the  square-law  detected  output  of  the  latter  system. 

Indeed  it  is  shown  that  the  square-law  detected  output  is  a  degenerate  case 

of  the  correlation  output.  An  example  of  the  improvement  obtainable  by 

using  the  cross-correlation  system  is  given  for  the  case  of  the  optimum 

40 

Dolph-Chebyshev  array  design 

The  generalization  of  a  two-antenna  correlation  system  to  one  containing 
N  antennas  is  analyzed.  It  is  shown  that  such  a  system  yields  information 
about  the  higher  order  moments  of  the  temporal  probability  distribution 
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of  the  sources-  The  methods  of  Fourier  analysis  used  in  ihe  single 
correlation  case  can  be  generalized  to  more  dimensions  in  the  multiple 
correlation  case. 

Finally,  a  new  type  of  radar  system  employing  cross-correlation  is 
described  and  analyzed  in  detail.  Just  as  a  cross-correlation  receiving 
system  employs  two  antennas,  thf*  radar  system  uses  two  antennas  for  both 
transmitting  and  receiving,  four  antennas  in  all.  It  is  shown  that  such 
a  system  has  distinct  advantages  over  conventional  systems  employing  a 
single  antenna.  An  improved  Dolph-Chebyshev  design  becomes  possible 
and  interference  noise  can  be  suppressed.  The  system  could  also  be 
used  to  determine  possible  correlation  between  the  returns  from  various 


targets . 
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2 .  LINEAR  ARRAY  THEORY 


A  linear  array  consists  of  a  number  of  elementary  antennas  which  are 
colinear,  i.e.,  are  located  on  a  straight  line  in  space.  If,  as  is  customary, 
the  array  elements  are  identical,  then  the  array*s  radiation  pattern  is 
factorable;  it  is  the  product  of  the  element  pattern  and  the  pattern  of  the 
array  when  its  elements  are  replaced  by  isotropic  antennas.  Consequently 
in  linear  array  theory  we  need  consider  only  arrays  of  isotropic  antennas 
knowing  that  for  any  specific  type  of  element,  e.g.,  a  half-wave  dipole, 
we  can  obtain  the  actual  array  pattern  by  multiplying  the  elemeht 
pattern  by  the  pattern  of  the  isotropic  array.  We  should  also  note 
that  the  polarization  response  of  the  array  is  the  same  as  that  of  its 
elements.  Thus  by  removing  the  element  factor  we  are  reducing  linear  antenna 
array  theory  from  a  vector  to  a  scalar  formulation. 

If  the  elements  of  an  array  of  finite  length  are  increased  in  number 
without  limit  and  the  element  spacing  approaches  zero,  we  obtain  a 
"continuous  array",  as  opposed  to  the  discrete  array  mentioned  above.  In 
what  follows  we  will  refer  to  both  continuous  and  discrete  arrays,  whose 
elements  are  located  on  a  straight  line,  simply  as  linear  antennas . 

Now  ns  is  shown  in  Figure  1,  two  linear  antennas  (with  isotropic  elements) 
are  located  on  the  x  axis  of  a  convenient  coordinate  system.  One  is  of  length  L 

A 

and  the  other,  located  a  distance  ^  to  the  right  of  the  first,  is  of  length  L  . 

It  is  well  known  that  except  for  a  constant  factor,  the  far  field  patterns  of  the 

antennas  operating  at  frequency  w  are  given  by 

o 

fV2 

A(U)  =  ~h~  a(x)eJUXdx 

j-l  n 

A 


(1) 


Figure  1.  The  Two-Antenna  Cross-Correlation  System. 
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and 


B(u>  =  m 


(2) 


where  u  =  B  sin  0, 

B  =  phase  constant 

0  =  angle  measured  from  the  normal  to  the  antenna  aperture, 
x  =  aperture  coordinate, 

a(x)  -  antenna  A*s  aperture  weighting  function, 
b(x)  =  antenna  B*s  aperture  weighting  function. 

Note  that  the  time-harmonic  far  fields  are  represented  by  complex  quantities 
according  to  the  usual  convention.  We  will  use  this  convention  throughout  the 
thesis  for  both  fields  and  voltages.  However,  in  cases  where  two  voltages  are 
being  multiplied  one  must  consider  the  product  to  be  the  product  of  two  real 
quantities.  The  real  part  of  the  product  of  the  corresponding  complex 
voltages  is  not  the  same  as  the  product  of  the  real  voltages.  Consequently, 
in  any  case  such  as  this,  the  notation  "Re.../1  will  be  inserted  to  describe 
the  operation  precisely. 

Since  a(x)  and  b(x)  are  identically  aero  far  |x-f>L  /2  ahd  Jx~*|>L./2  respectively, 

A'  B- 

it  can  be  seen  from  Equations  (1)  and  (2)  that  A(u)  is  the  Fourier  transform 
of  a(x),  while  B(u)  is  the  Fourier  transform  of  b(x>.  Consequently  we  can  use 
the  shifting  theorem  to  show  that  if  for  example,  a(x)  is  changed  to  a(x)  e  ^UsX, 
the  pattern  becomes 


A(u-u  ) 

s 


oo 

r 

*^-00 


r  /  -v  -JusX,  JUX  . 

[a(x)e  b  ]  e  dx 


(3) 
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The  pattern  can  be  shifted  or  scanned  by  introducing  a  progressive  phase  shift 

of  u  =  B  sin  6  radians  per  meter  across  the  aperture.  This  is  called 
s  s 

* 

electrical  scanning  and  will  be  the  only  type  considered  in  this  thesis. 

It  also  follows  that  the  inverse  Fourier  transforms  of  A(u)  and  B(v)  are  the 
weighting  functions 

a(x)  =  ~  f  A(u)  e“jUX  du  (4) 

J-co 


and 


b(y) 


(5) 


It  will  be  observed  that  these  antennas  are  essentially  "one-dimensional" 
since  they  are  described  by  one-dimensional  aperture  and  pattern  functions. 

In  physical  (three-dimensional)  space  the  angle  0  =  constant  describes  a  cone 
whose  apex  is  at  the  origin  and  whose  axis  is  the  x  axis.  When  the  antenna 
is  receiving,  all  incident  signals  whose  propagation  vectors  are  parallel 
to  any  generator  of  this  cone  are  indistinguishable  by  the  antenna.  Their 
phasor  sum,  as  "seen"  by  the  antenna,  is  taken  as  the  field  of  "a  source" 


'  Tt  should  be  noted  that  in  practice,  electrical  scanning  is  complicated  by 
a  number  of  factors*  For  example,  suppose  the  antenna  is  matched  to  the  load 
when  its  pattern  is  beamed  in  the  broadside  direction  (u  =s  0) ;  then  in  order 
to  scan  the  pattern  to  the  u  =  v  direction,  a  progressive  phase  shift,  which 
of  course  changes  the  aperture*s  weighting  function,  must  be  introduced.  This 
will,  in  general,  result  in  a  mismatch  between  the  antenna  and  the  load.  The 
mutual  impedance  effects  of  the  antenna  elements  also  complicate  the  analysis 
of  electrical  scanning.  Thus  in  Equation  (3)  we  have  hopefully  assumed  that 
these  impedance  problems  have  been  solved  and  we  can  perform  an  ideal  scan  of 
the  pattern  in  the  u  domain.  However,  when  considered  as  a  function  of  angle 
0  the  pattern  becomes  distorted  when  scanned. 
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In  the  6  -  constant  direction.  In  addition,  by  thinking  of  S  as  the  propagation 
vector  "3,  which  indicates  the  direction  as  well  as  the  frequency  of  the  signals, 

ii 

we  see  that  u  =  B  sin  ©  is  actually  the  three-dimensional  scalar  product  of  B 

A  * 

with  a  unit  vector  x  parallel  to  the  antenna  axis. 

•*  A 

u  -  B  -  x  =  B  cos  (ir/2~0) 


or 


u  —  B  sin  9  *  (6) 

The  two-dimensional  counterparts  of  linear  antennas  are  planar  antennas  whose 
elements  are  located^  for  example,  in  the  z  ~  0  plane  of  a  cartesian  coordinate 
system.  The  one-dimensional  Fourier  analysis  of  linear  antennas  can  be  extended 
to  these  two-dimensional  antennas  with  the  transform  coordinates  being  given  by 
x,  y,  and  u  =  B  sin  0  cos  <p,  v  =  B  sin  ©  sin  (p .  In  this  case  the  propagation 
vector  can  be  represented  as  a  point  on  a  sphere  whose  coordinates  are  (Pj3}cp) 
and  by  analogy  with  the  one-dimensional  case  we  have 

u  =  B  °  "x  -  B  sin  0  cos  <P  } 

(7) 

v  =  B  •  y  =  B  sin  0  sin  <P  . 


i 

I 


* 

In  Appendix  A  is  a  more  detailed  account  of  the  relation  between  the  three- 
dimensional  physical  space  and  its  projected  counterparts  in  both  one  and  two 
dimensions . 


i 
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Admittedly,  such  antennas  are  of  much  greater  practical  value  than  linear 
antennas.  However  the  method  of  analysis  is  essentially  the  same  for  both 
and  for  simplicity  of  notation  we  will  restrict  ourselves  to  the  one¬ 
dimensional  case.  Generalization  to  the  planar  antenna  of  the  results 
obtained  is  quite  straightforward. 


13 


3 o  DESCRIPTION  OF  THE  SOURCE  DISTRIBUTION 

In  almost  all  cases  of  practical  interest  the  exact  nature  of  the 

source  distribution  is  unknown;  if  it  were  known  exactly  a  priori  it 

would  never  be  necessary  to  build  an  antenna  to  measure  it.  However^  we  can 

make  some  reasonable  assumptions  about  the  distribution  which  are  deduced 

from  some  of  its  known  physical  properties*  Furthermore^  the  statistical 

3  4 

description^ which  proved  so  useful  in  communication  theory  9  will  be  used 

here.  The  basic  concepts  of  this  jbype.  of  approaehuare  embodied  in  the 

5  1 5  16 

theory  of  partial  coherence  as  developed  by  Born  and  Wolf  and  Par rent  9 
301  Assumptions 

It  will  be  assumed  that  the  source  distribution  possesses  the 
following  physical  and  statistical  properties: 

1)  The  sources  are  remote  from  the  receiving  antenna  and  fixed  in 
space  during  the  observation. 

2)  The  direction  of  a  source”  is  given  by  the  single  parameter 
u  =  P  sin  0, 

3)  In  the  time  domain  the  sources  are  stationary;  they  emit  random 

signals  whose  statistical  properties  are  invariant  under  a  shift 
42 

of  the  time  origin 

4)  The  sources  are  ergodic;  the  time  averages  of  various  quantities 

are  equivalent  to  their  mathematical  expectations  or  statistical 
42 

averages 

5)  The  signals  emitted  by  the  sources  are  described  by  zero-mean 
complex  random  variables  with  statistically  independent  real  and 


imaginary  parts. 
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6)  The  signals  from  sources  in  different  directions  are  not  necesr- 
sarily  statistically  independent. 

43 

7)  The  statistics  of  the  signals  are  gatissian 

8)  The  signals  are  quasi-monochromatic .  The 

bandwidth  2  a  w  is  much  less  than  the  carrier  frequency 


The 'above  properties  constitute  an., acceptably  .description  of  many 
sources  encountered  by  the  antenna  engineer.  For  example,  remote 
communication  transmitters,  radar  targets,  radio  stars,  and  even  earth 
satellites,  produce  signals  whose  general  description  is  given  above. 

Since  the  sources  are  remote,  the  field  due  to  each  is  incident  on 
-  the  antenna  system  in  the  form  of  a  plane  wave.  The  RF  field  at  time  t, 
at  the  origin  of  the  anteriha  system  (see  figure  1),  which  is  due  to  the 
source  in  the  u  direction,  can  be  written  as 

Jw  t 

e  (t,u)  =  V  (t,u)  e  u 

where  the  complex  scalar  function,  e  (t,u),  can  be  thought  of  as  one  of 
the  components  of  the  electric  field  vector  due  to  the  remote  source,  i«e., 
the  component  to  which  the  antenna  elements  respond.  V  (t,u)  is  the  com¬ 
plex  modulation  envelope  of  the  carrier  signal  at  frequency  oj  and  is  the 

o 

zero-mean  random  variable  with  independent  real  and  imaginary  parts  men¬ 
tioned  above.  This  signal  contains  the  desired  information  aDout  the 
source.  Accordingly,  lefcus  suppress  the  carrier  factor  in  Equation  (8)  and 
focus  our  attention  on  the  coherence  properties  of  V  (t,u). 
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3.2  Coherence  Theory 

5 

Following  Born  and  Wolf  ,  we  define  the  complex  degree  of  coherence 
between  the  fields  from  the  source  in  the  u  direction  and  the  source  in 
the  v  direction' as 


<  V  (t,u)  V*  (i:-T,v)>^ 


— 

^  \  W./  T 

■  n  7 _ 

<|v  (t,u)|2>  . 

<lv  (t,v)|2> 

1  '  1  N 

1  '  1 

(9) 


where  <  f(t)  .->  denotes  the  time-average  of  f(t),  i.e., 


<  f(t)  >  =  T 


lim 


oo 


f  (t)dt  o 


(10) 


One  signal  is  delayed  by  t  seconds  relative  to  the  other  before  their 

*  * 

Hermitian  product  is  formed  and  time-averaged.  The  normalized  form  of 
this  average  is  V  (t.u^v)  and  it  can  be  easily  shown  that 


0  < 


v  (T,u,v)  <:  1 


(ID 


"’’Both  u  and  v  are  independent  running  variables  in  the  {3  sin  0  domain,  i.e., 
we  could  write  u  -  (3  sin  9  and  v  =  f3  sin  0  .  They  should  not  be  confused 
with  the  orthogonal  coordinates  u  and  v  in  ¥he  two-dimensional  antenna 
case  . 

*  Jw  t 

More  specifically,  Born  and  Wolf  use  the  total  complex  phasor  V(t,u)e*  9 
rather  than  just  the  envelope  V(t^u).  Their  complex  degree  of  coherence 
differs  from  Equation  (9)  by  an  RF  phase  factor  * 
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If  \y  (TjUjV)!  =  0,  for  u  *  v,  the  two  sources  are  completely 
incoherent  while  if  \y  (t,  u,  v)i  -  1  they  are  completely  coherent. 

Note  that  y  (t^u^v)  is  a  function  only  of  the  difference  of  the 

* 

time  arguments  of  the  two  field  phasors  V  (t,  u)  and  V  (t-T, v).  This 
is  a  consequence  of  the  stationary  nature  of  the  sources.  Let  us 
define 

T  (t,u,v)  b  <  V  ( t,  u)  V*  (t-r,v)  >  (12) 

* 

to  be  the  mutual  coherence  function  of  the  source  distribution  . 

For  r  a  0  and  u  «  v,  this  function  is  real  and  is  proportional  to  the 
average  power  radiated  by  the  source  in  the  u  direction  toward  the 
antenna,  i.e.,  it  is  a  measure  of  the  temperature  brightness  of  the 
source. 

In  practice  it  is  the  temperature  brightness  of  the  two-dimensional 
distribution  of  remote  sources  which  we  wish  to  measure.  Appendix  A 
is  devoted  to  a  description  of  a  planar  array  which  can  measure  such  a 
distribution.  The  effect  of  the  element  pattern  and  of  changing 
variables  from  the  (u,  v)  domain  to  the  (0,c{>)  domain  is  considered. 
Nevertheless,  such  an  analysis  is  simply  a  generalization  of  the  one¬ 
dimensional  case  which  for  simplicity  of  notation  we  are  considering  in 
the  main  body  of  this  thesis. 

Since  we  have  assumed  that  all  of  the  sources  are  remote,  they  must 

* 

Born  and  Wolf  define  the  mutual  coherence  function  as  T  (t^u^v)  = 

T  (t,u,v)  eJ  o  . 
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lie  in  the  so-called  visible  range  of  the  antennas.  Consequently  we 
can  write 


V  (t,u) 


for 


u|  >  P, 


and 


■y  (t,u.v)  =  0,  for  u  >  P  or  v  >  P, 


(13) 


T  (t^UjV)  3  0,  for  u  >  (3  or  u  >  P 


44  45  46 

In  the  terminology  of  distribution  theory  *  ’  we  say  tl^at  "yCr,  u, v)  and 

T  (Tj  u,v)  have  a  finite  support  in  the'  uv  domain  while  V  (t,u)’s 
support!®  finite  in  the  u  domain. 

3.3  Fourier  Analysis 

3.3.1  Time  Domain 

Since  V  (t,  u)  is  a  sample  function  of  a  stationary  random  process  ith- 
inverse  Fourier  transform  in  the  time  domain  does  not  exist,  i.e.,  V  (t,  u) 
is  not  absolute-square  integrable,* 


T 


(14) 


However, the  complex  mutual  coherence  function,  T  (t,  u,  v),  which  may  be 


written  as 
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T  (T,u,v) 


lim  _L  /Ty  V*(t-T,v)dt 

T — ►  °o  2T  -T 


(15) 


does  exist.  This  means  that  there  is  a  finite  average  cross-power  be¬ 
tween  the  signals  from  the  u  and  v  directions,  as  observed  :*t  the  janic-iuia 
system's  origin. 

Thus  we  define  the  cross-power  spectral  density  or  power  spectrum 
to  bp  th^  inverse  Courier  -translatra;:  of  T  (T-uyv)^  :  . 


t  (^,u,  v)  = 


/ 

-  oo 


T  (T,u}v)e  ^WdT; 


It  also  follows  that 


T  (T,u,v) 


■  ,  ■  i(jJT 

u,  v)e  dto 


(16) 


(17) 


3.3.2  Space  Domain 

i 

Because  V  (t,u),  as  a  function  of  u,  has  a  finite  support 
(V  (t, u)  =  0  for  | u  |  >p  )  and  is  square  integrable,  it  is  possible  to 
define  its  inverse  Fourier  transform  in  the  space  or  spatial  fr<-.;  -mcy 
domain.  Thus  we  let  its  spatial  frequency  spectrum  at  time  t  be 
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oO 

v  (t^x)  -  — J  V(tju)e'^uxdu 
JljF  -  OO 

p 

*—  f  V(t,u)e-1UXdu  .  (18) 

mjp 

The  designation  "spatial  frequency  spectrum'1  is  used  as  a  term  which  is 

analogous  to  "temporal  frequency  spectrum"  in  the  time  domain* 

Physically,  v  (t9y)  is  the  complex  envelope  of  the  field  which  exists  at  a 

point  on  the  x  axis  which  is  y  units  from  the  origin  (see  Figure  1)*  This 

field  is  the  integral  of  all  the  plane  waves  from  the  remote  sources 

i  uy 

modified  by  the  phase  factor  e  On  the  other  hand  we  can  recover  the 

complex  field  due  to  the  source  in  the  u  direction  as  follows 


V 


(t,  U)  ==  — 

\TiF 


r  v(tjX)e  "*U“dx 


-  OO 


(19) 


Let  us  define 


t  (T,x,y) 


* 

V 


(20) 


to  be  the  mutual  coherence  function  between  the  field  envelopes  at  the 
points  x  and  y4  If  we  substitute  Equation  (19)  into  the  above  expression, 


there  results 
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(w»  -  £ v  <t’u,eJUXdu//  v*(t'T,v> 


le"Jvydv>  (21) 


,P.P 


1 

27T 


f  f  <^V(t,  u)  V*  (t-T;  v-)^>  e 

J-p 


j (ux-vy) 

dudv 


(22) 


or 

*PrP  j  (ux-vy) 

=  .  .  1  II  T(T  U  v)e  dudv.  (23) 

t  (T,x,y)  =  —  J  J 

211  -P  P 

Thus  for  a  given  time  delay  t  the  spatial  frequency  s.pectrura»  : 
of  T(t?u,v)  is  given  by  the  mutual  coherence  function  t(T^x^y)  as  defined 
by  Equation  (23) . 

3.3,3  Combined  Time  and  Space  Domain 

Using  Equations  (16)  and  (23)  we  can  immediately  define  the  combined 
temporal  and  spatial  spectrum.. of  T (t^  U^'v)  as  '  ' 


....  .  1  fP  „<T 

t(-,x>y>  =  —  8/  |  /  /  T(T,u,v)e 

(270  V  J  Jn  J„ 


(24) 


-oo  -j3  -(3 


and  of  course  the  inverse  is 
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T(t 


»  P  B  j  (ur-ux+vy) 

,  u,v)  =  - v  fff  t(v>,xty)e  do*ixdy 

C2ff)  2  J  J  . 

-p  -p 


(25) 


-OO 


Note  that  the  argument  of  the  exponential  function  in  each  transform  is 
t  4ux  -vy).  The  negative  sign  associated  with  vy  occurs  because  the 
mutual  coherence  functions  are  expended  values  of  a  Hermitian  product, 
i.e,,  the  product  of  one  complex  function  and  the  complex  conjugate  of 
another.  This  taking  of  the  conjugate  gives  rise  to  the  negative  sign  as 
can  be  seen  in  Equations  (20)  to  (23) . 

3.4  Special  Limiting  Cases 

3,4.1  The  Coherent  Limit 

The  limit  of  complete  coherence  between  sources  in  the  u  and  v  direc- 
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tions  has  been  defined  by  Par rent  to  mean  that 

[  v  (T,u,v)  j  =  1  (26) 

for  all  values  of  t.  However,  Parrent  has  shown  that  only  a  strictly 
monochromatic  source  distribution  will  satisfy  Equation  (26);  iTfte  mutual 
coherence  function  of  such  a  distribution  can  be  written  as 

T  (t,  u,  v)  =  V  (u)  V  (v) 


(27) 
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Where  V  (u)  is  a  fixed  amplitude  and  phase  factor  associated  with  the 
source  in  the  u  direction, 

In  practice  one  can  at  best  deal  with  quasi-monochromatic  signals 
for  which  the  power  spectrum  t  (co^u^v)  is  practically  zero  except  for 
|toj<  a  oj  where  A  w  «<  w  .  In  such  a  situation  one  defines  a  coherent 
quasi-monochromatic  source  distributionl5to  be  one  for  which 


*  '  12  ff 

T  (t^  u,v)  =  V  (u)  V  (v),  for  T  «  — 


(28) 


2V 

The  mutual  coherence  function  is  independent  of  T  (for  T  «  and 

* 

separable  into  the  cartesian  product  of  two  functions  V(u)  and  V  (v) 
which  are  complex  conjugates.  Since  T  (t^u^v)  is  separable  for  small  T 
it  follows  that  the  spatial  frequency  spectrum  is  also  separable. 


=  *  2  IT 

t  (t,  x,y)  -  V  (x)  V  (y),  for  T  «  — 


(29) 


Thus  for  small  values  of  T}  the  coherent  quasi-monochromatic  distribution 
behaves  as  one  which  is  strictly  monocromatic. 

3.4.2  The  Incoherent  Limit 

Complete  incoherence  between  the  sources  in  the  u  and  v 
directions  is  characterized  by 

T  (TjUjV)  =  T  (T,u,u)  6  (u-v)  (30) 

where  6(u)  is  the  Dirac  delta. 
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The  mutual  coherence  function  is  non-zero  only  when  u=v,  in  which  case 
we  obtain  the  self-coherence  or  autocorrelation  function  of  the  source 
in  the  u  direction.  The  spatial  frequency  spectrum  of  the  incoherent 
distribution  is 

,  #8-6  j(ux-vy) 

t  (t  x  y)  =  ~r  lFfF  T  (T,  u,  u)  6  (u-v)e  dudv  (31) 

211  Jp 


or 


t  (Tj x-y) 


ju(x-y) 

T  (t,  u,  u)e  du  . 


(32) 


The  spatial  frequency  spectrum  for  Incoherent,  sources  is  a  function  only 
of  the  difference,  x-y,  of  the  spatial  frequency  coordinates.  If  we  let 
x-y  =  z  we  can  define  a  one-dimensional  spatial  frequency  spectrum  as 


*  <T,*>  -  r= 

4  27 r 


juz 

(t,  u)e  du 


(33) 


where 


T  (t,u) 


T  (t.u.u,) 

fi  23r 


(34) 


in  the  incoherent  Cdse.  In  addition,  one  can  use  Equation  (12)  to  show 
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that  T  (o,u,u)  and  hence  T  (o, u)  are  real  and  non-negative  functions  of 
u.  Physically,  this  means  that  the  average  energy  flow  from  the  remote 
sources  is  non-negative. 


25 


4.  THE  GENERAL  THEORY  OF  CROSS  -OCRRELAT  I  ON 

Before  turning  our  attention  to  the  specific  case  of  the  two-antenna 
cross-correlation  system,  we  will  consider  the  cross-correlation  of  two 
arbitrary  narrow-band  RF  signals.  We  will  show  that  a  complex  cross¬ 
correlation  function  is  associated  with  these  two  real  signals.  Then  we 
will  consider  in  some  detail  a  practical  RF  cross-correlator;  this  device 
can  be  used  to  measure  the  real  and  imaginary  parts  of  the  complex  cross  - 
correlation  function. 

When  the  two  signals  are  indentical  the  cross-correlation  becomes 
an  autocorrelation.  If  the  two  signals  are  combined  in  reverse  order, 
another  cross-correlation  results  which  is  the  complex  conjugate  of  the 
original.  The  two  cross-correlations  and  the  two  autocorrelations  cair 
be  arranged  as  the  elements  of  a  2  x  2  matrix,  the  correlation  matrix, 
which  completely  characterizes  the  correlation  properties  of  the  two 

^  At  -  e 

signals.  This  matrix  is  analogous  to  the  well  known  coherency  matrix^ 
used  in  the  study  of  the  polarization  of  light. 

Finally,  we  will  return  to  the  cross-correlator  to  consider  the 
practical  problem  of  taking  finite  time-averages  instead  of  the  infinite 
time  averages  which  are  formally  specified  in  the  definitions  of  the 
correlation  functions, 

*  43 

See  Davenport  and  Roof  page  60. 
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4.1  The  ComplexiCfross-Correlation  Function 

Let  us  consider  two  RF  voltages  whose  cross-correlation  function  is 
to  be  determined.  We  can  represent  them  by 

jw  j 

a(t)  =  Re  (35) 

J»  t 

b(t)  =  R3B(t)e  (36) 

where  A(t)  and  B(t)  are  slowly  varying  complex  modulation  envelopes.  The 
cross-correlation  of  the  two  real  voltages  is,  by  definition, 


.  ! 


R 


AB 


(T) 


a(t)  b(t-T)  dt  t 


(37) 


In  complex  notation  this  becomes 


_  • 


R 


AB 


(T) 


4' 


(2t-T) 


A(t)  B(t-T)e 


+ 


<^A(t) 


*  t 

B  (t-T)e 


(38) 


The  first  term  on  the  right  side  of  the  above  equation  is  the  double 
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rab(t)  - 


AB 


(T)e 


■JW  T 

o 


(43) 


or 


RA3(t)  =<^A(t)  B*(t-r)  > 


(44) 


This,  of  course^  is  a  shift  of  the  cross-correlation  spectrum  from 

T 

to  0  and  removes  the  RF  M carrier”  factor,  e  .  from  the  correlation 
o  *  ’ 

function . 

We  see  that  a  complex  cross-correlation  function  is  associated  with 
any  pair  of  real,  narrow-band,  RF  voltages.  In  the  following  section  we 
will  show  that  the  real  and  imaginary  parts  of  this  complex  function  can 
be  measured  by  a  system  employing  synchronous  detection. 

4.2  The  Cross-Correlator 

We  will  now  describe  a  device  which  will  cross-correlate  the  two 

22 

real  RF  voltages  a(t)  and  b(t).  In  practice  the  Compound  Interferometer 

employs  this  type  of  cross-correlator.  Now,  as  is  shown  in  Figure  2,  one 

of  the  signals,  b(t),  is  delayed  by  the  variable  amount  T:  before  it  is 

fed  into  a  frequency  shifter.  The  frequency  shifter  increases  the 

carrier  frequency  of  the  signal  from  u  to  w  +  u>  with  co  >»  w  ,  The 

o  o  1  o  1 

output  of  this  device  can  be  written  as 


30 


j[W  T  -(u  +  W  )t] 

bCt-T,^)  =  Be  B  (t-T)e  °  (45) 

This  signal  along  with  a(t),  enters  a  mixer  which  forms  the  square  of 
their  sum. 


[  a { t)  +  b(t-T;co^)  ]  =1 


Hi* 


2  2  -j[2CV  “2W,  t]  -J[«  T  -“it] 

(t)  +  B  (t-T)e  +  A(t)  B(t-T)e  ° 


]  r°l] 


► 

4 


[  j[w  T  •*  W  t]  1 

+  -Re^jA(t)|  +  | B (t-T)  |Z  +  A(t)  B*(t-T)e  °  1  >  . 

I  J 


(46) 


The  first  term  on  the  right  side  of  the  above  equation  is  the  double 
frequency  component  of  the  mixer  output.  If  a  band-pass  filter  with 
center  frequency  at  w  =  is  used  to  reject  this  signal  as  well  as  the 
signal  whose  spectrum  is  centered  at  w  =  we  obtain  as  the  filter  out¬ 
put 


VBP(t> 


I  Re 


*  Jtv  -v] 

A(t)  B  (t-T)e  °  l 


(47) 


This  is  possible  only  if  u  >  A  «  where  A  «  is  the  bandwidth  of  the 
modulation  signals  A(t)  and  B(t). 
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This  filtered  signal  is  fed  into  a  synchronous  detector  along  with  a 
reference  signal  given  by 


VREF(t>  "  "e  6 


The  detector  output  is  the  low-frequency  component  of  the  product  of 
these  two  voltages.  It  can  be  written  as 

i  *  j°V 

VSD(t>  fe  A(t)  B  (t~T)e  ’ 


Finally,  this  signal  is  fed  into  an  averaging  filter.  If,  for 
simplicity,  we  assume  that  the  averager  has  the  following  ideal,  (but 
physically  unrealizable),  impulse  response 

Vo  ^  =  2¥~  ’  f°r  0  <  t  <  2  T, 

=  0  }  otherwise, 

the  averager  output  and  indeed  the  system  output,  after  2T’ seconds  of 
averaging,  is 


SAB(T’T> 


2T 


A(T-t1)  B* 


(T-t  -T)dt 


(48) 


(49) 


(50) 


o 


(51) 
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Letting  t  =  T-t^  this  becomes 


R 


AB 


(T,t) 


A(t)  B*(t-T)dte 


(52) 


In  the  limit  as  the  system  output  becomes 


lim 
T  — » 


Rat, 
oo  AB 


(T,*j 


A(t)  B*(t-T)dte 


(53) 


Comparing  this  result  with  Equations  (39)  and  (41)  see  that 


li.  S’  (T,r,  =  Rr  (T)  .  i  Re  5  (T)  (54) 

T  * — >  ^ 

where  R  ^(t)  is  the  RF  complex  cross-correlation  function.  Consequently 
by  time-averaging  the  output  of  the  synchronous  detector  we  obtain  half  of 
the  real  part  of  the  RF  complex  cross-correlation  function  which  is  defined 


by  Equation  (42). 


To  obtain  the  imaginary  part  we  simply  perform  the  same  cross¬ 
correlation  of  the  two  signals  a(t)  and  b(t)  except  that  this  time  they 
are  taken  in  phase-quadrature  *  More  specifically  we  cross-correlate 
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t 

a(t)  =  Re  A (t) e 


(55) 


and 


J(W  t  +  #2) 
b'  (t)  =  Re  B(t)e  ° 


jw  t 

as  Re  j  B(t)e 


(56) 


Letting 


B1  (t)  =  j  B(t) 


(57) 


the  correlator  output,  as  T“ is 


R 


ABT 


<^A(t)  B’*  (t-T) 


(58) 


which  becomes 


rr 


A(t)  B*(t-T)^> 


(59) 
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Consequently  we  can  form  the  complex  sum 


i 


=  2 

_rab(t)  +  jrab(t)J 

(60) 

jco  T 

=  <JV(t)B  (t-T)^>e 

(61) 

to  obtain  the  complex  cross-correlation  function  of  the  two  complex  RF 
j«  t  jw  t 

signals  A(t)e  and  B(t)e  .  As  before,  we  can  multiply  both  :  Ides 

t 

of  Equation  (61) by  e  and  obtain 


Rab(t)  =  <^A(t)  B*(t-T)> 


(62) 


which  is  the  complex  cross-correlation  function  of  A(t)  and  B(t). 
4.3  The  Autocorrelation  Functions 


The  complex  autocorrelation  function  of  say  A(t)  is  defined  to  be 


raa(t) 


lim  1  j  ^  A(t)  A* (t-T)  dt 
T — *oo  2T 


=  <^A(t)  A*(t-T)]^> 


(63) 


It  can  be  measured  in  the  same  manner  as  the  cross-correlation  function, 
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To  do  this,  we  divide  a(t)  into  two  equal  parts  and  “cross- 
AB  7 

correlate*'  the  two  signals  as  if  they  were  distinct.  The  complex 
correlator  output,  except  for  the  factor  1/4,  will  be  R  CO^  as  given 
by  Equation  (63)  above.  The  complex  autocorrelation  function  of  the 
second  signal  B(t)  is 


RBB(T)  =  <B(t)  B*(t~T)> 


(64) 


and  it  too  can  be  measured  by  the  correlator  which  is  described  in  the 
preceeding  section. 

4.4  The  Correlation  Matrix 

Let  us  consider  the  two  signals  A(t)  and  B(t)  as  the  elements  of  the 
following  two-dimensional  complex  vector 


A  (t) 


& 


(t) 


B(t) 


(65) 


The  correlation  matrix  associated  with  the  two  signals  is  defined  as 


— *  -Xt 

R(t)  =<C  (t)  &  (t-T)l> 


(66) 


where  the  dagger,  indicates  the  transpose  conjugate  of  the  vector. 
Using  Equation  (65)  it  can  easily  be  shown  that 
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R(t) 


<^A(t)  A*(t-T)]j]><^A(t)  B*(t-T)]> 


<J3(t)  A*(t-T)' 


'B(t)  B*(t-T)^> 


(67) 


and  in  view  of  Equations  (62),  (63),  and  (64),  this  becomes 


R(t) 


(T) 

rab(t) 

(T) 

rbb(t) 

(68) 


The  second  cross-correlation  function,  R^Ct),  is  related  to  RAg(T)  by 
Hermitian  symmetry 


kab<t) 


“  ba(-t) 


(69) 


The  correlation  matrix  has  as  its  four  elements  the  four  correlation 
functions  that  can  be  obtained  from  the  two  signals  A(t)  and  B(t).  It 
completely  characterizes  the  correlation  properties  of  these  signals. 


37 


4„5  The  Effect  of  Finite  Averaging  Time 

In  practice,  of  course,  one  cannot  go  to  the  limit  of  the  infinite 

averaging  time  which  is  formally  required  by  the  correlation  integral 

(see  Equation  (37)  )0  It  is  of  some  interest  therefore  to  investigate 

the  effects  of  a  finite  averaging  time  2T_,  on  the  system  output 0  The 

latter  generally  will  differ  from  that  which  results  in  the  limiting 

case  of  an  infinite  time  average.  We  will  obtain  expressions  for  the 

expected  value  and  variance  of  the  complex  cross-correlation  system  output 

(31 ) 

as  functions  of  T„  In  an  earlier  analysis  ,  Linder  has  dealt  with  the. 
effect  of  finite  averaging  time  on  the  statistics  of  the  real  output  of 
an  RF  correlation  detector.  Since  what  we  present  here  is  intended  as 
only  an  outline  of  the  problem  we  refer  the  interested  reader  to 
Linder* s  paper  for  a  more  detailed  analysis. 

Assuming  that  the  signal  voltages  a(t)  and  b(t)  are  ergodic  random 
processes,  we  can  write  the  expected  value  of  the  complex  system  output 
after  2T  seconds  of  averaging  as 


E 


A  (t)  B*  (t-T)dt  l 


H 


E  \  Aft) 


* 

B  (t- 


T)  >  dt 


(70) 


(71) 


=  E 


Aft) 


B*(t-T) 


(72) 


J 
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where  E  is  the  expected  value  (mathematical  expectation)  of  f(t). 

But  by  the  ergodic  hypothesis  we  have,  with  probability  1, 


E  <A(t) 


B*(t-T) 


=  <^A(t)  B  (t-T) 


(73) 


Consequently,  in  view  of  Equation  (62),  we  can  write 


E 


rab<t) 


(74) 


which  is  the  desired  output  obtained  after  an  infinite  time  average. 

Thus  the  output  error  has  an  expected  value  of  zero,  which  is  independant 
of  the  length  of  the  averaging  time. 

We  will  now  determine  the  system  output  variance,  a  real  positive 
number  which  can  be  written  as 

°2ab(T’t)  =  e||rab<t’t,I2 

It  can  be  thought  of  as  the  fluctuating  arc .  power  of  the  output  which 

2 

exists  along  with  the  desired  d-c.  power,  |R^g(T)j  •  We  can  write  the 
first  term  on  the  right  side  of  Equation  (75)  as 


f'  IMrabCT'tV  I 


(75) 
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•H. 


E 


A (tl )  B*(tx-T)  A*(t2)B(t2-T)dt1cit2 


(76) 


j/r-t 


A  (t^  )B*  (t^—'OA*  (t2 


■T)) 

J 


dtldt2 


(77) 


Until  now  the  probability  distribution  of  the  random  signals  has  not 

been  specified.  It  could  be  quite  arbitrary.  However  if  the  signals  have 
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a  g^ussian  joint  probability  distribution  it  can  be  shown  that 


E  <A(t1)  B*(tJ-T)  A*(t2)  B(t2-T) 


E  <A(tx)  B*(t1-T)J-  E  ^  A*  (tg)  B(t2-T) 


+  E  < A(tl)  A  (t2)  }  E 


i 


B  (t1-T)  B(t  -T)f  (78) 
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and  by  using  the  ergodic  property  of  the  signals  we  have 


A(tl)  B*(t  *-T)A*(t  )  B(t2-r)/ 


bab<t)  "  ab<t> 


BAa‘VV 


R  „„ (t, -t  ) 
BB  1  2 


(79) 


Substituting  this  result  back  into  Equation  (77)  gives 


BAA(VV  "  BB<t1-t2>dtldt2 


+ 


bab<t) 


2 


:  (86) 


Jn  view  of  Equations  (74)  and  (75),  the  variance  can  be  written  as 


I 


1  * 

=  — —  If  R  (t  -t  )  H  (t  -t  )  dt  dt 
2  I  f  AA '  1  2  BBV  1  2  12 

4T  -T  -T 


(81) 


Letting  t^-t 


t  it  can  be  shown  that  we  can  reduce  the  above  double 


integral  to 


41 


2 


a 


AB 


(T,t) 


Re 


(l-  5t>  RAA(t)  B*BB<t)dt 


Note  that  this  is  independent  of  the  time  delay  t. 

The  above  expression  is  clearly  dependent  on  the  nature  of  auto- 
uuireia i,iun  functions  K^tT)  and  for  analytical  purposes, 

we  assume  that  the  signals'  power  spectra  are  both  uniform  in  the  band 
|w|  <  A  w  and  are  zero  outside,  then 


RA  A  (T)  =  R_(T)  =  R  sin  Awt 

AA  BB  o  - - - 

Awt 


The  output  variance  is 


(sin  A  w  t^  2 

r  a  « ■  t  i  dt 


which  after  some  manipulation  can  be  shown  to  be 


2 


=  R 

o 


Si  (4  A  co  T)  + 
2  A  w  T 


2 

R 

o 

2  (A  ^  T)2 


Ci(4A  T)  -1  -y  -ln<2T)i-  2A  ^  cos  4A  T 

4A  «  T 


(82) 


(83) 


(84) 


(85) 
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* 

where  Si  (x)  and  Ci  (x)  are  the  sine  and  cosine  integrals  respectively^ 
and  y  =.5772...  is  Euler's  constant.  This  is  a  rather  complicated 
result,  but  for  large  T  we  have 


T  5S> 


27T 

AW 


o 


Si(4fla  T) 
2  T 


T 


(86) 


The  output  variance  is  inversely  proportional  to  T  for  T  large. 

Although  this  has  been  shown  for  the  particular  case  of  uniform  power 

spectra  it  can  be  demonstrated  that  even  for  arbitrary  power  spectra 

2  23T' 

C  ^g(T,  t)  is  inversely  proportional  to  T  when  T  » 

2  If 

On  the  other  hand  for  T  small,  (T  «  —  ),  the  correlation  functions 

*  * 

in  Equation  (82)  are  essentially  constant  over  the  entire  range  of 
integration  and  we  can  write 


(J2(T,t) 


R  (o)  R  ,  (o) 

AA  BB 


T  « 


27 T 

AW 


(87) 


which  for  uniform  power  spectra  in  the  interval  to gives 


See  for  example,  Schelkunoff,  "Applied  Mathematics  for  Engineers  and 
Scientists,* "  Chapter  18,  Van  Nostrand,  1948. 
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0-2{T,t) 


T  « 


2 U 
Au 


~  R 

—  o 


(88) 


* 

Since  it  can  be  shown  that 


|RAB('T>|  «  (R 


AA(o)  RBB 


<=))  + 


(8  9) 


it  follows  that  the  standard  deviation  of  the  system  output  for  small  T 
is  equal  to  or  greater  than  its  expected  absolute  value.  For  meaningful 
results  one  must  therefore  take  reasonably  long  averaging  times.  Although 
we  have  considered  only  the  statistics  of  the  output  of  the  cross¬ 
correlator  it  is  abvious  that  the  results  apply  to  the  autocorrelation 
outputs  also.  Thus 


and 


raa(t) 


(90) 


(T,t) 


R  (t) 

AA 


2 

dt 


(91) 


Similar  expressions  obtain  for  the  autocorrelation  of  b(t). 


See  Davenport  and  Root  p.  61. 
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5.  THE  TWO-ANTENNA  CROSS -CORRELATION  SYSTEM 

We  will  not  consider  the  cross-correlation  of  the  voltages  obtained 
from  the  two  linear  antennas  when  they  are  excited  by  an  arbitrary  dis¬ 
tribution  of  remote  radio  sources.  It  will  be  recalled  that  the  source 
distribution  can  be  described  by  its  mutual  coherence  function  T(t,u,v). 

We  will  show  that  the  two-antenna  correlation  system  can  measure  this 
function.  The  patterns  of  the  two  antennas  are  scanned  independently, 
one  in  the  u  direction  and  the  other  in  the  v  direction,  and  when  the 
voltages  from  the  antennas  are  cross-correlated  we  obtain  a  complex  cross¬ 
correlation  output  function  which  can  be  written  as 

RAB(T,UjV)  “  T(TJu>v)  *  Cab(t,u,v)  <92) 

where  *  is  the  symbol  for  convolution  (in  three  dimensions),  and 
(t , u, v )  is  the  complex  system  function  which  will  be  defined  in 
Section  5,2, 

A  Fourier  analysis  of  this  output  will  be  made  and  from  the 

T  (T  u  v) 

analysis  we  will  be  able  to  define  a  principal  solution  AB  }  }  for 

the  cross-correlation  antenna  system.  This  principal  solution  is  a 

generalization  of  the  one-dimensional  principal  solution  T  (u)  which  was 

o 

18 

first  proposed  by  Bracewell  and  Roberts  in  connection  with  radio 
astronomy. 

5,1  Description  of  the  Antenna  Voltages 

A  diagram  of  the  two-antenna  cross-correlation  system  is  shown  in 
Figure  1.  Before  entering  the  cross-correlator,  each  of  the  terminal 
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voltages  a  (t,u)  and  b  (t.u)  of  the  two  linear  antennas  is  passed 
o  o  ’ 

j|C 

through  narrow-band  RF  filter  „  The  two  filtered  signals  Can  be 
represented  by 


A(t,  u) 


B(t,v) 


V(t"W  A(u1"u)du1  fA<tl)  dtl 


V(t~W  B(v1-v)dv1  fBCt2)dt2 


(93) 


(94) 


where  fA(t)  and  f  (t)  are  the  envelopes  of  the  impulse  responses  of  the  two 

A  o  , 

JWot 

filters  and  the  RF  carrier  factor  e  has  been  suppressed.  Note  that  the 

patterns  of  antennas  A  and  B  are  scanned  independently  in  the  u  and  v 
directions  respectively.  The  source  density  function  V(tjU)j  weighted  by 
these  pattern  functions^  is  integrated  over  the  visible  range  to  yield  the 
terminal  voltages  which  then  are  convoluted  with  the  filter  response 
functions.  If  we  define  the  mirror  image^  or  reverse^  of  a  function  to  be 


g(x)  =  g(-x)  (95) 


we  can  write  the  above  equations  as 


We  have  tacitly  assumed  that  the  frequency  response  of  each  antenna  itself 
is  essentially  constant  over  the  bandwidth  of  its  RF  filter. 
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A  (t, u) 


V(t,u) 


fA(t)  A (u) 


(96) 


V 

B(t,v)  =  V(t,v)  *  fB(t)  B(v)  .  (97) 

The  convolution  sign  applies  to  both  variables,  t  and  u  in  Equation  (96), 
and  t  and  v  in  Equation  (97) . 

These  two  equations  show  that  each  antenna  with  its  associated  RF  filter 
acts  as  a  combination  spatio-temporal  frequency  filter  of  the  two- 
dimensional  space-time  signal  V(t,  u).  As  is  indicated  in  Figure  1  it  is 
the  two  output  voltages  from  these  filters  that  will  be  cross-correlated. 
5.2  Cross-Correlation  of  the  Antenna  Voltages 

Comparing  Equations  (35),  (36)  and  (96),  (97)  we  see  that  the 

complex  cross-correlation  output  function  :t>£  the  antenna  system  can  be 
written  down  by  replacing  A(t)  and  B(t)  by  A(t,  u)  and  B(t,  v)  respectively, 
in  the  various  correlation  equations.  Thus  the  complex  system  output, 
after  2T  seconds  of  averaging,  is 

RAB(t,W)  =  2 k~fT  A(t;u)  dt  (98) 


and  the  expected  value  is 


RAB(T’U,V)  =  <CA(t>u)  B*(t-T,v)> 


(99) 


Substituting  the  expressions  for  A(t^u)  and  B(t^v)  given  by  Equations 
(96)  and  (97)  into  the  above^  we  obtain 
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RAB(T>U'v) 


-<i 


V(t,  u)  * 


tA  (t> 
A 


A(u) 


)(' 


V  (t-T,v) 


*  V* 

fB  (t-T)B(v) 


1> 


(100) 


UX)V  (t-T-tg, 


.  v/  (101) 

♦  n/  $ 

f. (t  )f  (t  )dt  dt  A(u-u  )B  (v-v  )du  dv 

l\  1  t>  i  1  e,  1  111 


In  view  of  Equation  (12)y  we  can  write 


<V(t-VV  V+(t-r-t2,V]L^>  =  T(T-t1+  t£ 


u 


’V 


(102) 


and  Equation  (101)  becomes 


rab(t'u>v) 


\J  V* 

T(tjU)v)  *  fABtT)  A(u)  B(v) 


(103) 
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where 


f  (t)  - 
AB 


VT) 


V 

f 


(r) 


vj  04  ) 


Thus  we  can  write 


rab(t<u’,)  =  1'l,»u.v> 


cab<t’u'v) 


(105) 


where  we  have  defined  the  complex  cross-correlation  system  function  to  be 


CAB(T’U’V)  =  fAB(T)  ^(U)  ^(V)  * 


1106) 


Equation  (.105)  is  the  fundamental  equation  of  the  cross-correlation 

process.  It  shows  how  the  system  function  C Ct}UjV)  operates  on  the 

source  function  T(r?u,v),  (by  convolution  in  three  dimensions),  to  give 

the  output  complex  cross-correlation  function  R  (tsu, v)r 

A  o 

5»3  Courier  Analysis  of  the  System  Output 

In  a  manner  similar  to  that  which  led  to  the  definition  of  the  com* 

bined  spatio-temporal  spectrum  of  T(t?u;v)  as  given  by  Equation  (24)  of 

Chapter  3,  we  define  the  spatio-temporal  spectrum  of  the  complex  cross  - 

correlation  system  output  R*^(T,u,v)  to  be 

AB  9  } 
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rAB(^X’y)  = 


(2JT) 


3/2 


OO  oo  <30 

/// 

-on  ~oo  -eo 


-  j (UT- ux+vy ) 


RAB(T’U<V)e 


dTdudv  .  (107) 


Substituting  the  expression  for  R  (t u  v)  of  Equation  (105)  into  the  above 

AB 

gives 


rAB{U’*’y)  " 


(2TT) 


3/2 


CO  BO  00 

1  1 1 

-oo  -bo  -oo 


(108) 


We  can  invbk^  the  convolution  theorem  to  write 


3/2 

rAB(Uj,X,y)  “{2ir)  t(w,x,y)  cAB(wJx,y)  (109) 

where  cA0(w>x>y)  is  the  inverse  Fourier  transform  of  the  system  function 
CAb^Tj u> •  Using  the  expression  for  C  (T,u>v)  given  by  Equation  (106) 
we  obtain 

oo  oo  oo 

Wu'*’y>  -  I3/2  /  /  /  fA<T)  *^b(T)  ^'u)  ®*v) 

\  dZ  U  J 

-00  -00  -oo 


-  j  (CJT-UX+  vy) 

o 


dT  dudv  . 


(110) 
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By  again  using  tb^  convolution  theorem  we  can  write 


cab^> x> y)  -  fb*^  a  b*^ 


on) 


where  F,  (u)  and  F  f(j)  .  the  inverse  Fourier  transforms  of  f  (T)  and 
A  B  *  A 

respectively,  are  the  frequency  response  characteristics  of  the 
two  RF  filters,  and  where  a (x)  and  b(y)  are  the  aperture  weighting 
functions  of  the  antennas.  Thus 


rAB 


(u,x,y)  =  47 r  t^x^y)  FA(w) 


Fb  (w)  a  (x) 


b*(y) 


=  (270  t  (w,  x,  y)  eAB  (w,  x,  y )  . 


(112) 


Equation  (112)  is  the  fundamental  equation  in  the  frequency  domain. 

It  shows  that  the  spectrum  of  the  output  is  proportional  to  the  spectrum 
of  the  mutual  coherence  function  weighted  by  the  system's  frequency 
response  function  which  is  fact6rable. 

5.4  The  Principal  Solution 

Our  main  purpose  is  of  course  to  measure  T(T,u,v),  or  equivalently^ 

its  spectrum  t(u,x,y),  However,  if  the  bandwidth  of  the  RF  filters 

are  both  say  2 AtJ,  it  follows  that  F  (ui)  F  (w)  —  0  for  |td|  >  The 

A  B 

system  passes  only  those  temporal  frequencies  of  the  source  distribution 
which  lie  in  the  narrow  RF  band  Igj— cjq  I  <  To  investigate  the 

pass-bands  of  spatial  frequencies  of  the  system  we  note  from  Figure  1 
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that  a(x)  =  0  for  Ix|  >  LA/2  and  b(x)  =  0  for  |x— M  >  L^/2  where  and 

are  the  lengths  of  the  two  apertures.  Consequently,  it  is  clear  that 

the  Cartesian  product  a(x)  b  (y)  =  0  except  when  |x|  <  L^/2  and  ly^i  <  itg/2. 

In  Figure  3  this  region  in  the  xy  plane  is  shown  cross-hatched  and 

labeled  Q  .  Note  that  the  spacing  &  between  apertures  causes  the 

region  to  move  away  from  the  origin.  Finally  we  can  define  a  three- 

dimensional  region  which  delimits  the  spatio-temporal  "aperture”  of 

the  system*  This  region  is  shown  in  the  diagram  of  Figure  4.  Strictly 

t 

speaking,  since  the  functions  F  (c*0  and  F  (<j)  are  analytic  in  a  half 

A  D 

4 

plane  ,  they  are  not  identically  zero  outside  of  Q  .  However  in 

Ad 

practice  there  is  negligible  error  in  assuming  that  they  are  zero. 

Now  from  Equation  (112)  we  see  that  since  x,y)  is  zero  for 

points  outside  Q^,  the  output  spectrum  contains  information  only  about 
that  part  of  the  spectrum  t((J,  x,y)  which  lies  within  Q  .  This  being 

Ad 

the  case,  we  can  define 


rAB<<J'X'y) 

Wu>x’y}  =  7a - ; - ~  for  (w'x'y)e  qab' 

41T  CAB(U>*’y) 


=  0,  otherwise,  (113) 

where  (^,  x,y)  €  Q  means  that  the  point  with  coordinates  (u,x,y)  belongs 

Ad 

to  the  point  set  Q^g* 

The  Fourier  transform  of  tAr>(o/,x,y)  is  the  principal  solution 
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OQ  oO  co 


tab(t'u'v)  = 


g«3/2  /  //  tAB 


(w^x^y)  xu+yv 3 dudxdy 


(114) 


-JOQ  -CO  —CO 


It  is  a  filtered  version  of  the  true  distribution  T(t,  u,v)  with  no 
components  of  frequency  which  lie  outside  the  region  QAg.  However, 
those  components  which  are  present  are  identical  to  those  of  the  true 
source  distribution. 

5.5  Direct*  Measurement  of  the  Principal  Solution 

Equation  (109)  indicates  that  the  output  spectrum  x,  y)  is 

generally  a  distorted  version  of  the  true  spectrum  t(u,  x,  y) .  Because 
the  system  characteristics  F  (w)  F  (to)  a(x),  and  b(y)  are  not  only 
zero  outside  of  Q,„  but  can  also  take  on  arbitrary  values  inside  the 
region,  it  was  necessary  to  divide  the  output  spectrum  by  the  system 
spectrum  (see  Equation  (113)  )  to  recover  the  principal  solution  spectrum 

x>y)- 

However,  if  the  system  characteristics  were  all  uniform  within 


Q 


AB' 


±.e.}  if 


cAB(“>x>y)  =  q«Hx,y)  = 


AB 


a*2 


for  (u,x,y)  C  Qab, 


(115) 


=  0  ,  otherwise, 

where  q^(w,  x,  y)  is  defined  as  the  system  function  which  is  uniform  in 
Qab>  then  from  Equation  (113)  we  have 
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rAB(w>  x'y)  =  for  e  QAB*> 


(116) 


=  0  *  otherwise. 


The  Fourier  transform  of  q(w,  x,  y)  is 


where  we  have  used  the  "sine"  notation  of  Woodward  ^ 


(117) 


sine  x  2= 


sin  2T  x 
H  x 


Comparing  Equations  (106)  and  (117)  we  see  that 


2  §  \ 
f  (r)  =  2  £&  sine  (£prr\ 

B  I  \  ir  J 


(118) 


and 


V‘>  -  1B( t) 


2  sine 
7T 


(119) 


The  responses  fA(t)  and  fg(t)  are  Physically  unrealizable  since  they  are 

non-zero  for  t  <  0.  However,  if  we  are  willing  to  accept  a  delay  of  t 

23T 

seconds,  with  t_^»  —  ,  between  the  input  and  the  desired  output,  then  to 
a  good  degree  of  approximation  the  response  functions 

*  ^ 

See,  for  example,  Davenport  and  Root  ,  p,  174. 
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yt-V  =  fB(t-Tl)  =a(|^AW  sine  Ct-r^) ]  (U9. a) 


can  be  synthesized.  The  two  types  of  response  f  (t)  and  f  (t-T  )  are  shown 

A  A  1 

in  Figure  5.  Since  f  (t-T )  —  0  for  t  <  0  there  is  usually  negligible 
error  in  having  the  actual  response  envelope  identically  zero  for  t  <  0. 

Consequently,  if  the  system  impulse  response  envelopes  and  patterns 
are  given  by  Equations  (117),  then  the  spectrum  becomes  (approximately) 


45T2Fa(w)  F*(co)  a(x)  b(y)  =  47T2  - 


jC0(Ti-Ti) 
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-fL,for  (u,x,y)  e  , 


=  0  }  otherwise  . 


(120) 


The  phase  factor  e  is  due  to  the  time  delay  necessary  to  insure  a 

physically  realizable  system.  Since  we  have  the  Hermitian  product  of 

jcor  +  -JCOT 

the  two  filter  spectra  F  (w)e  and  F  (co)e  ,  this  phase  factor 

cancels  out  in  the  final  result  (a  rather  fortunate  occurrence) ,  It  is 

interesting  to  note  that  the  analogous  condition  of  physical  realizability 

does  not  obtain  in  the  space  domain.  The  antenna  pattern's  reverse, 

v 

A (u),  is  actually  the  response  to  a  point  source  (spatial  impulse)  in 

the  u  =  0  direction.  Since  most  patterns  are  real  and  have  even  symmetry, 

v  i 

it  is  clear  that  A(u)  ^  0  for  u  <  0.  Physically,  we  might  "explain" 
this  by  saying  that  we  can  go  forward  and  backward  in  space  but  only  forward 
in  time.  However  in  practice  there  is  a  very  important  physical  condition 
that  must  be  satisfied  in  the  spatial  domain,  namely,  the  antenna  must  be 
of  finite  length.  In  a  sense  this  is  an  even  more  stringent 
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Figure  5.  Impulse  Responses  of  Realizable  and  Unrealizable 
Time-Domain  Filters . 
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condition  than  that  in  the  time  domain.  It  allows  the  aperture  function 
to  be  non-zero  only  over  a  finite  region  (e.g.,  [x|<C  L^/2) }  whereas  the 
impulse  response  of  a  time  filter  can  be  non-zero  over  the  semi¬ 
infinite  region  t  ^  0„ 

Returning  to  the  system  with  a  uniform  frequency  response  in  the 

22 J 

region  Q  we  see  that  to  a  good  degree  of  approximation  (if  r.»  — ) 

^  AB  1 


rab(t'u'v) 


*  tab<t'u'v) 


T(t,u,v)  *  Qab(t,u,v) 


(121) 


which  is  the  principal  solution. 

5,6  Cross-Correlation  System  Outputs  for  Coherent  and  Incoherent  Source 
Distributions 

In  Equation  (28)  in  Chapter  3  it  was  indicated  that  a  coherent 
quasi-monochromatic  source  distribution  has  the  following  mutual  coherence 
function^ 


T(t,u,v)  =  V (u)  V*(v),  for  t  « 


where  is  the  bandwidth  of  the  signals  and  satisfies  the  inequality 

«  co  „  Equation  (29)  indicates  that  the  spatial  frequency  spectrum  of 
o 

the  above  function  is 
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t  (r,x}  y)  =  v(x)  v*(y),  for  T  «  2ir/^U  • 


Both  the  mutual  coherence  function  and  its  spatial  frequency  spectrum 
are  factorable  for  small  T. 

Now  the  system  output  in  the  presence  of  such  a  distribution  is 


rab  (t<^v) 


=  V(u)  V  (v) 


«? 


y/  \/4c 

*  fAB(T)  A<U)  B(V) 


(122) 


«  K  V  (u)  *  A(u)  V*(v)  *  B*(v) 


(123) 


where  K  is  a  constant.  The  output  function  is  independent  of  T  and 
separable  (for  T  «  2 tt/Aj)  . 

Consequently  we  can.  write  the  spatial  frequency  spectrum  of  the 
output  as 


r 

AB 


=  K  21T  v(x)  a(x)  v*  (y)  b*(y) 

T  «  2tt/£u 


(124) 


i.e.,  this  spectrum  is  also  separable  and  independent  of  T  for  small  t. 

If  the  source  distribution  is  incoherent  then  we  can  write  (see  Chapter 
3}  Equation  (30) ), 


T(T,u,v)  =  T(T,u,u)  S(u-v) 

as  its  mutual  coherence  function.  The  corresponding  spectral  density 


function  is 
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oO  J30  OO 

t(“,x,y)  — ^72“  J  J  J  T<T,U,V>  6(u-v)e:5(WT'UX+Vy) 

^.oO  — oO  iJbo 


dTdudv 


or 


t<%x-y>  -  — hys  f  f  1 

(2)7) 3/2  JJ 


T-00  ^00 


(125) 


which  is  a  function  only  of  the  difference,  (x~y),  of  the  spatial  frequency 
coordinates  *  If  the  system  has  a  uniform  weighting  function  in  the  region 
(?ABJ  then  its  output  spectrum,  and  the  principal  solution  spectrum,  will  be 

rAB(0J>x~y)  =  for  6  QAb> 

=  0  ,  otherwise,  (126) 

which  also  is  a  function  of  the  difference,  x-y,  of  the  spatial  frequency 
coordinates „ 

Combining  both  the  coherent  and  incoherent  cases  we  can  make  the 
following  statements: 

a)  If  the  source  distribution  is  coherent  the  spatial  frequency  power 

271/ 

spectrum  is  separable  for  t  « 

b)  If  the  source  distribution  is  incoherent  the  spatial  frequency 
power  spectrum  is  a  function  only  of  the  difference,  (x-y),  of 
the  spatial  frequency  coordinates « 


In  practice  one  can  observe  the  spectrum  only  in  the  "aperture" 

Q  Although  t(w,  x,y)  is  an  analytic  function,  and  in  theory  can  be 

nD 

extended  to  points  outside  Q  by  analytic  continuation,  it  has  been 

nD 
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shown  that  in  the  presence  of  measurement  errors  very  little  meaning¬ 
ful  continuation  is  possible.  Consequently  even  if  t(w,x,  y)  is 
experimentally  observed  to  be  "coherent11  or  "incoherent  within  one 
cannot  be  sure  that  it  will  continue  to  be  so  outside  the  aperture. 
However  from  a  practical  point  of  view,  if  the  spectrum,  when  measured 
within  Q.  satisfies  either  of  the  above  conditions,  it  is  reasonable 
to  infer  that  the  sources  are  coherent  or  incoherent  as  the  case  may  be. 
5.7  Simplification  of  the  Cross-Correlation  System  for  the  Case  of 
Incoherent  Sources 

If  the  sources  in  different  directions  are  incoherent,  it  is  of  no 
to  point  the  patterns  of  the  two  antennas  in  different  directions.  Only 
when  the  two  beams  are  pointed  in  the  same  direction  will  there  be 
appreciable  output  from  the  system.  Thus  if  we  set  u  =  v  and  substitute 
Equation  (30)  into  Equation  (105)  we  can  write 


RAB(T’U'V) 


=T  (t,  u,  u)  *  CAfi(T,u,u) 

u=v 


,oO  *oO  — CO 


II 

"too  *oo 


\/ 

T(T-t,u1,v1)  6(u^  Afu-Uj 


u-v^ )dtdu^dv^ 
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.00  „co 


■// 

-00-00 


V  V* 

T(T-t,  u^u^)  fAg(t)  A(u-u^)  B  (u-u.^)  dtdu^  • 


(127) 


Using  Equation  (34),  and  letting  Rab(t,  u,  u)  =  a[21T  Rab^TjU^  we  obtain 


V  V* 

*ab(t,u)  =  T(T,u)  *  fAB(T)  A(u)  B  (u) 


(128) 


as  the  two-dimensional  cross-correlation  output  for  the  case  of  an  incoherent 
source  distribution,, 

We  define  the  spatio-temporal  frequency  spectrum  of  this  output  to  be 

♦  vjc 

r ab  z )  =  27T  t(o^z)  FA(w)  ^  (to)  a(z)  *  b  (z)  (129) 


where  z  x-y,  (see  Equation  (33) )„  Since  the  system*  s  pattern  is  now  in 

the  form  of  a  product  of  the  individual  field  strength  patterns  A(u)  and 

B  (u)^  the  spatial  frequency  spectrum  is  in  the  form  of  the  convolution 

V* 

of  the  aperture  weighting  functions  a(z)  and  b  (z) *  Consequently^  if  the 


L.  L 

A  13 

aperture  functions  are  uniform  for  |z!^_~  and  ,  as  Is 

shown  in  Figure  6a,  then  the  spatial  frequency  spectrum  a(z)  *  b  (z)  is 


shown  in  Figure  6b.  The  spectrum  is  a  trapezoidal  function  and  weights 
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Figure  t>.  Aperture  Distributions  and  Spatial  Frequency.  Spectrum  (for 
Incoherent  Sources)  of  Two  Uniformly  Weighted  Apertures. 
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more  heavily  the  spatial  frequencies  near  the  center  of  its  "aperture" 

than  those  near  the  edges.  This  results  in  a  distortion  of  the  source 

spectrum  as  "seen"  at  the  system's  output.  Note  also  that  the  width  of 

the  "aperture"  of  the  system  is  given  by  the  sum,  LA  +  L^,  of  the  aperture 

widths  of  the  individual  antennas  and  is  centered  at  z  =-i. 

In  order  to  obtain  a  uniform  spatial  frequency  spectrum  one  must  use 

22 

a  compound  interferometer  ,  A  diagram  of  the  weighting  functions  and 

associated  spectrum  of  the  system  is  shown  in  Figure  7,  Instead  of  two 

uniformly  weighted  apertures  the  compound  interferometer  consists  of  just 

one  uniformly  weighted  aperture  and  a  simple  interferometer.  The  two 

isotropic  elements  of  the  interferometer  are  iocated  at  the  end  points 

of  the  second  aperture.  Mathematically,  they  can  be  represented  by 

Dirac  deltas  and  their  convolution  with  the  other  aperture  function  results 

in  the  uniform  spectrum  shown  in  part  b  of  the  diagram.  The  width  of  the 

"aperture"  is  also  given  by  L  +  L  and  its  center  is  located  at  z  =-i. 

A  B 

In  Chapter  3  it  was  shown  that  in  the  incoherent  case  T(t,  u)  is  a 

real  function.  Consequently,  its  spectrum  is  complex  symmetric,  i,e., 

* 

t(a>,  z)  =  t  (-a),  ~z).  This  means  that  by  measuring  t(«,  z)  on  the  interval 
2|z  +  f  l<  La  +  Lb  ,we  can  automatically  deduce  its  values  on  the  interval 
2]z-f|'^  1<A  +  Lg  .  This  region  is  also  shown  in  Figures  6b  and  7b. 

Finally,  we  can  define  a  system  "aperture"  Q*  for  the  incoherent 

Ad 

case  which  is  analogous  to  the  more  general  QAfi  of  the  partially  coherent 
case.  Thus  in  Figure  8  is  shown  the  combined  spatio-temporal  frequency 
plane  for  the  incoherent  case  with  the  region  Q1 A^  indicated  by  cross¬ 


hatching. 


Figure  7.  Aperture  Distributions  and  Spacial  Frequency  Spectrum  of  the 

Compound  Interferometer  in  the  Presence  of  Incoherent  Sources. 


66 


Figure  8*  Spatio-Temporal  Frequency  Plane  for  the  Case  of  an  Incoherent 
Source  Distribution  with  the  Region  Q*  Shown  Hatched, 
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In  summary ^  we  can  say  that  an  incoherent  source  distribution  has  a 

principal  solution  spectrum  which  is  a  function  only  of  the  difference 

z  -•  x-y  of  the  spatial  frequency  coordinates.  Consequently^  the  problem 

reduces  from  three  dimensions  to  two  and  by  pointing  both  antenna  patterns 

in  the  same  direction  (u  -  v)  the  system  output  also  reduces  to  a  function 

of  two  dimensions^  delay  and  beam  direction  u.  A  principal  solution 

can  be  defined  whose  frequency  spectrum  is  complex  symmetric (t  (o)^  z)  - 

L 

A  R 

and  for  the  above  antenna  system  is  identically  zero  for  U-lill  >  — 

Ci 

and  |w|<  Indeed^  it  was  this  type  of  principal  solution  which  was 

18 

first  defined  by  Bracewell  and  Roberts  in  connection  with  the  mapping 
of  the  incoherent  sources  which  are  encountered  in  radio  astronomy. 
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6.  THE  CORRELATION  MATRIX  OF  THE  ANTENNA  SYSTEM 

It  will  be  recalled  that  in  Section  4.4  we  introduced  the  correlation 
matrix  R(T)  which  is  associated  with  the  complex  envelopes,  A(t)  and 
B(t),  of  the  narrow-band  voltages  a(t)  and  b(t)*  We  will  now  show  that 
for  the  two-antenna  correlation  system  we  can  obtain  a  correlation 
matrix  R(t,u,  v)  which  yields  considerably  more  information  about  the 
source  distribution  than  the  original  cross-correlation  function 
6.1  Cross-Correlation  of  Signals  from  Two  Coincident^  Antennas 

Let  us  consider  one  of  the  antennas,  for  example  antenna  A.  If  the 
signal  from  each  of  its  elements  is  divided  into  two  parts,  with 
predetermined  weighting  coefficients,  one  obtains  two  sets  of  signals 
which  can  be  combined  adaitively  to  form  two  distinct  output  voltages 
from  the  antenna  aperture.  In  effect  they  are  the  signals  from  two 
coincident  antennas,  and  A^,  whose  patterns  can  be  independently 
Scanned  in  the  u  and  v  directions  respectively.  These  signals,  after 
each  has  passed  through  an  RF  filter,  can  be  cross-correlated  in  the 
same  fashion  as  the  signals  from  the  distinct  apertures  A  and  B.  The 
complex  cross-correlation  function  that  results  is 

RAA(T'U'V)  =  T(T>u>v)  *  Cm(t„;v)  (130) 


where 


caa(T’^v) 


\/*  V  \/  * 

f  (T)  *  f  (T)  A  (u)  A  (v)  . 
A1  A2  1  2 


(131) 
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A  similar  cross-correlation  function  of  two  outputs  ffom  antenna  B 
can  be  obtained  and  it  is  given  by 

Rbb(T,u,v)  =  T(T,u,v>  *  Cbb(T,u,v)  (132) 


where 


CBB<T'U'V)  =  fB  <T)  *  V(T)  Bl(u)  B2*(v)  ’  <133) 

1  2 

6.2  The  Correlation  Matrix  of  the  Two-Antenna  System 

By  analogy  with  Equation  (67)  we  can  define  the  correlation  matrix 
of  the  antenna  system  to  be 


R(T,u,  v) 


EAA(T’U’V)  rab<t'u''° 
RBfl<T>U’V)  Rbb<T'U,V)  - 


(134) 


As  before  the  cross-correlation  function  R^a(t,u,v)  is  related  to  Rat_{t.u.v) 

BA  AB  * 

by  the  formula 


rba(t'u'v)  =  rab  <-t'v'v)  • 


(135) 


Now  just  as  in  the  case  of  R  (T,u,v), 

AJd 


we  can  take  the  Fourier 


transform  of  this  correlation  matrix  and  we  obtain 
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Each  of  the  elements  of  the  frequency  domain  correlation  matrix  can  be 

* 

non-zero  only  where  the  corresponding  system  function,  e.g„,  F  (u)F  (u)b  (x)b  (yj 

B1  B2  1  2 

is  non-zero.  Let  us  assume  that  the  temporal  frequency  pass-bands  are  all 
equal  in  size.  We  can  then  turn  our  attention  to  the  spatial  frequency 
domain*  Figure  9  shows  the  regions  in  the  xy*  plane  where  the  four 
distinct  system  functions  are  non-zero.  Note  that  is  the  "aperture" 

of  the  original  cross-correlation  system  (see  Figure  3) .  The  three  other 
regions,  Q^,  and  QgA,  indicate  the  additional  "apertures"  of  the 

system* 

6.3  Principal  Solution 

By  analogy  with  the  principal  solution  for  the  original  cross- 
cori elation  system  giver,  by  Equation  (114),  we  can  define  a  principal  ’solution 


i 

\ 
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♦  * 


I 


Figure  9.  Spatial  Frequency_Plane  on  which  the  Cross-Correlation  System* 
"Apertures”  ^AA>  QAB,  QgA,  and  QBB  are  Shown  Hatched. 


72 


for  the  correlation  system  as  a  whole  to  be  the  Fourier  transform  of 


t0<w,x,y)  =  tAA(wJxj  y)  tAB(tJ'x,y) 


+  tBA(W,Xjy)  +  tBB(UjX,y) 


(141) 


where 


rAA(^^y) 

4712  F  (u)F*  (w)a(x)a*(y) 
A1  A2 


if  (cj^x^y)  £  ^ AA # 


=  0 


) 


otherwise , 


(142) 


Similar  formulas  obtain  for  *  >y)  ,  y),  and  t^Cw^x^y)* 

Thus 


T  (T,u,v)  = 
o 


jCtjT-xu+yv) 


dudxdy 


(143) 


is  the  principal  solution  for  the  correlation  system  as  a  whole*  Clearly 
this  solution  gives  us  much  more  information  about  the  source  distribution 
than  does  the  cross-correlation  system* s  principal  solution  as  given  by 
Equation  (114) . 

The  special  cases  of  coherent  and  incoherent  distributions  give  rise 
to  the  same  type  of  outputs  for  this  more  general  system  as  for  the  one 


described  previously. 


f 

I 

Finally  <,  It  should  be  noted  that  if  the  two  apertures  A  and  B  are 
j  considered  as  parts  of  a  single  aperture  and  if  two  distinct  signals 

|  from  this  larger  aperture  are  cross-correlated,  as  described  in  Section 

6„1,  then  the  output  from  this  system  will  have  a  spatial  frequency 
^aperture"1  which  is  the  same  as  that  of  the  correlation  matrix. 
Consequently,  the  principal  solution  for  the  aperture  C  is  the  same  as 
the  system  solution  given  by  Equation  (143),  Therefore,  when  mapping  a 
source  distribution,  it  is  sufficient  to  obtain  a  cross-correlation 
function  from  the  total  available  aperture.  There  is  nothing  available 
in  the  cross-correlation  of  signals  from  parts  of  an  aperture  that  is 
not  present  in  a  cross-correlation  of  the  signals  from  the  entire  aperture. 
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7 .  COMPARISON  OF  THE  CROSS -CORRELATION  AND  THE  CONVENT  I ONAL  ANTENNA  SYSTEM 

x 

In  this  chapter  it  will  be  shown  that  when  the  terminal  voltage  of 
a  conventional  linear  antenna  is  square-law  detected  and  time-averaged, 
the  resulting  output  is  just  a  special  case  of  a  cross-correlation  antenna 
system.  This  means  that  a  comparison  of  a  linear  and  a  cross-correlation 
system  should  be  made  between  the  square- law  detected  output  of  the  former 
and  the  cross-correlated  output  of  the  latter  system, 

7,1  A  Linear  Antenna  with  Square-Law  Detection 

A  diagram  of  a  conventional  linear  antenna  with  a  square-law  detection 
system  is  shown  in  Figure  10 ,  After  being  passed  through  the  RF  filter 
the  antenna  voltage  can  be  written  as 

jcj  t 

d(t,u)  =r  D(t, u)  e  (144) 


where  by  analogy  with  Equations  (96 )  and  {97)of  Chapter  5,  we  have 


D(t,u)  =  V(t,u)  *  fD(t)  K(u)  - 


(145) 


The  antenna’s  pattern  is  D(u)  and  the  filter  response  envelope  is  f  ^(t) , 
This  filtered  signal  is  fed  into  a  square-law  device  whose  output  is 


Vt’u)=  t Re 


r 


(t,u) 


i D(t  u) r 


(146) 


This  signal  is  then  time-averaged.  Since  the  double  frequency  component 
of  the  square- law  detector  output  has  a  zero  average  we  obtain 
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i 

(dSL(,'Ul) 


Figure  10.  A  Conventional  linear  Antenna  which  Uses  Square-haw 
Detection, 


<dsL<,^>  =  !<lD(t'u) 
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(147) 

Now  it  will  be  recalled  that  the  cross-correlation*s  complex  output 
is 


rab(T'u'v)  =  <3(t>u)  B*<t-T>v^>  *  0-48) 

Comparing  these  two  results  we  see  that  except  for  the  constant  factor. 

1/2 9 the  output  from  the  square-law  detection  system  is  simply  a  degenerate 
case  of  the  cross-correlator  output  which  occurs  when  A(t^u)  -  B(t,u)  =  D(t^u). 
It  is  half  Ihe  Autocorr e lat ion  function  of  the  signal  from  antenna  D'  evaluated 
at  t  -  o. 


2  <*Sl 


v  v* 

T(T  u  v)  *  f  CO  D(u)  D  (v) 
*  *  ©D 


T=0 

u=v 


(149) 


where 

fDD(T)  =  fD(T)  *  V(T)  <150) 

7.2  Disadvantages  of  the  Conventional  System 

Since  the  input  to  the  square-law  detector  is  the  terminal  voltage 
of  only  one  antenna  with  only  one  scan  parameter  u,  we  see’  that' when 
the  source  distribution  is  partially  coherent^  It  is  impossible  to-obtain  the 
three-dimensional  mapping  of  the  mutual  coherent;  function  of 
the  source  distribution*  As  indicated  by  Equation  (149)  one  can  only 
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obtain  the  values  of  the  antenna  voltages  autocorrelation  function 
along  the  single  line,  t  =  0,  u  =  v,  whereas  the  cross-correlation 
system  employing  a  single  aperture  yields  the  values  of  the  filtered 
function  for  all  points  in  the  three-dimensional  output  space  of  the 
system* 

In  the  case  of  incoherent  sources  the  problem  is  not  nearly  as 
serious  since  even  with  the  correlation  system  it  is  only  the  u  =  v  plane 
that  is  considered  when  the  source  distribution  is  scanned*  (See  Chapter  5, 
Section  7)*  Indeed  the  only  difference  in  the  space  domain  between 

the  two  outputs  is  that  the  power  pattern  of  the  conventional  system 

2  * 
is  |D(u)|  and  that  of  the  cross-correlation  system  is  A(u)  B  (u) 

(compare  Equations  (103)  and  (149)) 0  The  latter  pattern  is  clearly  more 

flexible  since  it  is  the  product  of  two  distinct  field  strength  patterns. 

It  degenerates  to  the  conventional  power  pattern  when  the  two  field 

strength  patterns  are  identical^  (A(u)  -  B(u)  =  D(u)), 

7„3  Comparison  of  the  Two  Systems  for  The  Case  of  Dolph-Chebyshev 

Array  Synthesis 

Let  us  first  consider  a  conventional  system  whose  antenna  is  in 

40 

the  form  of  a  linear  array.  It  was  shown  by  Dolph  that  the  optimum 
array  pattern  is  given  by 

D  (41)  =  T  (4>)  (151) 

n 

where  4*  -  ^  cos  sin  0), 

i  -  element  spacing, 

0  -  angle  measured  from  the  normal  to  the  array, 
n4-l  ~  number  of  elements  in  the  array. 
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T^{40  is  the  Chebyshev  polynomial  of  degree  n* 

The  value  of  4*  is  related  to  the  sidel&be  level  p  of  the  pattern  by 
o 

the  formula 


P  =  20  log  T  OM  .  (152) 

10  n  o 

The  pattern  is  optimum  in  the  sense  that  for  a  given  array,  of  all  patterns 
with  an  arbitrary  sidelobe  level  p,  the  Chebyshev  pattern  has  the 
narrowest  beamwidth,  or  conversely,  of  all  patterns  with  a  given 
beamwidth  it  has  the  lowest  sidejLobes.  0 

However,  if  the  terminal  voltage  of  such  an  optimum  array  is  square- 
law  detected,  then  the  detector  output  for  a  point  source  in  the 
direction  $  is  proportional  to 


IdoM*  )  I2  =  t  )  . 

s  ns 


(153) 


This  output  is  a  polynomial  of  degree  2n  but  is  not  the  Chebyshev 

polynomial  of  that  degree,  .  Thus  when  the  power  pattern  is 

considered  the  conventional  Dolph-Chebyshev  design  is  not  optimum. 

Let  us  now  turn  our  attention  to  the  cross-correlation  system  shown 

in  Figure  11.  The  array  is  given  the  weighting  coefficients  necessary 

to  produce  the  conventional  field  strength  pattern  T  OlO  *  However, 

n  7 

the  center  element  of  the  array  has  its  signal  split  into  two  parts;  the 
first  is  given  the  amplitude  -  (1  -  7?}  and  the  second  the  amplitude 
^  (1  +  «  The  signals  from  the  other  n  elements  are  combined  in  the 


*For  a  discussion  of^he  properties  of  these  orthogonal  functions,  see 
Courant  and  Hilbert  ,  pp0  88-90. 


CONVENTIONAL 
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Figure  11  *  System  Used  to  Obtain  the  Optimum  Factor  Patterns  for  the 
Dolpl.-Chebyshev  Product  Array. 
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usual  way  to  form  T^C^);  however,  since  the  central  element's  signal 
(normalized  to  1)  is  missing  we  actually  obtain  T^O^)  -  1.  This  signal 
is  split  in  two  halves  and  to  one  half  id  added  the  first  part  of  the 
center  element's  signal  to  form 

A<+)  =  1  (Tnd0  -  •  (154) 

Likewise  BOP)  is  formed; by  adding  the  remaining  two  signals. 

BQ\> )  =  -  (T(40  +  .  (155) 

Now  if  we  note  of  the  following  identity 

Trt  (^)  =  2T  20lJ)  -  1  (156) 

2n  n  > 

.and  if  the  above  pair  of  combined  signals  are  cross-correlated  with 
zero  time  delay  .It.  is  easy  -to i;  see :  that-  except  ■  for  a  constant,  factor  we^ 
obtain  the.  folldwilng- product  pattern 

R(40  =  A (+)  B*d0  »  T  (+)  .(157) 

2n 

For  the  point  source  in  the  direction  4*  we  have 

Rdo  =  T  (^  -  Y  )  (158) 

2n  s 

36 

as  the  system  output.  This  method,  due  to  Price  ,  yields  the  desired 
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result.  The  system  output  is  a  Chebyshev  polynomial  of  degree  2n  and  as 
such  it  has  the  lowest  sidelobes  of  any  possible  pattern  of  the  array 
when  a  certain  beamwidth  is  specified.  In  particular,  the  improvement 
In  the  sidelobe  level  of  this  cross-correlation  pattern  over  that  of 
the  conventional  Dolph-Chebyshev  pattern  (i.e.,  Equation  (153))  is 
shown  in  Figure  12 .  Both  patterns  have  the  same  beamwidths  measured 
to  the  first  null  and  the  element  spacing  is  V2.  The  improvement, 
plotted  as  a  function  of  the  sidelobe  level  of  the  convent  ibna'l  pattern, 
is  about  6  db  but  is  slightly  higher  for  small  arrays  and  relatively 
high  sidelobe  levels. 

It  is  worth  remarking  that  the  method  of  synthesizing  the  product 
pattern  T^Ol*),  which  we  have  just  described,  is  not  uniquel  To  show 
this  we  note  that  we  can  write 

2n 

1 

T2n(+)  =  22n_1  Tf  (^  -  +k)  (159) 

k=l 

where  ^  ,  the  kth  zero  of  the  polynomial , is  given  by 

1  2k  -  1 

=  C°S  -  *1*  k  =  h  2>  - ,  2n  .  (160) 

The  2n  zeros  determine  the  polynomial.  They  can  be  arranged  in  two  groups 
of  n  zeros  in 


W 


1  (2n) ! 

2  aryys 


(161) 


ways. 
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Since  one  group  of  n  zeros  will  determine  say  AOpK  and  the  other  will 

specify  B^),  we  see  that  there  are  W  possible  ways  to  synthesize  the 

product  pattern  T  (41)  =  A(40  B  (41)  .  It  is  natural  to  ask  if  there 

is  an  optimum  factorization  of  T  (40  ,  i.e. ,  which  of  the  W  possible 

2n 

pairs  of  factor  patterns  AC^)  and  3(40  is  best  in  some  way. 

If  the  system  is  being  used  to  receive  a  signal  from  some  known 
direction,  and  if  there  is  a  distribution  of  incoherent  noise  sources  in 
addition  to  the  signal,  then  the  system  output  will  be 


Rc,m(t,4o  =  [S  (T)8eM»  )  +  N(T,4»)]  *  f._CO  Ado  B  (+)  (162) 

brW  O  S  Ab 


where  S  (T)  )  and  N(T  4*)  are  the  signal  and  noise  coherence  functions 

o  s  9 

respectively.  It  is  clear  that  it  does  not  matter  how  the  zeros  of 

T  (40  are  shared  between  A(4o  and  B(4o  since  these  patterns  enter  the 
2n 

output  expression  simply  as  the  product  A(40  B*(40  =  T  (+)  „  Consequently 

2  ^ 

the  W  possible  ways  of  producing  the  optimum  product  pattern  are  all 

equivalent  when  the  sources  ar«  incoherent. 

This  is  not  the  case  for  coherent  and  partially  coherent  sources, 

In  particular,  let  us  consider  the  case  of  a  desired  signal  incident  from 

the  direction  4*  and  a  coherent  interfering  signal  from  an  arbitrary 

direction  4"  ,  In  practice  this  situation  occurs  when  there  are  two 
c 

distinct  paths  between  a  transmitting  and  a  receiving  antenna  (multipath 

x  51 

propagation) „  It  has  been  shown  elsewhere  that  if  the  relative 

amplitudes  of  the  desired  and  interfering  signals  are  1  and  Se^  respectively 

(§  is  the  relative  phase),  then  the  real  system  output,  when  T"=  O' and 

both  patterns  are  beamed  in  the  4?  direction,  is  proportional  to 

s 
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R  (0.+  )  =  1  +  2S 

J  s 


A  (t  -  )  +  B  0^  -  +  ) 

o  s  c  o  S  c  t 

- - — - - - cos  5 


+  S2  A  (+  -  +  )  B  <+  -  +  ) 

os  cos  c 


(103) 


where 


A  ob 
o 


a  eb 

A(0) 


B  cm  _  sob 

Bo(^)  “  B(0) 


(164) 


are  the  normalized,  real,  factor  patterns  of  T  (40  *  Notice  that  the 
second  term  in  the  above  expression,  which  is  due  to  the  cross-product 
of  the  desired  and  interfering  signals,  is  proportional  to  the  sum  of  the 
two  factor  patterns.  If  the  interfering  signal  is  low  level  (S  «  1), 
the  interference  rejection  will  depend  on  this  sum  pattern.  Now  if 
we  select  the  factor  patterns  in  the  manner  of  Figure  11  then  we  have 


1 

2 


[A  Ci'  -Y  )  +  B  (+-+)]  =  i 
o  s  c  o  s  c  2 


t  d*  -+  )  -  J 

c  JZ 


n 

t»<°>  - 


T  Q\>  )  + 


7S 


7? 


T  (0)  + 
n 


is 


(165) 


-  lf¥ 


T  <0)  _  1 

n 

2T  (0) 

n 


(166) 


If,  as  is  usually  the  case,  Tn(°)  »  we  can  write  the  sum  pattern  to 
a  good  degree  of  approximation  as 


~[  A  (4j-4')+B(4'*4i)]  = 


(167) 
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Consequently,  the  output  is 

_  /  T6Jj**Uj)  >■  2  T  fib  —  ill  ) 

R  (0,+  )  =  1  +  2S  nws  Yc  cos  5  +  S  '  2n,*B  Yc  (168) 

T„(0)  T2„«> 

and  we  see  that  not  only  does  the  system  reject  the  sdlf -product  of  the 
interfering  signal  in  the  optimum  Chebyshev  fashion 

but  the  cross-product  of  the  two  signals  is  also  rejected  in  the  optimum 

fashion  (T  (+  -4»  )).  This  occurs  only  if  the  factor  patterns  AOb  and 
n  s  c 

Bob  are  given  by  Equations  (154)  and  (155)  respectively.  In  this  sense, 

then,  this  factorization  is  the  optimum  one. 

Finally  we  should  note  that  since  the  pattern  of  any  uniformly 

52 

spaced  linear  array  can  be  represented  as  a  polynomial^  we  can  synthesize 
a  product  pattern  of  any  form  by  simply  determining  the  zeros  of  the 
pattern  and  distributing  them  equally  between  its  two  factor  patterns'?^* 

The  manner  in  which  the  distribution  is  made  could  be  determined  from 
considerations  of  interference  rejection  as  in  the  Chebyshev  case. 

The  patterns  of  two  coincident  arrays  of  n+1  elements  will  each  be  a 
polynomial  of  degree  n.  Their  product  pattern  will  be  a  polynomial  of 
degree  2n  with  2n  distinct  zeros.  On  the  other  hand  the  power  pattern 
of  a  conventional  array  of  n-f-1  elements  will  also  be  a  polynomial  of 
degree  2n  but  will  have  only  n  distinct  zeros.  Consequently  we  can 
say  that  the  product  pattern  has  twice  as  many  degrees  of  freedom  as  the 


conventional  power  pattern. 
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8,  MULTIPLE  CROSS-CORRELATION 


Having  studied  in  some  detail  the  cross>coi*relat;Lon  af  the  signals  of 
two  antennas  we  now  will  consider  the  effect  of  multiple  cross¬ 
correlation  of  the  signals  from  an  arbitrary  number  of  '  antennas.  By 

now  it  is  well-known  that  when  used  for  mapping  the  temperature  distribution 

20  3 1 

of  remote  sources  these  systems  are  basically  nonlinear.  *  ,  The  output 

of  a  four-antenna  system,  for  example,  will  contain  cross-product  terms  due 
to  sources  in  different  directions  even  when  their  signals  are  mutually 
incoherent.  The  explanation  for  the  failure  of  such  systems  is  simple. 

The  temperature  distribution  is  a  distribution  of  average  power  and  as  such, 
is  determined  by  a  second  moment  of  the  probability  distribution  of  the 
field  phasor  V(t,u),  i.e., 


T(T, u, v) 


(169) 


* 

where  p  [Vt  u,Vt_T  y]  is  the  joint  probability  distribution  of  the  field 

from  the  source  in  the  u  direction  at  time  t  and  the  field  from  the  source 

in  the  v  direction  at  time  t-T,  To  emphasize  the  fact  that  t,  T,  u,  and  v  are 

essentially  parameters  in  the  above  integral.,  w£  have ’for.  example, 'let  V  be  thfc 

value  of  V(t,u)  at  time  t  and  for  the  direction  u.  The  integration,  of 

course,  is  over  all  values  of  the  complex  amplitudes  V  and  V* 

t,u  t-T,v 

Now  by  going  to  a  higher  order  of  cross-correlation  we  are  no  longer 
dealing  with  second  moments  and  hence  by  definition  with  temperature 
brightness.  All  temperature  information  about  the  distribution  is  contained 
in  the  mutual  coherence  function  T(*",u, v).  The  information  in  the  higher  order 
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correlations  clearly  concerns  the  higher  order  moments  and  hence  can  be 
used  to  specify  the  joint  probability  distribution  of  the  sources^  e.g„, 

P{Vt-W  V  t-r2>u2}  Vt_T3jU3  Vt-Jn,un|  •  Conversely,  if  the 

probability  distribution  is  known,  all  of  the  moments  can  be  determined 
from  it  including  the  second  or  temperature  moment.  However,  it  will  be 
shown  in  Section  8*3  that  the  higher  order  moments  of  a  gaussian  distribution 
are  completely  determined  by  the  first  and  second  moments. 

841  A  Multiple  Cross-Correlation  Antenna  System 

A  diagram  of  the  system, whereby  the  terminal  voltages  of  N  antennas 
are  cross-correlated*  is  shown  in  Figure  13c  In  Appendix  B  a  detailed 
analysis  shews  that,  the  complex  system  output'  has  an  expected  value  of 
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where 


,(N) 
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H  *r 


with  rk  and  u^  being  the  time  delay  and  beam  direction  of  the  kth  antenna 
whose  pattern  is  A  (u)  and  i  liter  response  envelope  is  f  (t) 
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It  is  clear  from  Equation  (170)  that  a  Fourier  analysis  of  the  output 
will  be  simply  a  generalization  to  2N  dimensions  of  the  analysis  for  the 
case  of  a  single  cross-correlation.  Thus  the  system  output  spectrum  is 


_ An  <1*2f  “n'  XV  V - >  XN) 


vn  m 


*  1  V~'V  V*?  W  a3(X3J 


•4(v 


v  V* 

F  (u  )  F  (u  ) 
iv  1  2  2 


\/  * 
m  N 


(172) 


where 


CO  X 

'>  W  V 


V‘ 


a.  (x  ) ,  and  F,  (w  ),are  the  inverse  Fourier  transforms  of 
k  k  *  k  k  J 


(N) 

RA  .  .  .A  ^T1>T2^  "  *  '>rW  Ul,***>  UN^  *  T  ^^2' ‘  ‘ '>TN,lV  U2>  ' 
1  n 


W*  and  fk(T  t 


respectively. 

Letting  Q  be  the  support  of  the  system  function  a  (x  )  a*{x  ) . . 0 . 

o  o  o  ^  X  1  4  ^ 

5fC  V  +  \J  & 

aN^XN^  P2^W2^  *  *  ’  ’  *  * "  ‘  FN^N^  We  de^±ne  +he  Principal  solution  for 


the  Nth  order  correlation  system  as  the  Fourier  transform  of 
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=  0  ,  otherwise  , 


Thus  the  principal  solution  is 
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which  is  a  filtered  version  of  the  following  Nth  moment  of  the  joint 
probability  distribution 
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Note  that  for  N  even,  t  -  T  -  t  -  .  .  *=  T  =0  arid  u  -  u  =  u  =  . 

J  1  2  3  N  12  3 

=  uN  =  u  we  get 

(176) 


which  is  the  Nth  moment  of  the  probability  distribution  of  the  source  in  the 
u  direction 

It  should  be  remarked  that  for  the  above  results  we  have  tacitly  assumed 
that  N  was  cven„  If  N  is  odd  the  output  in  most  practical  cases  will  be 
zero  since  most  sources  have  even  probability  distributions  (e.g.,  gaussian) . 
However^  more  general  distributions,  are  possible  and  their  odd  moments 
can  be  measured  by  using  an  odd  number  of  antennas.  In  Appendix  B  we  give 
the  modification  of  the  correlation  system  that  will  perform  a  correlation 
of  an  odd  number  of  signals  . 

8.2  Generalized  Correlation  Matrix 

Just  as  it  was  shown  in  Chapter  6  that  a  cross-correlation  of  the  voltages 
of  two  antennas  yields  only  a  part  of  the  available  information  from  the 
system,  so  in  the  generalized  case  it  can  be  shown  that  the  multiple 

N 

correlation  output.  R  (T  , T  ,  .  .  .  ,T  ,  u  ,  u  .  .  * . ,  u  )  is  only  one  of  N 

1  ‘  '  *  N  1  2  "12  N 

possible  outputs,  all  of  which  are  distinct  functions  of  the  delay  and 

scan  variables.  For  example,  we  could  form  N  distinct  outputs  from  a 

single  antenna  and  obtain  an  N-dimensional  cross-correlation  just  as  we 

obtained  a  two-dimensional  cross-correlation  of  the  signals  from  antenna  A, 

(see  Chapter  6) .  By  analogy  with  that  output  we  could  let  R,  *  .  (t  ,  # 

A  A,  .  .  1-  *  r 

11  1  y 

*  - —  its  ^“dimensional  counterpart.  The  2x2 
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correlation  matrix  will  also  have  counterpart  and  the  typical  element  of  this 
generalized  correlation  "matrix”  is 


R 


A  A^A  ;  J  “A  ^  1* T 2’  *  s  * 9  Ui*  U2^  *  *  *  * 9  uj^  = 
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V*  V  V* 


•fq(V  \<Ul)  Ai(V  W  •‘•“VV 


(177) 


where  each  of  the  N  indices,  k,  m^  .  . . ,  q,  assume  values  from  1  to  N. 

Each  element,  RA  A  A  (T-, T?}  •  •  • ,  TM,  u., ,  .  .  „ . ,  uN),  of  the  generalized 

‘  k  i .  «  (N) 

correlation  matrix  has  a  principal  solution,  T  (T  u), 

A„  »  •  •  •  A  1  ft 

k  q 

associated  with  it.  We  define  the  principal  solution  of  the  Nth  order 
correlation  system  as  a  whole  to  be  the  sum  of  these  NN  principal  solutions, 
i.e.  , 
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where 
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is  the  Fourier  transform  of 
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(N) 
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and  rA  a  (W-  ,  *  . , .  ,  y  x  )  is  the  inverse  Fourier  transform  of  the  output 

A  4  O  «  9  ft  -i-  IN 

k  q 

cross-correlation  function  R,  ,  (T  , . . „ D 0 , u  ) . 

A,  .  . . .  .  A  1*  *  N 

k  q 

8,3  Sources  with  Gaussian  Statistics 

If  the  sources  have  a  gaussian  joint  probability  distribution  with  zero 
54 

means  it  can  be  shown  that 


v(t_TijUi)  v  v(t"Vu3),**,,y  Ct“' 
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is  even, 


~  Qy  if  N  is  oddj,  (180) 


where  a„d„p„  means  all  distinct  pairs  and 


V'  (t-Tifu^)  =  V(t-T^,u  ), for  i  odd^ 

*  „ 

=  V  (t-T  .  ,u.  ),  f  or  i  even 
i  i  ’ 


(181) 


94 


Thus  the  higher  order  moments  of  a  gaussian  distribution  are  all 
expressible  in  terms  of  the  sum  of  products  of  second  order  moments. 
For  example  if  N  =  4  we  can  write 


VU-W  V  (t-T^u2)  V(t~T3ju3)  V  (t-T4,u4) 


T(W  V  V  T  <VV  V  V 


*  T(VT1- 


u  ,  u  )  T  (T  -T  , 
l’  4  V  2  3 


u  ,  u  ) 
3J  2 


+  T’ ("T1  +  V 


O  T*  (-T+T 


U4) 


(182) 


where  TfT^u^v)  is  given  by  Equation  (19)  and 

Tf  (T.u.v)  =  E  jv(t,  u)  V  (t-T, 

If  the  ,  complex  signals  ^  have  zero-mean  statistically  independent  real 
and  imaginary  components  and  the  covariances  of  the  components  are  the  same, 

47 

then  T*(TjU,  v)  =  0,  Such  signals  represent  circular  complex  gaussian  processes  . 
This  result  is  simply  a  generalization  of  the  well-known  fact  that  if 

the  mean  arid  variance  (first  and  second  moments)  of  a  gaussian  random 

♦ 

variable  are  known,  then  all  higher  moments  can  be  deduced  from  them  . 

We  can  infer  that  if  such  is  the  case  no  further  information  about  the 
source  distribution  is  available  in  the  higher  order  correlations. 

Consequently,  if  we  knew  a  priori  that  the  sources  had  gaussian  statistics 
a  single  cross-correlation  of  the  signals  would  be  sufficient. 

♦See,  for  example,  Davenport  and  Root  ,  pp.  146-147. 
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Conversely,  if  the  measured  higher  order  correlations  are  related  to 
the  measured  second  order  correlations  in  the  manner  indicated  by 
Equation  (180),  then  one  can  infer  that  the  source  distribution  has 
gaussian  statistics. 
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9.  THE  APPLICATION  OF  CROSS -CORRELATION  ANTENNA  SYSTEMS  TO  RADAR 

In  this  chapter  it  will  be  shown  that  antenna  cross-correlation 
techniques  can  be  used  to  advantage  in  radar  systems.  In  the  radar 
case;  since  one  is  both  trahsmitting  '  hndu  receiving  signals,  it  is 
possible  to  use  two  antennas  for  transmitting  as  well  as  two  for  re¬ 
ceiving,  four  antennas  in  all.  To  distinguish  the  two  transmitted 

signals,  one  is  shifted  in  frequency  by  the  amount  w  .  with  «<o>.  The 

r  i  o* 

receiving  antenna  system  is  designed  to  respond  only  to  the  cross- 
product  of  these  two  signals  which  return  after  striking  the  targets. 

Such  a  cross-correlation  radar  system  has  distinct  advantages  over 
the  conventional  system.  It  will  be  shown  that  the  directivity  for  a 
given  side  lobe  level  can  be  increased  by  extending  the  Dolph-Chebyshev 
synthesis  technique  to  the  composite  pattern  of  the  four  antennas.  Further¬ 
more,  since  the  system  is  designed  to  respond  to  the  cross-product  of  the 
two  returning  signals,  it  will  be  shown  that  the  effect  of  remote  active 
noise  sources  is  virtually  eliminated  in  the  time-averaged  system  output. 
9,1  The  Cross-Correlation  Radar  System 

A  diagram  of  the  system  is  shown  in  Figure  14.  Thete  are  four 
antennas,  two  of  which  are  used,  for  transmitting.  :(A^.  and  A^),  and 
two  of  which  are  used  for  receiving,  (B^and  B2).  The  pulsed  transmitter 
signal  is  split  in  two  equal  parts  and  that  part  which  is  fed  to  antenna 
A2  iS  ta&Sed  with  an  additional  frequency  u  i.e0,  the  carrier  frequency 

is  shifted  from  to  to  w  +  co  wi  th  oj  »>  tu  .  This  distinguishes  the  two 
o  o  1*  o  1 

signals  which  are  transmitted  by  the  antennas. 
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After  striking  the  target  (or  targets)  the  signals  return  and  are 

received  by/  the  other  pair  of  antennas,  B  and  B  In  practice  of 

course  the  four  antennas  need  not  be  physically  separated.  One  could 

use  four  coincident  antennas  which  share  the  same  aperture  (see 

Section  9.2).  The  cross-correla lion  of  the  terminal  voltages  of  these 

antennas  is,  in  principle,  the  same  as  with  any  correlation  system, 

except  that  the  various  filters  ^uSeid  touat  :be  wide-hand  .in  ■ 

* 

order  to  pass  the  video  pulses  of  the  radar.  A  detailed  analysis  of 
the  processing  is  given  in  Appendix  C. 

The  system  output  in  the  radar  case  is  usually  displayed  on  an 
oscilloscope.  Thus,  a  saw-tooth  wave,  synchronized  with  the  pulse 
repetition  rate,  is  used  to  produce  a  repeating  image  on  the  'scope. 

The  abscissa  of  the  display  is  usually  the  range  coordinate  s,  while 
the  vertical  deflection,  at  each  s,  is  proportional  to  the  system 
output  and  hence  to  the  strength  of  the  return  from  the  target  at  that 
range.  In  the  appendix  it  is  shown  that  for  a  target  at  the  range  r 
and  in  direction  v,  the  average  deflection  on  the  1  scope,  as  a  function 
of  antenna  beam  direction  u  and  range  coordinate  s,  is 

For  simplicity  of  analysis  we  have  assumed  in  what  follows  that  the 

filters1  frequency  characteristics  are  uniform  in  the  band 

|to-to  |<  £  to .  We  also  have  tacitly  assumed  the  case  of  t  =0, 
o 
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*  J*v4 

R(s  u)  =  P  (s-r)  M  (v-u)B  (v-u)A  (v-u)B  (v-u)e 

O  "IX  i  £>  £ 


J 


(184) 


where 


p  (s)  =  !  P(2  a/c)  |2 

O 


(185) 


is  equal  to  the  absolute  square  of  the  video  pulse,  P(2s/'c),  and  A^(u), 

A  (u),  B  (u),  and  B  (u)  are  the  field  strength  patterns  of  the  four 

U  1  Ci 

antennas . 

J2w  r/c 

The  phase  factor  e  i  <  depends^on  the  target  distance,  but  for  all 

targets  within  the  maximum  practical  range  r  % ,  one  can  choose  to  low 

max*  ’  1 

enough  so  that  this  factor  is  almost  unity.  Consequently,  we  can  rewrite 
Equation  (184),  to  a  good  degree  of  approximation,  as 


R(s;u)  = 


r 


P  (s-r 

o 


’t 


.  (v-u^  (v-u)A^  (v-u)B2  (v-u) 


(186) 


Now  with  a  conventional  radar  system,  where  quite  often  the  signals 

* 

are  transmitted  and  received  on  the  same  antenna  ,  the  analogous  output  is 


D(s,u)  =  P  (s-r)  |  D(v-u)  | 
7  o 


(187) 


* 

A  duplexer  is  used  to  isolate  the  transmitted  and  received  pulses. 
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where  D(u)  is  the  antennal  field  strength  pattern  and  we  have  assumed 
that  squafe-law  detection  is  used.  Comparing  Equations  (186)  and  (187) 
we  see  that  mathematically,  the  latter  is  just  a  degenerate  case  of  the 
former  which  occurs  when 

A  (u)  «  A  (u)  =  B  (u)  =  B  (u)  =  D (u)  -  (188) 

1  2  1 


Clearly  the  cross-correlation  system  is  preferrable  since  its  four 
distinct  patterns  can  be  specified  independently;  in  a  sense,  this  system 
has  four  times  as  many  "degrees  of  freedom"  as  the  conventional  one.  We 
will  now  consider  the  implications  of  this  in  the  design  of  optimum 
Dolph-Chebyshev  systems. 

9.2  Optimum  Dolpb-Chebyshev  Design 

We  recall  from  Section  7,3  that  the  pattern  of  an  optimum  Dolph- 
Chebyshev  array  is  given  by  Equation  (153..)  shown  below 

D  (40  =  T  (40  • 
n 

However,  when  such  an  array  is  used  for  radar,  the  system  output  is 
proportional  to 

D(s,40  =  P  (s-r)  T4  (4j  -4)  (189) 

*.  o  no 


where,  by  analogy  with  Equation  (187),  we  have  assumed  a  point  target  at 
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range  r  and  in  the  direction  4  .  As  a  function  of  4  “  *4  the  system 

o  o 

output  is  a  polynomial  of  degree  4n  but  is  not  the  desired  Chebyshev 
polynomial  of  that  degree^  (4  -  4). 

Now  let  us  consider  the  special  case  of  the  cross-correlation  radar 
system  when  the  four  antennas  A  A  B  and  Bn  are  in  the  form  of 
coincident  linear  arrays.  This  is  indicated  in  Figure  15  where  an 
aperture  weighting  network  is  shown  which  takes  the  signals  to  and 
from  the  elements  of  the  array  and  combines  them  in  a  predetermined  way 
to  produce  the  four  distinct  patterns  A^(4),  A^(4),  B^(4),  and  B^C^lO- 
In  particular^  if  we  let 


A  (4)  =  T  (4)  + 

1  n 


1 

2 


B  (4)  =  T  (4)  - 

a  n 


theiij  except  for  a  constant  factor^  we  can  show  that  the  correlation 


system  output  is 
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Figure  15.  The  Aperture  Weighting  Network  Associated  with  the  Four 
Coincident  Radar  Arrays. 
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R(s,40  -  P  (&-r)  T  (iji  -  +)  •  (191) 

7  o  4n  o 

This  is  the  desired  output  since  the  system* s  over-all  pattern  is  given 
by  the  Chebyshev  polynomial  of  degree  4n. 

A  quantitative  analysis  of  the  sidelobe  level  of  this  pattern  and 
of  that  of  the  conventional  case,  as  given  in  Equation  (189) ,  shows, 
that  if  both  have  the  same  beamwidth  measured  to  the  first  null,  then 
for  large  arrays  (e„gtJ  f  =  ,  n  >  9),  the  cross-correlation  patterns 

sidelobe  level  is  about  7  or  8  db  lower  than  that  of  the  conventional 
pattern.  This  is  shown  in  Figure  16  where  the  sidelobe  level  improvement 
in  db  is  plotted  as  a  function  of  the  conventional  pattern’s  sidelobe 
level,  also  in  db.  The  second  curve,  indicated  by  "Mattingly  System", 
is  the  improvement  obtained  when  two  distinct  antennas  are  used,  one  for 
transmitting  and  one  for  receiving.  If  the  two  patterqp  are  given  by 


a(4o  =  t  (4j)  -  — 

n  4a” 


(192) 


B(40  -  T  4*  +  — 
n  ^ 


(193) 


then  the  system  output,  corresponding  to  Equations  (189)  and  (191), 


is 
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Figure  16.  Sidelobe  Level  Improvement  of  the  Mattingly  and  Cross-Correlation 
Radar  Patterns  as  a  Function  of  the  Conventional  Dolpb-Chebyshev 
Sidelobe  Level. 
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R„(s,40  -  P  (s-r)  T*  (+  -  +)  (194) 

M  *  o  2n  o 

and  the  pattern  is  better  than  that  of  the  conventional  system  in  that  its 

sidelobes  are  lower  by  about  5  or  6  db  as  indicated  in  Figure  16.  This 

53 

system  was  first  proposed  by  R.  L,  Mattingly  in  1960 

It  can  be  concluded  that  the  conventional  radar  antenna  system 

employing  a  Dolph-Chebyshev  design  is  inadequate  due  to  the  multiplicative 

effect  of  the  transmitting  and  receiving  patterns  and  of  the  square-lav/ 

detection.  Thus,  instead  of  the  optimum  Chebyshev  polynomial  T  (40,  one 
4 

obtains  T^OjO.  It  was  shown  that  by  properly  selecting  the  four 
individual  patterns  of  the  cross-correlation  system,  (Equation  XL90)), 
one  can  s  ynthesize  the  over-all  pattern  T^(40  which  of  course  is 
optimum  in  the  sense  mentioned  above. 

9d3  Noise  Suppression 

Since  the  transmitted  and  hence  the  received  signals  are  in  the 
form  of  two  distinct  pulse  trains  whose  carrier  frequencies  differ  by  the  fixed 
amount  the  receiving  section  of  the  system  has  been  designed  to 

respond  to  the  cross-product  of  these  signals  (it  therefore  is  a  true 
cross-correlation  system). 

On  the  other  hand,  given  a  background  distribution  of  active  noise 
sources,  the  random  noise  due  to  this  distribution,  which  is  incident  on 
the  receiving  antennas,  will  in  general  not  have  such  a  dual  nature.  The 
noise  distribution  can  be  given  mathema tically  bv  N(t^u)  which  gives  the 
field  at  the  origin  of  the  antenna  coordinate  system  at  time  t  due  to  the 
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noise  source  in  the  direction  u.  The  terminal  voltages  of  the  two 

antennas, due  to  the  noise^can  be  written  as  N  (t)  and  N  (t) 

**1  2 

where 


(t) 


N(t,u)  B  (u)du 
*  k 


(195) 


for  k  =  1,  2. 

Since  the  noise  and  signals  are  uncorrelated,  the  time-averaged 
effect  of  cross-products  of  noise  and  signals  in  the  cross-correlator 
output  will  reduce  to  zero.  In  Appendix  C  it  is  shown  that  the  "self- 
product'*  of  the  noise  in  the  output  can  be  written  as 


- jco  t 

n  (t)  =  N  (t)  N*  (t)  e  '  1 
o  B,  B„ 


(196) 


The  average  value  of  n  (t)  sampled 

o 


2s 


for  N  samples  is 


N-l 

E 

k^O 


(—  +  k  t  )N* 
c  o  B2 


k  T 

O 


C 


(197) 


where  we  have  changed  variables  from  time  t  to  range  s, 


with 


This  is  proportional  to  the  average  vertical  deflection  of  the  oscilloscope 
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at  range  s  and  in  the  appendix  it  is  shown  that  the  expected  value  of 
this  deflection  is 


E 


sin  T 
_ 1  o 

N  sin  &  T 
1  o 


(198) 


where 


But 


and 


B1S2 


(199) 


E  <n  (s) 
o,N 


B1B2 


sin  N  to  T 
1  o 

N  sin  to  t 
1  o 


(200) 


1  lm 

N  — >  °o 


I 

i 

lim 

sin  N  to  t 

E  < 

L°'  J 

> 

^  N->°o  I 

V 

1  o 

1 

I 

N  sin  <o  t 

1  o 

=  0 


(201) 
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except  for  ujt  =  2q  IT  .  with  q  integral.  Consequently  the  time-averaged 
1  o  J 

output  due  to  the  noise  is  inversely  proportional  to  the  number  of  pulse 

repetitions  N,and  in  the  limit  as  N  approaches  infinity,  the  noise  output’s 

average  valpe  approaches  zero.  We  can  explain  this  by  noting  that  the 

system  is  designed  to  respond  to  a  pair  of  signals  with  a  fixed  frequency 

difference  .  Since  the  noise  does  not  have  this  dual  nature  its  effect 

in  the  system  output  averages  to  zero.  The  special  cases  u>  t  ~  2q  7 T  need 

1  o 

not  concern  us  since  we  can  always  design  our  system  to  avoid  these  values. 

This  result  contrasts  with  the  noise  output  of  a  conventional  system 
employing  square-law  detection.  In  such  a  situation  the  part  of  the 
terminal  voltage  of  the  antenna  which  is  due  to  the  noise  is 


v*> 


N(t,  u) 


D(u)  du 


(202) 


where  S(u)  is  the  antenna’s  field  strength  pattern.  The  low-frequency 
output  of  the  square-law  detector  due  to  this  noise  voltage  is 


“si/1’  *  i  I  Vl>  |2 


(203) 


This  output  has  a  positive  average  which  we  define  as 
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<nSLU)> 


N 

D 


(204) 


Consequently,  in  the  presence  of  background  noise  the  conventional  system^ 

display  will  have  a  positive  vertical  deflection  whose  average  value  is 

independent  of  range  and  is  proportional  to  the  average  noise  power  , 

However,  if  the  antenna  is  scanned  this  averaged  noise  output  will  vary 

with  scan  angle  (except  in  ’•he  rather  unrealistic  case  when  the  noise 

sources  are  uniformly  distributed).  On  the  other  hand, the  cross- 

correlation  system  output  due  to  the  noise  sources  has  a  time-average  of 

zero  which  in  general  is  independent  of  both  range  and  scan  angle. 

9P4  Suppression  of  Jamming  Signals 

In  some  situations  there  might  be_,  in  addition  to  the  target  echoes^ 

a  series  of  jamming  pulses  which  are  received  by  the  antenna.  They  will 

be  from  the  directions  of  the  jamming  transmitter  (or  transmitters)  and 

of  any  targets  which  can  reflect  these  j«.r.ming  pulse-  towards  the 

receiving  antenna.  If  the  carrier  frequency  ^  and  pulse  repetition  rate 

o 

t  of  this  interference  is  the  same  as  that  of  the  desired  echoes,  a 
o  * 

conventional  system  has  no  choice  but  to  accept  both  the  signal  and  the 
jamming  noise. 

The  cross-correlation  system  differs  from  this  since  it  is  designed 

to  respond  only  to  a  pair  of  target  echoes  at  frequencies  to  and  to  +  w  . 

o  o  1 

Therefore  unless  the  jamming  signals  are  also  in  this  dual  form  the  time- 
average  of  the  interference  will  approach  aero  as  fhe  number  of  pulse 
repetitions  increases  The  proof  of  this  statement  is  essentially  the  same 
as  thai  for  the  suppression  of  random  background  noise  of  Section  9  3. 
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10.  THE  MAPPING  OF  A  TARGET  DISTRIBUTION  WITH  A  CROSS -CORRELAT ION  RADAR  SYSTEM 

An  arbitrary  distribution  of  radar  targets  will  generally  produce 
echoes  which  are  partially  coherent.  In  this  chapter  a  general  theory  of 
mapping  such  a  target  distribution  is  presented.  One  must  use  a  cross¬ 
correlation  system  the  details  of  which  are  described  in  Chapter  9,  The 
resolution  of  the  system  is  limited  in  two  ways.  First,  the  non-zero 
beamwidths  of  the  antenna  patterns  smooth  the  observed  target  distribution, 
as  the  patterns  are  scanned .  Second,  due  to  the  transmitted  pulses  not 
being  infinitely  sharp,  the  measurement  of  the  range  of  the  targets  will 
also  be  a  smoothed  version  of  the  true  range  distribution. 

Let  us  consider  a  system  which  uses  just  one  aperture  of  length  L. 

Its  four  patterns, A  (u),  A (v) ,  B  (u)  and  B  (v)  are  generated  by  the  aperture 
weighting  network  of  Figure  15,  We  will  assume  that  the  radar  pulses  have 
a  spectrum  of  width  AP.  Under  these  restrictions,  an  aperture  of  length  L 
and  a  bandwidth  of  it  is  pertinent  to  ask  what  is  the  best  one  can  do 

in  mapping  ,the  target  distribution.  We  will  now  show  that  we  can  obtain  a 
principal  solution  which  is  similar  to  the  one  described  in  Chapter  5  for 
active  sources. 

10,1  Description  of  the  Target  Distribution 

The  targets  are  assumed  to  be  located  within  the  range  of  the  system 

2  . 

and  also  in  the  far-field  of  the  antennas  f2L  /\  r  <  r  )„  We  associate 

max 

with  the  target  in  the  u  direction  and  at  range  r,  the  complex  number 
j w  r/c 

g(t, u,  r)e  o  Note  that  it  is  a  function  of  time;  this  allows  for 

target  scintillation,  rotation  and  translation  as  time  progresses.  The 


Ill 


r  ■'  c 

phase  factor  e  is  present  because  our  reference  point  is  located  at 

jw  r/c 

the  center  of  antenna  aperture  Thus  g(t,u, r)e  represents  the  target 

distribution  as  observed  from  the  antenna. 

* 

10,2  System  Output  in  the  Presence  of  the  Target  Distribution 

Let  us  suppose  that  we  are  interested  in  the  amount  of  correlation  be¬ 
tween  the  echoes  from  the  targets  at  (u,i)  and  (v, r) .  If  we  beam  one  of 
the  transmitting  antenna  patterns  in  the  u  direction  and  the  other  in  the 
v  direction  and  if  we  delay  the  set  of  pulses  sent  to  the  nearer  target  by 
s/c  ~  (i-r).'c  seconds,  then  the  two  targets  will  be  struck  simultaneously 
by  the  pulses.  The  delay  of  the  pulses  sent  to  the  nearer  target  compen¬ 
sates  for  the  delay  imposed  by  the  extra  distance  s  on  the  pulses  sent  to 
the  other  target.  Since  the  returning  echoes  also  have  a  different  distance 
to  travel,  the  output  signals  will  also  have  a  relative  delay,  i.et?  the 
signals  from  the  more  distant  target  is  delayed  in  returning  by  an  addi¬ 
tional  s/c  seconds.  Accordingly  the  output  of  the  receiving  antenna  whose 
pattern  is  beamed  at  the  nearer  target  is  delayed  by  the  same  relative 
amount.  This  will  insure  that  the  echoes  from  the  two  targets  enter  the 
cross-correlator  at  the  same  time. 

The  two  transmitted  signals  are 


N-i  jtu  |/c 

A.n-I'c  u  -u)  =  Z  PCr.i/c  -  2kT  »  A  In  -ule  ° 

1  l  .  oil 

k-  o 


(205) 
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A2(t+r/c,  u-v) 


N-l 

2 

k=o 


j 

P(t+r/c  -2kT  )  A  (u,-v)e 
o  2  1 


r/c  + 


*] 


(206) 


We  have  advanced  one  signal  by  r/c  seconds  and  the  other  by  i/c  seconds* 

The  second  signal  is  tagged  with  the  constant  phase  modulation  .  Note 

that  the  pulse  repetition  rate  is  2t  where  t  =  2r  /c.  This  permits  a 

o  o  max 

delay  of  one  pulse  relative  to  the  other  of  as  much  as  t  with  practically 

o 

no  interference  from  preceding  or  succeeding  pulses. 

After  striking  the  targets  the  pulses  return  and  are  received  by 

antennas  B  and  C  Their  two  output  voltages  are 
1  £ 


2w 

oo  oo  N-l  I  2(i-r,  )  j — (f-r  ) 

le'  C  A  (u  -u)  + 

_  _  k=e  ’ 

00-00 


00  oo  N-l  2(i-r  ) 

B  (t  i  r  u  v)  =  I  f  2  J  P(t-: - - 2kr  )i 

1  J  j  k=0  1 


-  %i.  2kTo)eJ 


CO 


i+r-2r  *1 

—  (i+r-2r  )  +  w  (t.+  - - -) 

c  1  1  c  J 


a2  (Vv) 


i-r. 


2co  r 

3  o  1 


B1(u1-u)  g(t  +  r  )  e  c  du  dr 


(207) 
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B2(t,i,r,u,v)  = 


OO 

/  / 


-  GO  -  OO 


oo  N-l 
2 
k=o 


fP(t+ 


0  +  r 


2r 


-2kT  )e 
o 


o 

7T 


«4r-2i2) 


Al(u2-u)  + 


2  (r-r  ) 
P(t+; - — 


j 

-  2kT  )e 
o 


g 


200  r 
J  o  2 


e  c 


du  dr 
2 


2  * 


(208) 


Note  that  the  voltage  (t^  Jt}  rf  u,  v)  has  been  advanced  by  the  amount  i/c 
and  Bg  (t,  i  9  r,  u}  v)  has  been  advanced  by  r/c.  In  Appendix  D  it  is  shown 
that  when  these  signals  are  cross-correlated  with  t=0,  the  system  output 
can  be  written  as 


-  G(o,f,r,u,v)  *  C^fu)! 


(v)  P(i ) P(r) 


(209) 


where 


C^u)  =  Ai(u)B1(u) 


(210) 


C2(V)  =  A2(V)B2(V) 


(211) 


P(r)  =  P(2r/'c) 


(212) 
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and 

G(t,  f,r,  u,  V)  =  </g(t,  u,f)  g  (t-T,  v,  r^>  (213) 

is  the  mutual  coherence  function  of  the  target  distribution. 

Equations  (210)  and  (211)  show  that  in  effect  we  have  two  independent 
radar  systems  each  of  which  uses  distinct  patterns  for  transmitting  and 
receiving.  One  system  beams  its  two-way  pattern  in  the  u  direction  and 
the  other  system  beams  its  pattern  in  the  v  direction.  Equation  (212) 
indicates  the  role  of  the  pulse  shape  in  the  mapping  of  the  target  distri¬ 
bution.  P(r-r  )  =  p[—  (r-r  )  is  the  response,  as  a  function  of  range  r, 
to  a  point  target  at  r^.  The  pulses  P(t)  in  the  radar  system  correspond 

to  the  RF  filter  response  functions  f  (t)  and  f  (t)  in  the  cross- 

A  B 

correlation  receiving  system  described  in  Chapter  5. 

10.3  Fourier  Analysis  of  the  System  Output 

The  system  output  R(i,r,  u, v)  is  a  smoothed  version  of  the  mutual 
coherence  function  of  the  target  distribution  (evaluated  at  t=0) .  If  we 
take  the  4-dimensional  inverse  Fourier  transform  of  Equation  (209),  the 
convolution  theorem  can  be  used  to  obtain 


r<M>x,y)  =  (2*)2  g(o,S,n,x,y)  p(&>  ^*(n)c(x).  c*(y^>  (214) 


where 


,  g(o,S,ti,x,y)#  p(£),p<n),  djCx), 


and  c2  (y)  are  the 
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inverse  Fourier  transforms  of  R(£,  r,  u,  v),  G(o,  i,  r,  u,  v),  P(r),  C^(u), 

and  C_(v)  respectively.  Now  since  in  Equation  (210) 

A 

C1(U)  =  A1(U)  B1(U) 


we  can  again  invoke  the  convolution  theorem  to  obtain 


6i 


(X) 


a^jc)  *  bx(x) 


(215) 


where  a^ (x)  and  (x)  are  the  inverse  Fourier  transforms  of  A^(u)  and 
(u)  respectively.  In  a  similar  fashion  we  can  write 


c2(y) 


aij(y)  *  b2(y). 


(216) 


c  (x)  and  c:  (y)  can  be  thought  of  as  the  spatial  frequency  spectra  of  the 

X  A 

radar  system’s  composite  patterns  C  (u)  and  C  (v),  respectively, 

X  A 

Since  the  four  apertures  coincide  and  are  all  of  length  hf  it  follows 

that 


c^(x)3  0^  for  |  x |  >  ly 


c2(y)  £  0,  for  |y|  >  L  . 


(217) 
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Thus  the  system  has  a  spatial  frequency  response  which  is  non-zero  only 
within  the  region  of  the  xy  plane  for  which  |x|  ^  L  and  |y|  ^  L, 

Similarly  we  define  p{£)  as  the  spatial  frequency  spectrum  of  the 
system's  response  in  the  radial  direction  t t  This  spectrum  is  related 
to  the  temporal  spectrum  of  the  transmitted  signal  pulse.  Thus  if 


(218) 


then 


(219) 


and 

5  «>  -  f  p  <§e>  . 


(220) 


But  we  stipulated  that  p(c*>)  was  practically  zero  for  |io|  >  (see  page  110); 
consequently  we  have 


p(£:>  « 


for  |£(>  ML  . 


c 


(221) 
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Thus  the  antenna  system  will  reject  those  radial  spatial  frequencies  which 

e  2AW 

lie  outside  the  interval  |5|^  — —  .  Clearly  it  will  also  reject  those 

frequencies  ^  such  that  „ 

It  has  been  shown,,  therefore,  that  the  correlation  antenna  system  acts 
as  a  low-pass  filter  of  the  spatial  frequencies  x  and  y  which  correspond  to 
the  bearing  directions  u  and  v  and  it  acts  as  o  low-pass  filter  of  the 
spatial  frequencies  £  and  'H  which  correspond  to  range  coordinates  i  and  r. 
The  "aperture11  of  the  system  is  a  four-dimensional  region  centered  at  the 

C  +  2A^ 

origin  and  with  boundaries  given  by  the  hyperplanes  5  =  -  j 

+  +  4* 

^  =  _  — —  ,  x  =  _  L,  y  =  _  L.  Furthermore,  within  this  region  the  system 
function  weights  the  target  spectrum  in  a  more  or  less  arbitrary  fashion 
which  depends  on  the  pattern  and  pulse  spectra.  The  output  spectrum  is 
usually  a  distorted  version  of  the  target  spectrum  because  of  this. 

i 

However  we  can  use  Equation  (214)  to  define 

g  (Oj^n^y)  =  - -5 — ,r,(K!L*ir> - — -  r  fnr  (222) 

°  <W)2d1<*)«>*(y)  P(6)  ?<1) 

=  0,  otherwise, 

where  we  have  defined  the  four-dimensional  point  set  Q  to  be  that  for  which 

Ic  (x)  c*(y)  p(£)  p*  (*1)  I  £  0.  In  our  case  it  is  the  above-mentioned  four- 
1  2 

dimensional  volume  with  center  at  the  origin.  The  Fourier  transform  of 
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g^Co,  £,  'H ,  x,  y)  will  be  defined  as  the  principal  solution, Go<o,^,r,u,v)  . 

The  spatial  frequency  spectrum,  g^to, 5, T, x, y) ,  is  identical  to  that  of  the 
required  mutual  coherence  function  G(o,-^  ,r,u,  v),  but  only  over  the  region  Q. 
For  points  outside  of  Q  it  is,  by  definition,  equal  to  zero.  Consequently 
G  (o,^,r,u,v)  is  a  smoothed  version  of  G(o,^,r,u,v)  whose  non-zero  spatial 
frequency^  components  are  identical  to  those  of  G(d>,^,r,u  v)  ,  In  this  sense 
then,  G  (o,-^,r,u,v)  is  the  best  or  principal  solution  to  the  problem  of 
mapping  a  stationary  target  distribution  with  a  finite  antenna  and  a  finite 
bandwidth* 
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11.  CONCLUSIONS 

In  this  thee  is  we  have  considered  the  application  of  cross-correlation 
techniques  to  linear  antenna  arrays.  Probably  the  most  importan'  result  of 
this  work  was  the  demonstration  that  a  two-antenna  cross-correlation  system 
can  be  used  to  measure  a  source  distribution's  mutual  coherence  function. 
Since  coherent  and  partially  coherent  signals  are  not  uncommon  in  practice, 
©„g.,  radar,  multipath  propagation,  and  scattering  from  nearly  objects,  the 
results  obtained  here  should  be  of  interest,  to  the  antenna  engineer.  The 
demonstration  tl$at  a  conventional  antenna  system  with  square-law  detection 
is  just  a  special  case  of  the  cross -correlation  system  makes  it  possible  to 
evaluate  in  a  straightforward  manner  the  relative  merits  of  a  conventional 
and  a  cross-correlation  system.  An  interesting  result  is  that  the  analysis 
of  multiple  correlation  systems  is  quite  simple  Cat  least  in  theory)  since 
their  outputs  are  ammenable  tc  Fourier  analysis  just  as  in  the  two-antenna 
case.  Finally, it  is  felt  that  the  new  cross-correlation  radar  system,  in 
spite  of  its  involved  circuitry,  has  a  number  of  advantages  over  the  con¬ 
ventional  system  which  should  recommend  i t  in  some  applications. 
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APPENDIX  A 

THE  CASE  OF  TWO-DIMENSIONAL  SOURCE  DISTRIBUTIONS  AND  PLANAR  ANTENNAS 

In  the  main  body  of  this  thesis  we  considered  a  one-dimensional  source 
distribution.  More  realistically,  the  direction  of  a  remote  radio  source  can 
be  specified  as  a  point  on  a  sphere  of  "infinite"  radius  with  the  antenna  at 
its  center.  This  leads  to  a  two-dimensional  source  distribution  which  can  be 
measured  by  a  planar  antenna  system.  We  will  now  outline  this  generalization 
from  one  to  two  dimensions.  At  the  same  time  we  will  consider  the  practical 
problem  resulting  from  the  directive  properties  of  the  elements  of  the  array. 
Then  the  problem  of  changing  variables  from  the  (u, v)  domain  to  the  (0,  (0  do¬ 
main  will  be  studied.  Finally  we  will  show  that  for  electromagnetic  fields  the 
Fourier  transform  relation  between  the  far-field  pattern  and  the  aperture 
distribution  is  a  valid  relation  when  the  far  electric  field  and  the  aperture 
current  density  are  considered.  The  relation  between  the  far  electric  field 
and  the  aperture  electric  field  is  not  as  simple. 

A . 1  The  Output  of  a  Planar  Antenna  Cross-Correlation  System  and  the  Effect  of 

Element  Pattern 

Let  us  consider  a  planar  antenna  located  in  the  z  =  0  plane.  The  total 
pattern  of  the  antenna  can  be  written  as 


=  E(u  v)  A  (u  v) 
>  a  ’ 


A(u,  v) 


(A,  1 ) 


where 


u  =(3sin0  cos  (p  v  =  |3sin0  sin^ 


(A. 2) 


and  E(UjV)  is  the  element  pattern  which  as  a  function  of  u  and  v  varies  much 

more  slowly  than  the  array  pattern  A  (u,v).  The  terminal  voltage  of  the  array 

a 

in  the  presence  of  the  source  distribution  V(t^u^v)^  is 


ACt,  u.  v) 


V(t,  u *  v *  ) 


it  it  it 

E(u  v  )  A  (u  -u,v  -v)  du  dv. 
3  a  3 


(A.  3) 


Note  that  the  array  pattern  can  be  scanned  but  the  individual  elements*  pattern 
cannot  be  scanned. 

If  we  let 


VgCt^u^v)  =  V(t_,u,v)  E (u}v) 


(A.  4) 


then  Equation  (A. 3)  becomes 
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A{t.U.v)  = 


ff 

JJ 


V,(t.u  v  )  A  (u  -u.  v  -v)  du  dv 
£  -  *  a 


(A. 5) 


=  V  (t  u.  v)  *  A  (u.  v) 

ci  a 


(A.6) 


From  Equation  (A, 6)  we  see  that  the  effect  of  the  element  pattern  is  to  change 
i;he  source  distribution  as  "observed"  by  the  array.  This  change  is  given  by 
Equation  (A, 2).  It  can  be  easily  shown  that  the  mutual  coherence  functions*;/ of 
the  "observed"  source  distribution  is 


:tr<Vi«V>)  "  <^v*(t'ui>vi)  VE(t-T>Vv2^> 


(A. 7) 


=  E(ui>Vl)  E*(u2,v2)  T(x,upv;I 


VV 


(A. 8) 


where  T(t? u  f v  , u  v  )  is  the  true  mutual  coherence  function. 

X  i*  jL 

The  output  from  a  second  planar  array  can  be  written  as 


B(t,  u,v) 


V  (t, u, v)  *  B  (u.v) 
E»  a  > 


(A.  53) 
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and  when  the  two  output  voltages  are  filtered  and  cross-correlated  we  obtain 


HC^VVVV  =  VT'ur  WV 


*  fAB(T)  ®*CVV“ 


(A. 10) 


Equation  (A. 10)  indicates  that  the  cross-correlated  output  of  the  system  is  a 
filtered  version  of  T  (t;u  v  ,u  )  and  not  of  the  true  mutual  coherence 

J-*  i.  -L  ^5  & 

function  T<-r,u1>v1,u2,v2'>' 

The  inverse  Fourier  transform  of  R(t u, , v, , u„, v  )  is 

’  V  1-’  2’  2 


r((ji’Xl’yl’X2’y2)  =  tE(fa>^Xl'Wy2)  CAB(^VyrVy2>  (A.  11) 

where 

=  FACW)  FB(W)  a(xpyi}  b+(Vy2>  CA-12) 


We  define 


r(^Xl’y1’Vy2) 


AB 


£  QAB  (a . i 3 ) 


0 


otherwise. 
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The  Fourier  transform  of  t  (to,  x  y  x  ,  y  )  is  the  principal  solution.  But 

X  I7  £t 

in  view  of  Equation  (A,8)  we  can  write  the  principal  solution  as 


(A. 14) 


TEAB<T'VVVV  =t(t^upWv2)  e(VV  e  (VV  * 


where^  by  analogy  with  Equation  (117)^  Q(t^  u  y  v  }  u  ,  v  )  is  the  Fourier  trans- 
form  of  a  system  function  which  is  uniform  within  Q  j  the  system's  "aperture". 

Note  that  the  element  patterns  are  included  in  the  convolution  integral  of 
Equation  (A.14)  and  strictly  speaking  cannot  be  separated  from  the  desired 
mutual  coherence  function  TCTjU  v  u  v  ).  However^  in  most  cases  the  antenna 
elements  are  much  smaller  than  the  array  itself  and  the  element  patterns  are 
slowly  varying  functions  when  compared  to  QCt^u  v  u  ,v  ),  If  this  is  the 

-L  X  £§  & 

case  we  can  write 


(A. 15) 


TEAB(T'VV1'VV2)  =  E(Vvi}  E  (u2'V  \  t(t'Vvi-Vv2)  *  q(t>WVv2) 


7 


to  a  good  degree  of  approximation.  Consequently  we  obtain 
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T  (t  u  v  u  v  )  = 

ABV  ’  V  1>  2'  2 


teab(t<  VVVV 

E<VV  E*<VV 


(A. 16) 


as  the  principal  solution  corrected  for  the  effect  of  element  pattern.  This 

result  is  valid  only  in  regions  where  Je(u^  v)J  is  non-zero. 

A, 2  The  Mutual  Coherence  Function  in  the  Angular  Domain 

In  practice  it  is  more  meaningful  to  express  the  mutual  coherence  function 

as  a  function  of  the  angular  variables  0  and  rather  than  of  u  =  psin0  cos cp 

and  v  =  Psin0  sin(p.  Since  T(r^u  ,v  ,u  ,v  )  is  a  power  density  function  we 

Know  that  T(Ofu  ,v  }u  v  )  du  dv  du  dv  is  the  cross-power  incident  on  the 
X  X  a  a  X  X  2  2 

antenna  system  from  the  elementary  regions  located  at  (u  . v  )  and  (u  }  v  ) . 

x  1  2  2 

Transforming  this  element  of  cross-power  from  the  (u;  v)  to  the  (8^^)  domain^ 
we  write 


T<°)'V  WV  dV''ldU2dVi 


T  (Oj 


P  s  i  n0  ^  c  os  <p^ ,  P  s  i  n0  ^  s  i  n(p^  ?  p  s  i  n0  ^  c  o  s(p^ ,  p  s  i  nQ 


s 


inV 


SCUj, 

5^7 


WV 


d9  d<£  dQ  d<p 

1  2  r2 


(A, 17) 
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where 


VvvV 

8«i 


=  P  cos0coe0osin0  sin0r 


(A.  18) 


is  the  Jacobian  of  the  transformation. 

Now  let  us  define  the  mutual  coherence  function  in  the  angular  domain 

t 

to  be  T(Tj0,  9(p  9  •  The  cross-power  incident  on  the  system  from  the 

JL  J-  & 

elementary  regions  located  at  (0  f  cp. '  )  and  (0  )  is 

-L  1  ^  2 


T  (0,  ei^<Pl’92^2)  sin91sin92  d(pl  d<^2  Cl0l  d02 


(A. 19) 


If  we  equate  this  with  the  cross-power  given  by  Equation  (A. 17)  we  obtain 


T  (0,  0^,^,  co&01.cos02,  T(O_tPsin0icos^1?Psin01sin<p1)Psin02cos^2. 


psin02sin^2> 


(A.  20) 


as  the  relation  between  the  mutual  coherence  functions  in  the  (u^v)  and 
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(Q3<p)  domains. 

This  result  also  applies  when  we  consider  the  source  distribution  to 
be  weighted  by  the  antenna  element  pattern. 


TE(T,01’^1j02^2)  =  P  c°sS1cos02  E  J  psin0icos(plipsineisint/>i 


E  <  psinG^cros^^psinQ^sin <p^  >  T(t^ psinG^os^P sinG^^sin^ <,Ps inQ^cos^, 


psinQ  sin<#  ) 

M  M 


(A. 21) 


and  if  the  element  patterns  are  slowly  varying  we  can  write 


T  (r  4  G  rfl  0  (p  )  - 
AB  v  ^°1'T1,D2-  2 


teaB 

EfPsinS^cos^^psinO^sin^  ) 


(A.  22) 


£  ' 

E  (P s in02 cosr/^ 9 p s ln^p  ) 


A. 3  Fourier  Analysis  for  Electromagnetic  Fields 

Let  us  assume  +hat  the  electromagnetic  field  is  known  on  the  plane 
z  -  o„  This  plane  can  be  thought  of  as  the  extended  aperture  of  an  antenna. 
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Any  electromagnetic  field  in  the  plane  can  be  decomposed  as  the  sum  of  two 
modes,  the  TE  and  TM  modes,  each  of  which  can  be  represented  by  a  scalar 
function.  In  particular  let  us  consider  the  TE  mode  and  a  field  which 
varies  only  in  the  x  direction.  This  is  analogous  to  the  one -dimensional 
apertures  which  were  considered  in  the  thesis.  The  electric  component  is 


E(x,  z) 


(x,  z)  . 


(A. 23) 


The  scalar  field  E  (x, o)  in 

y  ’ 

A*/ 


the  aperture  plane  has  as  Fourier  transform 


Ey(psin0),  which  represents  a  spectrum  of  plane  waves.  Considering  just 
one  of  these  plane  waves  we  can  write 


A  'v 

y  E  (Psin0) 

y 


e 


w  -j[(3sine  x  +  Pcos0  z] 

y  E^([3sin0)  e 


(A. 24) 


which  represents  the  field  at  any  point  (x,  z)  due  to  the  plane  wave.  The 
total  field  is  the  integral  of  all  such  plane  waves  over  all  possible 


directions „ 
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E  (x  z )  - 

y 


QO 


^  -  j[  psin0  x  +  (3 cos0  z] 

E  (p  s  i  n0  )  e 

y 


dp  s  i  n0  , 


(A , 25 ) 


The  field  at  any  point  can  be  recovered  from  its  plane  wave  expansion  on 
the  plane  z  =  o. 

Equivalent  Current  Sheet 

If  we  again  consider  the  plane  wave  given  by  Equation  (A, 24)  we  can  use 
Maxwell’s  Equation 


H(x3z) 


;url  E(x}z) 


(A. 26) 


to  write 


H^Cpsine) 


EqSsinG) 


fA.27) 


as  the  x  component  of  the  plane  wave’ s  magnetic  field.  But  it  also  follows 
from  Maxwell- s  Equations  that 


H  (x  o)  -  ~  J  (x) 
x  2  y 


(A. 28) 
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where  J  (x)  is  the  density  function  of  an  equivalent  sheet  of  electric 

y 

current  flowing  ir  the  z  =  o  plane.  Taking  Fourier  transforms  of  both 
sides  of  Equation  (A. 28)  gives 


H  ((3 s in0 )  si  J  ((3sinQ)  (A. 29) 

x  2  y 


and  in  view  of  Equation  (A. 27)  we  can  write 


J  (psinQ) 

y 


-goose  -Scpsme). 


(A. 30) 


Substituting  this  result  into  Equation  (A. 25)  gives 


E  (x,z) 

y  ' 


oo 


J  (PsinB)  e'JtPsl"9  51  +  I3*®0 

y 


(A. 31) 


which  represents  the  electric  field  at  any  point  in  terms  of  the  spectrum 
of  the  current  density  in  the  aperture  plane. 
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Far-Field  Approximation 

Let  us  rewrite  Equation  (A, 31)  as 


r  -j(3r  cos?  (0-0  ) 

E  (r  0 )  =  V  I  J  (|3sin0)  e  °  d0  (A.  32) 

y  o’  °  — —  /  y 

-oO 


where 


If  the  observation  point  (r  ,0  )  is  in  the  far-field  r  is  very  large  and  the 

o’  o  o 

integral  on  the  right  side  of  Equation  (A. 32)  can  be  evaluated  by  the  method 
of  stationary  phase.  We  can  write 


E  (r  0  )  Od  -  J(V 
y  o’  o  — — 


J  (3  s  in0  ) 
y  o 


2 


(6  -  e9)2 

*  d0 


nj 


IT 


H 


(2) 


(|3r  ) 
o 


J  (Psin0  ) 

y  o 


2 


o 


(A  33) 
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where  H ^ ^  (3 r  )  is  the  Haukel  function  of  the  second  kind.  Thus  we  see 

o  o 

that  the  far  electric  field  at  a  fixed  radius  r  from  the  origin  is  propor- 
- - - - — - —  o 

tional  to  the  Fourier  transform  of  the  current  density  function  in  the 
aperture  with  the  transform  variables  being  [3sm0  and  x.  This  is  not  the 
case  for  the  relation  between  the  aperture  electric  field  and  the  far 
electric  field.  Starting  with  Equation  t'A  2d;  we  can  show  that  in  the  far 
field 


f  o  i  ^ 

E  fr  a0  )  IT  (By  )  BcoS0  E  C'PsinQ  >  .  (A. 34) 

y  o>  o  o  ^  o  r  o  y  o 


Note  the  additional  factor  cos9 

o 
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APPENDIX  B 

OUTPUT  OF  THE  MULTIPLE  CROSS -CORRELATION  SYSTEM 


The  signal  from  each  of  the  N  antennas  is  passed  through  an  RF  filter 
and  delayed  a  certain  amount „  Thus  the  signal  from  the  kth  antenna  after 
the  filtering  and  the  delay  is 


-JW 


V'-W  "  V'-W 


O  k 


(Boi) 


where 


8  DO 

V*’V  -  J  /  v<t-t1.“)  \<u-ukHu  fk(tl)dt1 

-B  o 

=  V(t,uk)  *  fk(t)  Ak(uk)  . 


(Bv2) 


As  is  shown  in  Figure  13,  the  signal  from  the  first  antenna  is  fed  into  a 

frequency  shifter  (phase  modulator)  which  increases  the  carrier  frequency 

from  u  to  (j  +  (J  with  u  »>  u ,  -  Thus  the  signal  can  be  written  as 
o  o  1  o  1 

-j(w  VtJ.t) 

Vl(t~VW“  Ai(t“VUi)e  *  (B-3) 


j 

I 

I 


This  signal  is  fed  into  a  mixer  (a  squaring  device)  along  with  the  signal 

from  antenna  2,  V  (t-T  , u  )a  ^here  is  a  component  in  the  mixer  output 
2  2  2 

at  frequency  u  and  it  can  be  written  as 


J’i2<t'Ti,VvV=  V'-VV  V* 


r  ,  ^'VY 

-T  ill  e 
2'  2 


J^t 


(J3  *4) 


m 
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In  a  similar  fashion  the  signal  from  antenna  3  is  passed  through  a  frequency 
shifting  device  with  modulation  frequency  and  after  it  is  mixed  with 
the  signal  from  antenna  4  the  mixer  output  contains  a  -rignal  whose  carrier 
is  at  frequency  u i  <«  cj^  It  is  given  by 


A  w  C r  )  JU  t 

Vt'Yl'W  =  V'YY  V^W  "  3  1  e 


(Bv5) 


As  can  be  seen  in  Figure  13  the  signals  of  each  adjacent  pair  of  antennas  are 

correlated  in  this  manner  and  each  mixer  output  has  a  distinct  carrier 

frequency  u  }  with  k  =  1 , 3,5^ . * c 9 N-l * 

K 

* 

If,  after  passing  through  a  band-pass  filter^  each  of  these  single 
correlation  outputs  is  mixed  with  the  adjacent  one,  i,e5!  p  with  p  , 
p^_6  with  Pyg  i  etc.,  and  the  process  continued  until  a  single  output  results, 
we  can  write  the  expression  for  the  signal  in  the  output  at  frequency 
CJ  u  +  4*  GJ  as  follows 

1  J  J  ll**  i. 


123 _ N 


<V 


N’W 


-V 


(B„6) 


.  T  ,  r  .  ■w'i'VVV-  ■  -  T«w W’  ■**.1!?] . 
Vi"  K-i'Vi1  \"-YVe  1  I 


Then  if  the  output  is  fed  inTo  a  synchronous  de:ector  along  with  a  reference 
signal  at  frequency  ^ ,  the  detector  output  before  time-averaging 

will  be  the  following  real  voltage 


*For  simplicity  of  analysis  we  assume  that  the  filter  has  no  effect  on  the 
output  signals  in  the  narrow  band  ^  A  (J  and  completely  rejects  the 

double  frequency  RF  output  signals  and  the  signals  m  the  low  frequency 
band  t CJ j  ^  A  <j. 
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V2...*„(t’  'rVV-'V  V  V-V  -R.  A2(t-T2-V 


-i(J  (T  -T  +T  - .  .  .  -T  )) 
*  *  J  e>  1  ?T  T  NJ‘ 

A4(t-V  ^-V^V  .,0  12  3 


(B-j7) 


If  one  of  the  signals^  say  A^(t-TjjU^)e  f  is  taken  in  phase  quadrature 

then  the  above  result,  becomes 


R 


A  A  W-’-VV-'-V  =  zm  A2(t'W’* 

i-  M  IN 


•"•AN(t“V  v  e 


■juo[VVV-“V 


(Bv8) 


Thus  we  define 


R. 


1A2"  ‘ J TNj U1  J  '  "*UJ?  ~ 


l^Al‘'AN(t’Tl*  *’,Un^  ^  ^  ^Ai*  'ANCt,Tl'  *  '  ‘,UN^1 


(Bv9) 


as  the  complex  system  output  for  the  case  of  N  correlations.  Using 
Equation  (B.2)  we  can  write  the  expected  value  of  the  time-average  of  this 
output  as 


6  8  8  oo  oo  oo 


R 


L1A2**AN(T1’"‘,TN,U1 - 'V  "  E]/  /  ’ ’  /  /  /••*  *  /  V(t-T1“W 

tB  -6  -BOO  0' 

Vu-WV  v<t-VVu3> 


f  T  t 


1 


Al(UiUl)  W  f2(t2>  •-Wda1dU2”dUSdtl”dtN 


(B-10) 
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B  oo  oo 


■  / . /  /••/  T<N>,:V*rVt2 . VV  “i*  “a 

-B  -BOO 


. » \ 

>  *  •  *  >  ■'vs' 


V  V* 

\<V“I>  Vv 


'J* 

-u*) . , . .A  (u  -u*) 
2  N  N  N 


(B.H) 


VV  VV 


.f  (t  )du* . . .du'  dt. 
N  N  1  N 


.  .  dt 


N 


where 


"  •  >  UN> 


(V<t-T1>Ul) 


(t-T2Ju2) 


)>  (B-'12) 


Since  the  source  distribution  is  stationary^  this  function  is  independent 
of  the  exact  time  t  at  which  the  evaiuation  occurs.  Inspection  of  the 
integrand  of  Equation  CB-11)  will  reveal  that  we  can  write 


V2;"V<Ti'v"v'Vv'"u») 

\j  \/  V  *  V  V  *  V  y/  * 

.  [VWWV-'-'W  W  W . W1  <B-13> 


CN) 

which  is  the  convolution  of  T  (T^  ,  .  .  .  ,TN,  ,  u^)  by^tlie  reversed  the  system 

function  W  A2*(u2)...A;(„h)  VV'^VVV . 4<V  ' 

In  the  above  derivation  we  have  assumed  that  N  was  even.  In  most 
practical  cases  this  will  be  the  only  pertinent  situation  since  most  of 
the  source  distributions  will  have  even  probability  distribut ionsJ  e.g.., 
gauss ian,  and  hence  all  odd  moments  will  be  zero.  However,  this  will  not 
always  be  the  case  and  in  order  to  measure  the  odd  moments  ol  a  distribution 


we  must  use  an  odd  number  of  antennas  The  process  is  the  same  as  for  an 
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even  number  except  that  the  last  antenna  signal  is  shifted  in  frequency  by 

U)  -u  rather  than  by  u  •  It  is  then  mixed  directly  with  the  output 
No  N 

signal  resulting  from  the  mixing  of  the  signals  from  antennas  N-2  and 
N-l.  The  rest  of  the  process  is  the  same  as  before. 

Finally  we  note  that  the  two-antenna  system  is  just  a  special  case 
of  the  above.  The  output  of  the  single  mixer  is  fed  immediately  into 
the  synchronous  detector  along  with  a  signal  at  frequency  The  output 

iS  RA  A  and  if  We  let  Tl=  T2  =  T’  U1  =  U'-  u2  =  v>  A!  =  A> 

1  2 

and  we  have 


Bvl  (V  WV  “  RAB<T'U’V) 


(B-.14) 
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APPENDIX  C 

CORRELATION  RADAR  SYSTEM 

A  diagram  of  the  system  is  shown  in  Figure  14 .  The  pulsed  trans¬ 
mitter  signal  is  divided  into  two  equal  parts,  one  of  which  is  shifted 

in  frequency  from  w  to  w  +  oj  ,  with  w  »>  co  .  This  slight  change 
o  o  1  7  o  1 

in  frequency  can  be  accomplished  by  inserting  dir  the  feed  .line  a  rotary 
phase  shifter  ■  whose  phase  change  is  at  the  constant  rate  of  radians 
per  second. 

These  two  distinct  signals  are  then  fed  to  antennas  and  and 
the  far  fields  of  the  two  antennas  can  be  written  as 


N-l 

A  (t  v)  s*  S  P(t-k  T  )A  (v)  (C.  1) 

l  7  .  o  1 

k=o 


N-l  jco  t 

A2(t*V’<°l)  =  2  P(t"k  1fo>  A2(v)e  (C.2) 

k=o 


where  P(t)  is  the  video  pulse,  and  is  the  pulse  repetition  period.  We 
have  assumed  that  the  transmitter  is  turned  on  at  t  =  0  and  the  pulse  is 
repeated  N-l  times.  A_L  (v)  and  Ao(v)  are  the  far-field  patterns  of 

and  respectively.  The  absence  of  RF  filters  in  the 


a  n  te  nn  as 
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feed  lines  of  the  transmitting  and  receiving  antennas  of  Figure  14  is  due 
to  our  assuming^  for  simplicity^  that  the  actual  filters  have  uniform 
pass-bands  and  hence  pass  the  signals  unaltered. 

Now  for  simplicity  let  us  consider  the  case  of  a  single  target  in 
the  v  direction  and  at  a  range  r.  The  returning  echoes  are  received  by 
the  other  two  antennas^  B  and  .  The  terminal  voltage  of  antenna  B 
can  be  written  as 

CC.3) 

N-X  9  f  (t-  ~i  ju  (t-  — ) 

b  (t)  -  S  -  k  t  )|A  (v)  +  A  (v)e  JB  (v)e  °  c 

1  col  2  1 

k=o  W 


Note  the  time  delay  term  —  }  where  c  is  the  velocity  of  light  in  the 

medium*  This  is  the  time  required  for  the  pulses  to  reach  the  target  and 

after  reflection  to  return  to  the  antenna.  It  is  common  practice  to 

choose  the  pulse  repetition  period  so  that  if  the  maximum  practical 

2r 

target  distance  is  r  ,  then  t  =  max  .  Thus  all  echoes  due  to  one 

max'  o  - 

c 

pulse  will  have  returned  before  the  succeeding  pulse  is  transmitted.  This 
avoids  any  ambiguity  in  range  measurement.  If  for  the  present  we  assume 
that  the  target  is  stationary,  each  pulse  will  return  at  the  same  time 
relative  to  its  corresponding  transmitted  pulse.  Consequently,  we  can, 
for  simplicity  of  notation,  consider  the  time  interval  0  ^  t  <  knowing 
that  the  returns  of  succeeding  intervals  will  be  identical  to  the  first 
whose  response,  with  the  carrier  frequency  factor  suppressed,  can  be 


written  as 
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B1(t)  =  P(t-tr) 


(v)  + 


A  (v)e 

a 


ja)  (t-t  ) 
a  1  r 


] 


(v)e 


-ju  t 
or 


(C.4) 


for 


OX  t  <  r 

^  ^  o 


) 


and  where  we  have  let 


(C.  5) 


The  voltage  of  antenna  is 


Vl> 


-  Pft-tj;) 


(v)  + 


A  (v)e 
2 


•n 

r 


(y)e 


•Junt 
o  r 


(C,  6) 


for  0^  t  <  tq  y  where  B^Cv)  is  the  receiving  field  strength  pattern 
of  antenna 

Now  in  order  to  distinguish  the  two  output  voltages,  b^Ct)  is  shifted 

in  frequency;  it  is  fed  into  a  square-law  device  along  with  a  second  signal 

at  frequency  to  ,  The  frequency  to  is  chosen  to  be  much  larger  than  the  RMS 
2  2 

spectral  width,  A  to,  of  the  pulse  P(t),  while  w  is  very  much  smaller.  For 

co 

example,  if  the  operating  frequency,  —  ,  of  the  radar  is  10  Gc,  and  the 

A  to 

pulses  are  one  mi ci osecond  in  length  then  =  1  Me,  and  we  could  let 

oj  to  ** 

iol  =  .lKc,  and  — 2  =  1  Gc.  A  filter  is  located  at  the  frequency  shifter's 
22T  2  71 

output  which  passes  only  the  signal  at  the  difference  frequency,  to 
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This  can  be  written  as 


B^t)  = 


P(t-tv)  [a, 
r  j  1 


(v)  + 


A2  (v)e 


n(t-v 


1 


Bi (v)e 


j[w  (t-t.)- 
6  r 


V1 


(C.  7) 


This  voltage,  along  with  Bg(t),  is  fed  into  a  mixer  which  forms  the  square 
of  their  sum,  i.e.. 


■  Re(vt>  ♦  Y«)J 


[  Re  [  \ 


JW,  (t-tp-l  j[w  (t-t-)-W  t] 

(v)  +  Ag  (v)e  j  B1(v)e  °  r  * 


) 

J 


(C.8) 


jco  (t-tr)! 

+  Re  ^P(t-t-;)  [A1(v)  +  A2(v)e  1  J  B2(v)e  ° 


[Ai 


-i  2 


If  an  RF  filter,  whose  pussband  is  centered  at  frequency  u>  , 

£m 

at  the  mixer^output  it  will  pass  the  following  RF  signal 


is  located 
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C(V)  =  |P(t-t-’)|2  I  A  ( v )  1  2  +  |  A  (v)|2 
r  1  2 


+  A  (v)A  (v)e  +  A^(v)  A^(v)e  JB^  (v)B^  (v)  e  .  (C,9) 


If  this  signal  is  fed  into  a  synchronous  (phase)  detector  along  with  a 

reference  voltage  **t  frequency  w  +  to  the  detector  output  can  be 

2  1 

written  as 


?<t)  = 


|p(t-tr>r 


C2  (v)e 


,J2“i  (“  T-)j 


(C.10) 


where 


C  (v)  =  A  (v)  B  (v)  A* (v)  B* (v) 
o  112  2 

C](V)  =  ||A1(V)|2  +  iA2(v)'2j  B1(V)  B2(V) 
C2(V)  =  A1(V)  B1(V)  A2(V)  B2(V) 


(6. 11) 


As  mentioned  earlier^  this  result  is  valid  for  0  ,<  t  <  The 

output  for  any  value  of  t  can  be  written 
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N-l  r  jw  t  ju)  t 

"N(t)  =  S  I P(t-t£-fc  TQ)|2  lC0(v)e  1  ^  C!(v)  1 

k-o  L 

t 

j2co  (t~  ~  )“|  (C,  12) 

+  C2(v)e  2  J 


In  practice,  a  linear  saw-tooth  voltage  wave,  synchronized  with  the 

pulse  repetition  rate,  is  often  used  to  produce  a  visual  display  of  this 

output  on  an  oscilloscope.  Thus,  the  abscissa  nf  the  display  is  the 

range  coordinate  s  and  the  vertical  defleetion  is  proportional  to  the 

strength  of  the  target  return.  For  g' given  range  s  this  deflection  is 

,  s  2s 

given  by  the  values  of  r„(t)  sampled  at  t  =  k  t  —  .  The  average 

n  o  c 

value  of  this  deflection  after  N-l  repetitions  of  the  pulse  is 


Vs> 


1 

N 


P  (s-r) 
o 


N-l 

2 


t 

r 


+  C  (v)e 


J  /2s 
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j2w 


C2(v)e 


-f  k  t  - 
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r) 


'I 

► 

J 
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f  -jco 

P  (s-r)  <C  (v)  e 
o  ~ 


1° 


A  sin(N  w  Tb  |  .  /n-1  2s  \ 

i*r  „  , .  I  ‘-s ri  J“il—  v  - 

+  c  (v ) - : — — V  e  '  > 


v  *  '  *  »  \  1 
N  sin^o^  1  o j 


f 

sin 

+  S(v)  ra/Vo) 


V»)  [(B-1)To  *  iHr£] 


(C.13) 


where,  in  changing  variables  from  t  to  s,  we  have  redefined  the  pulse 
function  as 


Po(s)  -  !  P  f~j|2  .  (C.  14) 

It  will  be  shown  later  that  it  is  to  our  advantage  to  require  that 

caT  be  small,  say  0.5  radians  or  less.  If  this  is  the  case, 

1  A  7 

sin  w  t  —  co  t  and  the  absolute  value  of  the  second  and  third  terms  of 
1  o  1  o  } 

Equation  (C.13)  can  be  written  as 


P  (s-r) 
o 


sin|N  u  o I  r/.T  ,y 

i  *-?/  j“i  t»+ 

i',V^  f  T  \e  L\  2  / 

N  sin  [u  j 

\  X2  J 


2s 

c 


(C. 15) 


+  c2(v)  r-  c 

1  O 


r««-i>T  * 45  4  2i 

P  (s-r) 

2icfv)!  |c_(v)| 

o  C  J 

-  °  < 

1  2  ) 

L 

^  CJ  T 

!  N  4  N 

1  o  | 

L  J 
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Note  that  as  N  increases,  the  absolute  value  at  these  terms  approaches  zero 
and  in  the  limit  Equation  (C.13)  becomes 


lim 
N  — * 


jw  t 

R  (s)  =  P  (s-r)  C  (v)e  r 
N  o  o 


(C, 16) 


or 


jco  t 

R(s)  =  P  (s-r)  A  (v)  B  (v)  A  (v)  B*(v)  e  1  r 

O  1  i  Ci  £i 


(C 1 17 ) 


A  not  uncommon  case  is  when  all  the  patterns  are  reai^  and 
one  can  write  the  actual  real  system  output  as 


R(s) 


A 


R 

"1 


(V) 


A  <v) 
2' 


n 

2 


(v) 


(C.  18) 


But  if  2w^r/c~  71/2  the  output  is  zero,  a  rather  undesirable  result.  The 
cos  (2w^r/c  )  factor  is  due  to  the  time  delay  2r/c  between  the  transmitted 
pulse  and  that  which  is  received*  As  a  consequence  the  reference  signal 
in  the  synchronous  detector  leads  the  desired  signal  by  2^r/c  radians. 
ThiSj  in  general^  will  cause  a  reduction  in  the  system  output  and  for  real 
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patterns  it  is  given  by  cos  (2u>^r/c).  New  if  we  select  that 

radians,  then  the  reduction  in  the  system  output  will  be  to^cos \ ( . 5)^  100 

=  88$  of  the  maximum  possible  output.  This  will  occur  for  targets  at 

the  range  limit  r  *  For  targets  which  are  closer,  the  reduction  is  less* 
max  7 

A  better  and  still  rather  simple  scheme  would  be  to  set  the  reference 

signal  with  a  delay  so  that  it  is  exactly  in  phase  with  the  response  from 

targets  in  the  center  of  the  range  at  r  /2.  Then  it  would  lead  the 

max 

responses  from  ancl  lag  the  responses  from  r=0, by  0.25  rmdians,  The 

reduction  in  the  output  in  these  cases  would  be  to  96*8$  of  the  maximum 
possible  value.  Thus  we  see  that  by  keeping  the  phase  modulation 
frequency  low  enough  we  can,  to  a  reasonable  approximation,  obtain  the 
following  average  response  at  the  position  s  on  the  oscil loseope' s  range 
coordinate 


R(s)  =  P  (s-r)  A  (v)  A  (v)  B  (v)  B0(v) 
o  12  12 


(C. 19) 


However  it  can  be  shown  that  by  introducing  a  continuously  varying 

phase  compensation  in  the  reference  voltage  of  the  synchronous  detector 

one  can  completely  remove  the  phase  factor  2u>^r/c  for  all  values  of  r.  It  can 

also  be  shown  that  in  this  case  w  t  need  not  be  restricted  to  small  values 

1  o 

and  there  is  a  certain  value  of  OJ-^T0»  namely  4  Zf/3  ,  which  will  cause  the 
other  two  terms  of  the  detector  output  to  average  cut  to  zero  at  a  maximum 
rate.  This  scheme,  however,  is  rather  complicated  and  at  might  have  nn 
significant  advantage  over  the  one  mentioned  above  in  many  practical  situ¬ 


ations  * 
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For  this  simpler  system  it  is  clear  that  w  r  and  r  are  not 
K  J  1*  o'  max 

independent.  In  usual  practice  one  is  given  the  value  of  r 

150  kilometers.  From  this,  one  can  use  the  delay  formulas  to  obtain 


T 

O 


2r 

max 

■  '  • 
c 


(C. 20) 


“3 

In  this  case  t  =  2 (150)/300, 000  =  10  seconds.  The  pulses  are  repeated 
o  ? 

every  millisecond.  Finally,  we  require  that  C01To  23  .5;  hence  in  this  case 
-3 

^l"  ~  500  radians  per  second  which  is  approximately  80  cps. 

If,  as  is  commonly  .  done,  the  composite  patterns  of  the  antenna 
system  are  scanned,  then  for  a  scan  direction  u,  we  can  write  the  system 
output,  as  a  function  of  u  and  s,  as  follows 


R(s, u)  =  P  (s-r)  A. (v-u)  A* (v-u)  B  (v-u)  B*(V-ti)  .  (C.21) 

9  o  1  2  1  2 


This  is  the  system  response  to  a  point  target  in  the  v  direction  and  at 
range  rn 


SUPPRESSION  OF  BACKGROUND  NOISE 


Let  us  suppose  that  in  addition  to  the  distribution  of  passive  targets 
there  is  a  background  dis trihut i-  "-  of  inder»r,d*-nt',  n-.ioo  "urcco. 
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If  we  assume  that  the  antenna  has  an  RF  filter  which  passes  only  a  narrow 

band  of  frequencies  centered  at  to  we  can  give  the  following  description 

o 

of  the  noise  as  a  function  of  time  and  direction 

jw  t 

n(tjU)  =  N(tjU>e  °  (C.22) 

where  N(t^u)  is  the  narrow  band  complex  noise  phasor  at  time  t  due  to  the 
noise  source  in  the  u  direction.  It  is  a  random  function  of  time  and  we 
assume  it  is  stationary;  its  statistical  properties  are  invariant  under 
a  shift  of  the  time  origxn. 

The  system  output  will  contain  components  due  to  the  target 
reflections  alone  (signals)^  the  noise  alone;  and  cross-products  of  the 
reflections  and  the  noise.  But  since  the  noise  is  random  and  the  reflections 
are  not,  these  cross-products  terms  will  average  out  to  zero.  It  remains 
to  investigate  the  output  components  due  to  the  noise  alone.  The  two 
output  noise  voltages  arc 
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where 


(t) 


N(tj u)  (u)du  } 


(C. 25) 


and 


(t)  =  N  (t)eJ  ° 
B2 


(C, 26) 


where 


(t)  = 


CO 

/  N(t,u)  B 
-  oo 


(u)du  • 


(C. 27) 


The  mixer  forms  the  product  of  the  sunt  of  these  two  voltages  and  the  filter 
at  the  mixer  output  passes  the  component  at  frequency  to  It  is  given  by 


\(t)  =5  [NBl(t)  NB2(t) 


t 


(C, 28  ) 


The  output  of  the  synchronous  detector  corresponding  to  this  voltage  is 


pi opoi tiunal  to 
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n  (t) 
o 


N„  (t) 


* 


-  j^,  t 

(t)  e 


(C.  29) 


The  average  value  of  this  output,  sampled  at  t  =  2s/c  +  k  T 
for  N  samples,  is  given  by 


n  ,.(s) 
o,  N 


.  i  "s  BB  (21  ♦  „  T  V*  (2a  *  k  O  +  k  T=)  . 

K  k=o  B1  °1  %  1°  ’ 


(C,  30) 


This  is  proportional  to  the  oscilloscope1 s  average  vertical  deflection  at 

the  position  s  of  the  horizontal  range  coordinate.  Since  the  noise  is 

stationary,  and  essentially  uncorrelated  between  pulses,  we  can  write  the 

expected  value  of  n  __(s)  as 

°,N 


E  < 


no,  K 


1 

N 


N 


B1B2 

N 


2s 

c 


N-l  -jfro  t 
^  1  o 

e 

k=:0 


=  N 


B  B 
1  2 


sin  N  ca>  t 
1  o 


N  sin  t 
1  o 


e 


-n 


N-l 

2 


1 


(C. 31 ) 
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where 


N. 


B  B 
1  2 


(C. 32) 


is  the  mathematical  expectation  or  statistical  average  of  N  (—  +  k  T  | 

\  B,\c  i 

*  f  2b  \  1 

nb  ^  T|  f  and  due  to  stationarity  is  independent  of  k.  From 


Equation  (C.3I)  we  can  write 


E  4 

n  M(s), 
I  o.  N  [ 

|nb,b  1 

sin^N 

*  > 

N  sin 

i 

L  '  J 

1  2 

(C. 33) 


and  except  when  w  t  =  2qff,  for  integral  q,  we  have 


lim 

r  i 

sinl  N  to  t  ) 

*tVJ| 

<< 

li» 

1  \v‘ 

V  i  / 

N — >oo 

< 

i-N  sin 

/ 


(C. 34) 


=  0  • 
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APPENDIX  D 

RADAR  MAPPING  OF  A  TARGET  DISTRIBUTION 

If  the  output  voltage  B.^ .  r.  u,  v);  as  given  by  Equation  (207)  of 
Chapter  10  is  first  shifted  in  frequency  to  ^  ,  and  then  mixed  with 

B^(t u;v)  as  given  by  Equation  (208)^  the  output  of  the  band-pass 
filter  (at  frequency  )  can  be  written  as 


m(t,  t ,  r,  u,  v)  = 


0O  oo 

lb'1-'- 


u)  B1(u1-u)  Aj(u2-v)  B2(u2-v) 


CO  fto 

// 


N-l  N-l 

vr* 

£-t  <L-i 

k— o  i"0 


P(t-  -(i-r,  )-2kT  ) 
cl  o 


j 


P*(t~  — (r~r  )-2iT  )e 
c  2  o 


* 

g  (t+ 


r"r2 
c  ’ 


-  j  (w1+w2>  t 


U2»r2)dUldU2dridr2e 


(D 1 1 ) 


Note  that  the  usual  correlation  delay  t,  has  not  been  introduced.  Con¬ 
sequently  the  system  output  will  be  the  cross-correlation  of  the  two 
antenna  voltages  evaluated  at  t  =  oo  We  have  also  assumed  that  the  RF 
filters  of  the  system  are  capable  of  passing  the  video  pulses  of  the 
radar  with  no  distortion.  However  even  with  these  simplifications  the 
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mathematical  analysis  is  still  quite  involved. 
Letting  the  reference  voltage  be 


VREF(t)  =  6 


r  *  i  [ 

-J[(Wl4W2)t  +~[T- 


r 

max  | 
2 


(D.2> 


the  output  of  the  synchronous  detector  becomes 


RN(t,i,r,  u,v)  = 


oO  oo 

IP 

-00  -CO 


v  N-l  N-l 

1<U-U1)  V'"V  „S 

k--o  i-o 


GO  CO 

II 

-oo  -OO 


P  ^t+  -  (f-r  )-2kr 

c  1  of 


(D.3) 


2<jo 

*12  l  J~T  [i_r]  <“ri  *  r"r? 

P  <t+  -  (r-r  )-2ir  >e  "  y~- 

c  2  of 


gftn — gr— ,  u2r2)dr1-dr2du1du2 


where  we  have  defined 


(u)  =  (u)  B^  (u) 


(D.4) 
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C9  (v)  =  A^Cv)  (v)  - 


(D.5) 


There  are  other  terms  in  the  output,  just  as  in  the  single  target  case 
(see  Equation  C.10),  but  since  they  all  have  a  zero  time-average  we  have, 
for  simplicity  of  presentation,  omitted  them  in  the  above  output 
expression.  Since  the  pulses  are  essentially  non-overlapping  i  (orthogonal) 
we  can  write  the  output  as 


-  V  V 

RN(sji,r?u  v)  =  C^CuKgCv)  * 


CO  oO 

2  I  I  p  fs+i  "r  -2kr  je 

k=o  J  J  v  1  max' 

o  o 


.  2C0  *  .  ' 

j—  («+#) 
c 


(D.  6) 


P  ( s+r-r-2kr  ]e  c 

^  2  max 


i 


2co  "  >  ■' 

j  — £  (s+r) 


g 


(fitili'  u>ri)"*  ( 


2s+r-r2J  y,r2Ur1dr-2 


where  s  =  ct/2„ 

For  a  target  at  a  given  range,  for  example  r  ,  we  will  assume  that 
its  reflection  coefficient,  g(t,u,  r  ),  is  slowly  varying  compared  to  the 
pulse  that  strikes  ±t„  This  occurs  when 
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s  =  i-r  -2kr  =  0 
1  max 


(D.7) 


or 


r,  =  s+f  -  2kr  .  D,8) 

1  max 


Similarly  the  pulse  strikes  the  relatively  unchanging  target  at  range  r 
when 


r  =  s  +  r  -  2kr  (0.9) 

2  max  • 


Substituting  these  results  into  Equation  (D.6)  gives 


_  \J  V* 

RN(s,i,r,  u,v)=  Cju)  C2(v) 


00  00 

.*2  ff 

-  J  J 

o  o 


•  * 

P  I  ~2kr  1 

J  V 


Pfs+r-r  -2kr  j  e  C 
1  2  mar 


2oj 

j  —2  (i_r) 


t  max 


(D.10) 


J 


g( s+2kr  ,  u,  r.  s+2kr  v  r  )dr  dr 

_ max’  *  1  I  I _ max-*  3  2/  1  2 
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1 


A<u) 


C2<V) 


N-l 

2 

k=o 


P 


fs+i  -2kr  ]  p+|j 
^  max /  \ 


s+r«2kr 


max 


g f s+2kr  . u, s+i -2kr 

max7  7 


max* 


*  )  s  h 

max  /  ( 


s+2kr  ,vfs+r-2kr 
max7  3  i 


max 


(D.ll) 


Letting  the  output  be  sampled  at  s  =  2kr  (i.e,, 

max  3 


average  output  is  given  by 


at  t  =s  2kT  )  the 
o  7 


.  2w  , 

N-l  J — «a-r)j 
£  e  C  <C  (a) 

k-o  [ 


v/*  —  —  * 

C2(v)  P(f)  P  (r)  * 


(D ,  ]  2) 


Now  if  the  target  function  g(t^u^r)  satisfies  the  ergodic  hypothesis, 


then  we  have 
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E*g(t,u,  O  g  (t-r,v,r) 


>  =  G  (T,ijr,U,v) 


(D.13) 


where  G(T^l^r,u, v)  is  the  mutual  coherence  function  of  the  target 
distribution. 

The  expected  value  of  the  system  output  after  F  pulses  is 


1 

N 


N 

2 

k=o 


(u)  C*(v)  P(i)  P*(r) 


*  G  (o,f,  r,  u,v)^l 


(D.14) 


or 


2u> 


R(f>rJuJv)  = 


5 
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(l-  r'i 


—  .  .  — * 


E  <u)  C*(v)  PCO  p”(r)  *  G(o^  t,  r}utv) 


(D.15) 


'2w  ' 

-J  (t-T) 


e 


Multiplying  both  sides  of  the  above  equation  by- 


gives 
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K(i,r,u,v)  *  e  c  R  (i>r#u>v) 


(D, 16) 


a  C^fu)  C ^(v)  P(£)  P*(r)  *  G  (ofi}rt  u,v) 


(D .  17  ) 


as  the  required  system  output 
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Attn:  Dodge  Library 
Boston  15-  Massachusetts 

North  American  Aviation,  Inc. 

Attn°  Technical  Library  (M/F 
Engineering  Dept.) 

4300  E  Fifth  Avenue 
Columbus  16  Ohio 

North  American  Aviation,  Inc, 

Attn  Technical  Library 
(M/F  Dept.  56) 

International  Airport 
Los  Angeles,  California 


Northrop  Corporation 
NORAIR  Division 

1001  East  Broadway 

Attn:  Technical  Information  (3924-3) 
Hawthorne,  California 

Ohio  State  University  Research 
f  ouiiua.  t  ion 

Attn:  Technical  Library 

(M/F  Antenna  Laboratory) 

1314  Kinnear  Road 
Columbus  12 f  Ohio 

Ohio  State  University  Research 
Foundation 

Attn:  Dr.  C-  H-  Walter 
1314  Kinnear  Fcoad 
Columbus  12,  Ohio 

Phil co  Corporation 

Government  &  Industrial  Division 

Attn:  Technical  Library 

(M/F  Antenna  Section) 

4700  Wissachickon  Avenue 
Philadelphia  44,  Pennsylvania 

Westinghouse  Electric  Corporation 

A i  r  A rms  Division 

Attn:  Librarian  (Antenna  Lab) 

P,  0-  Box  746 
Baltimore  3,  Maryland 

Wheeler  Laboratories 

Attn:  Librarian  (Antenna  Lab) 

Box  561 

Smithtown,  New  York 

Electrical  Engineering  Research 
Laboratory 
University  oi  Texas 
Box  8026,  University  Station 
Austin,  Texas 

University  of  Michigan  Research 
Institute 

Electronic  Defense  Group 
Attn:  Dr,  J,  A,  M  Lyons 
Ann  Arbor  Michigan 

Radio  Corporation  of  America 
RCA  Laboratories  Division 
Attn  Technical  Library 

(M/F  Ante nn a  Sc ct iGa) 
Princeton,  ^ew  Jersey 
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Radiation,  Inc, 

Attn:  Technical  Library  (M/F) 
Antenna  Section 
Drawer  37 

Melbourne,  Florida 

Ramo-Wooldridge  Corporation 
Attn-  Librarian  (Antenna  Lab) 
Canoga  Park,  California 

Rand  Corporation 

Attn:  Librarian  (Antenna  Lab) 

1700  Main  Street 

Santa  Monica,  California 

Rantec  Corporation 

Attn:  Librarian  (Antenna  Lab) 

23999  Ventura  Blvd . 

Cal abas as ,  California 

Raytheon  Corporation 
Equipment  Division 
Library  -  J.  Portsch 
P«  0~  Box  520 
Waltham  54,  Massachusetts 

Republic  Aviation  Corporation 
Applied  Research  &  Development 
Division 

Attn:  Librarian  (Antenna  Lab) 
Farmingdale,  New  York 

Sanders  Associates 
Attn:  Librarian  (Antenna  Lab) 
95  Canal  Street 
Nashua,  New  Hampshire 

Southwest  Research  Institute 
Attn*  Librarian  (Antenna  Lab) 
8500  Culebra  Road 
San  Antonio .  Texas 


Sperrj'  Gyroscope  Company 
Attn.  Librarian  (Antenna  Lab) 
Great  Neck,  L »  I » ,  N ew  York 

Stanford  Electronic  Laboratory 
Attn:  Librarian  (Anterma  Lab) 
Stanford,  Calnornia 

Stanford  Research  Institute 
Attn:  Librarian  (Antenna  Lab) 
Menlo  Park,  California 

Sylvania  Electronic  System 
Attn:  Librarian  (M/F  Antenna  & 
Microwave  Lab) 

100  First  Street 
Waltham  54,  Massachusetts 

Sylvania  Electronic  System 
Attn:  Librarian  (Antenna  Lab) 

P .  0 .  Box  188 

Mountain  View,  California 

Technical  Research  Group 
Attn:  Librarian  (Antenna  Section) 
2  Aerial  Way 
Syosset,  New  York 

Ling  Temco  Aircraft  Corporation 
Temco  Aircraft  Division 
Attn:  Librarian  (Antenna  Lab) 
Garland,  Texas 

Texas  Instruments,  Inc. 

Attn:  Librarian  (Antenna  Lab) 

6000  Lemmon  Avenue 
Dallas  9,  Texas 

A,  S*  Thomas,  Incr 

Attn:  Librarian  (Antenna  Lab) 

355  Providence  Highway 
Westwood,  Massachusetts 


H  R  B  Singer  Corporation 

Attn  Librarian  (Antenna  Lab)  New  Mexico  State  University 

State  College.  Pennsylvania  Head  Antenna  Department 

Physical  Science  Laboratory 
Sperry  Microwave  Electronics  Company  University  Park,  New  Mexico 

Attn  Liorarian  (Antenna  Lab) 

P  O  Box  1828 
Clearwater^  Florida 
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Bell  Telephone  Laboratories,  Inc, 
Wliippany  New  Jersey 
At  +  n*  Technical  Reports  Librarian 
Room  2A-165 

Robert  G  Hansen 
Aerospace  Corporation 
Box  95085 

Los  Angeles  45;  California 

Dr,  Pi  chard  C-  Becker 
10829  Berkshire 
Westchester^  Illinois 

Dr  W  .  M  .  Hall 
Raytheon  Company 
Surface  Radar  and  Navigation 
Operations 
Bus  lull  ?um  Ro^d 
Wayland ,  Massachusetts 

Dr~  Robert  I  Carrel 
Collins  Radio  Corporation 
Antenna  Section 
Dallas^  Texas 

Dr  A  K  Chatterjee 
Vice  Principal  and  Head  of  the  De¬ 
partment  of  Research 
Birla  Institute  of  Technology 
P  0  Mesra 

Di st ri ct -Ranchi  (Bihar)  India 

University  of  Dayton 

Research  Institute 

Attn  Professor  Douglas  Hanneman 

300  College  Park 

Dayi on ,  Ohio 

Techn'j  srhe  Hochschule 
Attn  H  H  Meinke 
Munich  Germany 

NASA  Goddard  Space  Flight  Center 
Attn:  Antenna  Branch 
Mr  Lantz 

G  r  e  e  n  be  1 1  M  a  r  y  1  and 

Aeronautic  Division 

Ford  Motor  Comp  an  y 

Ford  Poid  -  Attn  Mr  J  M  Black 

Newport  Beach  California 


Professor  A«  A,  Oliver 
Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,  New  York 

Uo  Se  Naval  Ordnance  Laboratory 
Attn:  Technical  Library 
Corona,  California 

Avco  Corporation 

Research  and  Advanced  Development 
Division 

Attn  Research  Library/T »A „  Rupprecht 
201  Lowell  Street 
Wilmington,  Massachusetts 

Raytheon  Company 
Missile  and  Space  Division 
Attn:  Research  Library 
Bedford,  Massachusetts 

Teledyne  Systems,  Incorporated 
Attn:  Technical  Library/Antenna  Section 
1625  F  .  126th  Street 
Hawthorne,  California 

National  Research  Council 
Attn*  Microwave  Section 
Ottawa  2:  Canada 

Sicbak  Associates 
Attn:  W  Sichak 
518  Franklin  Avenue 
Nut  ley,  New  Jersey 

W,  T.  Patton 
2208  New  Albany  Road 
Cinn  Township 
Riverton  Post  Office 
New  Jersey 

Radio  Corporation  cf  America 
Missile  and  Service  Radar  Division 
Attn:  Manager,  Antenna  Engineering 
Skill  Center 
Moo res town.  New  Jersey 
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Advanced  Development  Laborator  los ,  Inc. 

Haines  Street 

Nashua  New  Hampshire 

Cornell  Aeronautical  Lab 
Attn  Research  Library 
Buffalo,  New  York 

Fairchild  Stratos  Corporation 
Aircraf t -Missi les  Division 
Attn*  Technical  Library  (Antenna  Section) 
Hagerstown  10 }  Maryland 

General  Electric  Company 
Light  Military  Electronics  Department 
French  Road 
Utica  New  York 

General  Electric  Company 
General  Engineering  Laboratory 
Bldg  371,  Room  478 
Schenectady,  New  YoYk 

Goodyear  Aircraft  Corporation 
Antenna  Department 
Litchfield  Park 
Phoenix.  Arizona 

Laboratory  for  Electronics,  Inc. 

Antenna  Department 
1073  Commonwealth  Avenue 
Boston  15,  Massachusetts 

Lockheed  Aircraft 

Electronic  and  Armaments  System  Office 
Burbank  C  alif ornia 

Mo^orol^i  7  nc 

Western  Military  Electronics  Division 
P  0  Box  .41^ 

Scottsdale,  Arizona 

Phiico  Corporation 
Ante" a  Sec*  ion 
38-5  Fabi  way 
Palo  Al*o  California 


Avco  Corporation 

Electronic  and  Ordnance  Division 
Attn;  Technical  Library 
Cincinnati  41,  Ohio 

Bendix  Corporation 
Research  Labs  Division 
Attn:  G.  M.  Peace 
Southfield  (Detroit),  Michigan 

Douglas  Aircraft  Co,,  Inc. 

Attn:  MSSD  Technical  Library 
(Antenna  Section) 

3000  Ocean  Park  Blvd, 

Santa  Monica,  California 

Emerson  and  Cuming,  Inc. 

Attn:  E.  J.  Luoma 
869  Washington  Street 
Canton,  Massachusetts 

Radioplane  Company 

Attn:  Librarian  (M/F  Aerospace  Lab) 
in  15  Kancho  Conelo  Blvd. 

Newbury  Park.  California 

Hughes  Aircraft  Company 
Attn:  Antenna  Section 
Full  rton,  California 


Spice  rec~  ology  Laboratory 
At +  ”  "tsea^ch  Library 
P  0  vex  95100 
Los  Angeles  45,  California 


ANTENNA  LABORATORY 


TECHNICAL  REPORTS  AND  MEMORANDA  ISSUED 


Contract  AF33 (616^-310 

^Synthesis  of  Aperture  Antennas/'  Technical  Report  No.  I,  C,  T.  A.  Johnk, 
October,.  1954.  *  AD-52134 

nA  Synthesis  Method  ior  Broadband  Antenna  Impedance  Matching  Networks, 
Technical  Report  No,  2J  Nicholas  Yaru,  i  February  1955=  *  AD  61049. 

f,The  Assymmetncal  ly  Excited  Spherical  Antenna/'  Technical  Report  No.  3, 

Robert  C.  Hansen,  30  April  1955:  *  AD-66220. 

"Analysis  of  an  Airborne  Homing  System,”  Technical  Report  No.  4,  Paul  E. 

Mayes,  i  June  1955  (CONFIDENTIAL).  AD-80147. 

Coupling  of  Antenna  Elements  to  a  Circular  Surface  Waveguide,"  Technical 
Report  No.  5,  H.  F„  King  and  R *  H.  DuHamel ,  30  June  1955.  *  AD-75733. 

"Input  Impedance  of  a  Spherical  Ferrite  Ar.tenna  with  A  Latitudinal  Current/1 
‘technical  Report  No  6?  W.  L.  Weeks,  20  August  1955  AD- 75732. 

"Axially  Excited  Surface  Wave  Antennas, "  Technical  Report  No,  7,  D.  E,  Royal, 
10  October  1955  *  AD-75731. 

Homing  Antennas  for  the  F-86F  Aircraft  (450-2500  me),  Technical  Report  No,  8, 
P.  E.  Mayes,  R.  F  Hyneman,  and  R  C.  Becker,  20  February  1957, 

Ground  Screen  Pattern  Range,”  Technical  Memorandum  No,  1,  Roger  R.  Trapp, 

10  July  1955.  * 


Contract  AF33 (616) -3220 


'Effective  Permeability  of  Spheriodal  Shells,"  Technical  Report  No,  9, 

E,  J.  Scott  and  R.  H  DuHamel,  16  April  1956.  AD-98615. 

"An  Analytical  Study  of  Spaced  Loop  ADF  Antenna  Systems,"  Technical  Report 
No.  10r  D.  G.  Berry  and  J.  B,  Kreer,  10  May  1956. 

"A  Technique  for  Controlling  the  Radiation  from  Dielectric  Rod  Waveguides, 
Technical  Report  No.  11,  J.  W ,  Duncan  and  R.  H.  DuHamel,  15  July  1956,  * 
AD-107344 

'Direction  Characteristics  of  a  U-Shaped  Slot  Antenna,"  Technical  Report 
No  12.  Richard  C.  Becker,  30  September  1956.  **  AD-116104. 

‘Impedance  of  Feirite  Loop  Antennas,"  Technical  Report  No,  13.  V.  H.  Rumsey 
and  Wt  L  Weeks 4  15  October  1956  AD  119780 
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Closely  Spaced  Transverse  Slots  in  Rectangular  Waveguide/'  Technical  Report 
Mo.  I4f  Richard  F  Hyneman,  20  December  1956  AO-127455. 

Distributed  Coupling  *o  Surface  Wave  Antennas/'  Technical  Report  No.  15, 

R.  R  Hodges,  Jr.,  5  January  1957.  AD-127454, 

"The  Characteristic  Impedance  of  the  Fin  Antenna  of  Infinite  Length, 1 
Technical  Report  16.  Robert  L  Carrel,  15  January  1957.  *  AD~1274530 

On  the  Estimation  of  Ferrite  Loop  Antenna  Impedance/'  Technical  Report  No,  17, 
Walter  L,  Weeks  10  April  1957,  *  AD-143989. 

A  Nute  Concerning  a  Mechanical  Scanning  System  fur  a  Flush  Mounted  Line 
Source  Antenna/  Technical  Report  No„  18,  Walter  L.  Weeks,  20  April  1957. 
AD-143990 

Broadband  Logarithmically  Periodic  Antenna  Structures/  Technical  Report  No,  19, 
R  H  DuHamel  and  D  F.  Isbell,  1  May  1957 .  AD-140734 „ 

’Frequency  Independent  Antennas/  Technical  Report-  No,  20,  Vfc  H0  Rumsey, 

25  October  1957.  AD-153260,  '  ” 

"The  Equiangular  Spiral  Antenna/  Technical  Report  No.  21.  J.  D0  Dyson, 

15  September  1957.  AD-145019. 

7  Experimental  Investigation  of  the  Conical  Spiral  Antenna,  '  Technical  Report 
No _ 22,  E,  L  Cartel,  25  May  1957  +  *  AD-144021. 

Coupling  De^een  a  Parallel  Waveguide  and  a  Surface  Waveguide/’  Technical 
Report  No  E  .i  3 coir,  10  August  195/.  AD- 

Launching  Efficiency  of  Wires  and  Slots  for  a  Dielectric  Rod  Waveguide/ 

Technical  Report  No. _ 24.  Jm  W  Duncan  and  R  H  DuHamel,  August  1957, 

AD-1 41030 . 

The  Char acter i s r ic  Impedance  of  an  Infinite  Bicnnical  Antenna  of  Arbitrary 
Cross  Section/  Technical  Report  Vo,  25,  Robert  L  Carrel,  August  1957. 
AD-144121  ~ 

Cav i ty -Backed  Slot  Antennas/  Technical  Report  No,  26.  R.  J6  Tector^  30 
October  1957,  AD-149451 

'Coupled  Waveguide  Exci'ation  of  Traveling  Wave  Slot  Antennas,  Technical 
Report  No  27,  W  L  Weeks.  1  December  1957  AD-155178 

’Phase  Velocities  in  Rectangular  Waveguide  Partially  Filled  with  Dielectric/' 
Technical  Report  No,  28.  W  L  Weeks,  20  December  1957  AD-155177, 


Measuring  rhe  Capa,:rance  per  Unit  Length  of  Bi coni cal  S+ructuies  of  Arbi¬ 
trary  Cross  S t c r ion  Technical  Report  No  29,  J  D  Dyson  10  January  1958 
AD  153561  ^  ~  - 


Non-Planar  Logarithmically  Periodic  Anrenna  S-rucTure,.  Technical  Report 
No  30  r  D  F  Isbell  20  February  1358  AD-156203 
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Electromagnetic  Fields  in  Rectangular  Slots/'  Technical  Report  No,  31, 

N.  J.  Kahn  and  P  E .  Mast.  10  March  1958,  AD-156709. 

"The  Efficiency  of  Excitation  of  a  Surface  Wave  on  a  Dielectric  Cylinder," 
Technical  Report  No.  32,  J,  W„  Duncan.  25  May  1958.  AD-162052. 

"A  Unidirectional  Equiangular  Spiral  Antenna/1  Technical  Report  No.  33, 

J.  D(  Dyson,  10  July  1958,  AD-201138. 

‘Dielectric  Coated  Spheroidal  Radiators,"  Technical  Report  No.  34,  W.  L . 
Weeks^  12  September  1958,  AD-204547. 

"A  Theoretical  Study  of  the  Equiangular  Spiral  Antenna/'  Technical  Report 
No.  35,  P.  E.  Mast,  12  September  1958,  AD-204548 


Cont  ract  AF33  (61 6 ) -6079 


Use  of  Coupled  Waveguides  in  a  Traveling  Wave  Scanning  Antenna/*  Technical 
Report  No.  36,  R.  MacPhie,  30  April  1959.  AD-215558. 

'On  the  Solution  of  a  Class  of  Wicncr-Hopf  Integral  Equations  in  Finite  and 
Infinite  Ranges/*  Technical  Report  No  37,  R„  Mittra,  15  May  1959. 

AD-220543 „  ~ 

"Prolate  Spheriodal  Wave  Functions  for  Electromagne tic  Theory/*  Technical 
Report  No.  38,  W„  L  Weeks,  5  June  1959.  AD-226727.  ~~ 

*Log  Periodic  Dipole  Arrays/’  Technical  Report.  No  39t  Dt  E,  Isbell^  1  June 
1959  AD-220651 

'A  Study  of  the  Coma-Corrected  Zoned  Mirror  by  Diffraction  Theory/*  Technical 
Report  No.  40,  S,  Dasgupta  and  Y.  T.  Lo,  17  July  1959.  AD-225273, 

1  The  Radiation  Pattern  of  a  Dipole  on  a  Finite  Dielectric  Sheet/'  Technical 
Report.  No,  41,  K.  G.  Balmain,  1  August  1959. 

"The  Finite  Range  Wiener-Hopf  Integral  Equation  and  a  Boundary  Value  Problem 
m  a  Waveguide/*  Technical  Report  No,  42,  Rs  Mittra,  1  October  1959. 
AD-22y279. 

*  Impedance  Properties  of  Complementary  Mul ti terminal  Planar  Structures/' 
Technical  Report  Np^  43,  G.  A„  Deschamps,  11  November  1959.  AD-230625. 

On  the  Synthesis  of  Strip  Sources,"  Technical  Report  No.  44,  R,  Mittra, 

4  December  1959  AD -230624 

''Numerical  Analysis  of  the  Eigenvalue  Problem  of  Waves  in  Cylindrical  Wave¬ 
guides  ,  cal  Rep ji  t  Nu.  45.  C„  H .  Tang  and  Y#  T,  Lo,  11  March  1960. 

AD-234536* 
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New  Circularly  Polarized  Frequency  Independent  Antennas  with  Conical  Beam 
or  Ommdi  rect  1  onal  Pa  *  terns/’  Technical  Heport  No,  4b,  J  JJ.  Hyson  and 
P.  E  Mayes,  20  June  1960.  AD-241 321. 

Logarithmically  Periodic  Resonant-V  Arrays/’  Technical  Report  No,  47 ;  P.  Et 
Mayes  and  R  L..  Carrel,  15  July  1960,  AD-246302. 

A  Study  of  Chromatic  Aberration  of  a  Coma -Cor rect ed  Zoned  Mirror  T  Techni cal 
Report  No  48,  Y.  T#  Lo,  June  1960.  Ad-252492.  ~~ 

Evaluation  of  Cr oss -Cor rel at l on  Methods  m  the  Utilization  of  Antenna  Systems/ 
Technical  Report  No  49,  R„  MacPhie,  25  January  1961.  AD-252493. 

Synthesis  of  Antenna  Products  Patterns  Obtained  from  a  Single  Array,”  Technical 
Report  No,  £0,  R,  H  MacPhie,  25  January  1961 „  AD-252494.  ~~  “ 

On  the  Solution  of  a  Class  of  Dual  Integral  Equations,  Technical  Report  No,  51_, 
R„  Mittva,  1  October  1961.  AD-264557. 

Analysis  and  Design  of  the  Log-Periodic  Dipole  Antenna/1  Techni  cal  Report 
No«  .52,  Robert  Ln  Carrel,  I  October  1961.  AD-264558. 

A  Study  of  the  Non-Uniform  Convergence  of  the  Inverse  of  a  Doubly- Infinite 
Matrix  Associated  with  a  Boundary  Value  Problem  in  a  Waveguide, *  Technical 
Report  No.__53,  R_  Mittra,  i  October  1961.  AD-264556.  — ^ 


Contract  AF33 (616 ) -8460 

xhe  Coupling  and  Mutual  Impedance  Between  Balanced  Wire-Arm  Conical  Log-Spiral 
Antennas,”  Technical  Report  No.  54,  J,  D  Dyson,  June  1962.  AD-285678. 

An  Investigation  of  the  Near  Fields  on  the  Conical  Equiangular  Spiral  Antenna, 
Technical  Report  No.  55.  0.  L,  McClelland,  May  i962,  AD-276463. 

Input  Impedance  of  Some  Curved  Wire  Antennas.”  Technical  Report  No.  56. 

C.  H.  Tang,  June  1962..  ATI-276461.  ~ 

Polygonal  Spiral  Antennas,”  Technical  Report  No.  57,  C.  H  Tang,  O.  L. 
McClelland,  Tune  1962.  — -  — 


On  Increasing  the  Effective  Aperture  of  Antennas  by  Data  Processing,  ' 
Technical  Report  No.  58,  R  H  MnrPhip  Tnlv  i  pro. 

Theoretical  Study  of  a  Class  of  Logarithmically  Periodic  Circuits,” 
Technical  Report  No _ 59.  R  Mittra,  July  1962.  AD-282847. 

Backward  Wave  Radiation  from  Periodic  Structures  and  Application  to  the 
Design  of  Fr equency- Independent  Antennas/  Technical  Report  No.  60, 

P  E  Maves .  G  A  Deschampsj  and  W  T  Patton,  December  1962  AD-291765. 
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’’The  Backfire  Bifilar  Helical  Antenna,  '  Technical  Report  No.  61,  W.  T.  Patton, 
September  1962,  AD -28 908 4 

"On  the  Mapping  oy  a  Cross-Correlation  Antenna  System  of  Partially  Coherent 
Radio  Sources,"  Technical  Report  No,  62,  R.  H„  MacPhie,  Octobe  r  1962. 

AD-289085. 

'On  a  Conical  Quad-Spiral  Array,"  Technical  Report  No.  63,  J.  D.  Dyson, 
September  1962.  AD-288252. 

"Antenna  Impedance  Matching  by  Means  of  Active  Networks,"  Technical  Report 
No.  64,  S.  Laxpati.  R.  Mittra,  November  1962,  AD-291766. 

"On  Maximizing  the  Signal-To-Noise  Ratio  of  a  Linear  Receiving  Antenna  Array," 
Technical  Report  No.  65,  R.  H.  MacPhie,  January  1963. 

"Application  of  Cross-Correlation  Techniques  to  Linear  Antenna  Arrays", 

Technical  Report  No.  67,  R.  H,  MacPhie,  March  1963. 

"A  Study  of  Wave  Propagation  of  Helices",  Technical  Report  No.  68,  Paul  W.  Klock, 
March  1963. 

"Relative  Convergence  of  the  Solution  of  a  Doubly  Infinite  Set  of  Equations," 
Technical  Report  No.  69,  R.  Mittra,  April  1963. 

"Theoretical  Brillouin  (k~P)  Diagram  for  Monopole  and  Dipole  Arrays  and  Their 
Application  to  Log-Periodic  Antennas,"  Technical  Report  No.  70,  R.  Mittra  and 
K,  E,  Jones,  April  1963. 


*  Copies  available  for  a  three-week  loan  period. 


**  Copies  no  longer  available 


